Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,703,251 A→T 7.5% intergenic (+17/+17) add → / ← ydgJ adenosine deaminase/putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,703,2510AT 7.5% 166.1 / 2.7 62intergenic (+17/+17)add/ydgJadenosine deaminase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (23/34);  new base T (2/3);  total (25/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCATCACGAATAGCCG  >  NC_000913/1703103‑1703398
                                                                                                                                                    |                                                                                                                                                   
tataCCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAAGGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAACCGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCaa                                                                                                                                                     >  1:403035/1‑149 (MQ=255)
    ccGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATCGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGc                                                                                                                                                 <  2:84632/149‑1 (MQ=255)
     cGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATCGATGGTGCAAATGCa                                                                                                                                                >  2:507685/1‑149 (MQ=255)
       ttGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCAcc                                                                                                                                              <  1:1114521/149‑1 (MQ=255)
          ccGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCATTTGCAccatc                                                                                                                                           <  2:278636/149‑1 (MQ=255)
            gcgcCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAACAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCAccatcca                                                                                                                                         >  2:232935/1‑149 (MQ=255)
               ccAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCAttt                                                                                                                                      <  1:845224/149‑1 (MQ=255)
               ccAGCTGCTGGGTTATCCCGCGAGAAAATCCGACATGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCAttt                                                                                                                                      <  2:317817/149‑1 (MQ=255)
                 aGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCTGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCCAAGTAACAAAATGGTTGGTGCATTTGCACCATCCAttttg                                                                                                                                    <  1:256662/149‑2 (MQ=255)
                 aGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCAttttt                                                                                                                                    >  1:49571/1‑149 (MQ=255)
                 aGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCAttttt                                                                                                                                    <  2:156722/149‑1 (MQ=255)
                   ctgctgGGTTATCCCGCGAGCAAATCGCCCAGGCACATCATAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGCTGCAAATGCACCATCCATTTTTCa                                                                                                                                  <  1:251332/149‑1 (MQ=255)
                     gctgCGTTATCCCGCGAGCAAATCCGCCAGTCACAGTTTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATg                                                                                                                                <  1:808539/149‑1 (MQ=255)
                           ttATCCCGCTAGCAAATCCGCCAGCCACAGATTAATGGTCACGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCTCGAAGTAACAAAATGTATGGTGCAAATGCACCATCCATTTTTCATGCAAGGc                                                                                                                          <  1:369756/149‑1 (MQ=255)
                            tATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGca                                                                                                                         <  1:947083/149‑1 (MQ=255)
                              tCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCTCTGAGGAACAACGCGCACTGCGAGAAAACGTCGCCGCGAAGTAACAAAATGGATGGTGGAAATGCACCTTCCATTATTCAGGCAAGGcaca                                                                                                                       >  1:669885/1‑149 (MQ=255)
                                             cgacAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTTAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGtt                                                                                                        <  2:560527/146‑1 (MQ=255)
                                                cAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTgg                                                                                                     <  2:259849/149‑1 (MQ=255)
                                                 aGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGc                                                                                                    <  1:345986/149‑1 (MQ=255)
                                                 aGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGc                                                                                                    <  2:810930/149‑1 (MQ=255)
                                                       aGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTc                                                                                              >  2:88892/1‑149 (MQ=255)
                                                       aGATTAATGGGCTGGAAATGGCTTTCCTCCGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCCCCCTCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTc                                                                                              >  2:479676/1‑149 (MQ=255)
                                                          ttAATTGTCTGTAAATGGCTGGCCTCAGCGCTGAGGAAAAACGCGCACTGCTAGAAAAAGTCGCCTCGAAGTAACAAAATGGATGGTGAAAATGCACCATCCATTTTTCATGCAAGTCACAAAGTCGCGCGATGTTTGGCGGATTCGAt                                                                                           <  1:526615/149‑1 (MQ=255)
                                                              tGGTCTGGAAATGGCGTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGccc                                                                                       <  2:46984/149‑1 (MQ=255)
                                                               ggTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGcccc                                                                                      >  1:297865/1‑148 (MQ=255)
                                                                  cTGGAAATGGCTTTCCTCATCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCt                                                                                   <  2:548094/149‑1 (MQ=255)
                                                                   tGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCt                                                                                   >  2:710777/1‑148 (MQ=255)
                                                                   tGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCt                                                                                   <  1:710777/148‑1 (MQ=255)
                                                                     gAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCtcaa                                                                                >  1:1152388/1‑149 (MQ=255)
                                                                      aaaTGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCtcaat                                                                               <  1:1127668/149‑1 (MQ=255)
                                                                          ggCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCtcaatcacc                                                                           >  2:521585/1‑147 (MQ=255)
                                                                           gCTTTCCGCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTTACAAAATGGATGGTGCAAATGCACCATCCA‑TTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGCTGCCCAGCTCAATCAACTc                                                                         >  1:382535/1‑149 (MQ=255)
                                                                            cTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTc                                                                         <  1:865680/149‑1 (MQ=255)
                                                                                  tCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTAc                                                                   <  1:189140/149‑1 (MQ=255)
                                                                                  tCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTAc                                                                   >  2:858006/1‑149 (MQ=255)
                                                                                    aGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGTAGTAACAAAATGGATGGTGCAAATGCACCATCCATTATTAATGCAAGGCACAAAGTCGCTAGATGTTTGTCAGCTACCATGCCCCGCTAACTCGACTCCTTTagcg                                                                 >  2:1092989/1‑146 (MQ=255)
                                                                                      cgcTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTgg                                                               <  1:610085/149‑1 (MQ=255)
                                                                                       gcTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGa                                                              <  1:137849/149‑1 (MQ=255)
                                                                                           tggAAATACGCGCACTGCGATAAAAATTCGGGGAGAACTAACAAAATGGTTGTTGCATAAGAAACATCCATTTTTCATGCAATGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGc                                                          <  1:281250/148‑1 (MQ=255)
                                                                                           atgAAAAACGCGCACTTCGAGAAAAAGTCGCCGCGAAGTAACAAAATTGATGGTGCAAATTCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGc                                                          <  2:1063360/147‑1 (MQ=255)
                                                                                           aGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGc                                                          <  1:516806/149‑1 (MQ=255)
                                                                                              aaaaaCGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTg                                                       <  1:415729/149‑1 (MQ=255)
                                                                                              aaaaaCGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTg                                                       <  1:464750/149‑1 (MQ=255)
                                                                                               aaaaCGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTgg                                                      <  2:381561/149‑1 (MQ=255)
                                                                                                      gcACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCg                                               >  1:866841/1‑149 (MQ=255)
                                                                                                        aCTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGa                                             >  2:570225/1‑149 (MQ=255)
                                                                                                         cTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGaa                                            <  1:1098356/149‑1 (MQ=255)
                                                                                                          tGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAAc                                           >  1:990720/1‑149 (MQ=255)
                                                                                                           gCGAGAAAAAGTCGCCGCGAAGTAACAAACTGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAAcc                                          <  1:232935/149‑1 (MQ=255)
                                                                                                             gagaAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCgg                                        >  2:730817/1‑149 (MQ=255)
                                                                                                                aaaaaGTCGCCTCGAAGTAACAAAATGGATGGTGCAAATGCACGATCCATTTTTCATGCAAGGCAGAAAGTCGCGCGATGTTTGGCTGATTCTATGCCCAGCTCAATCATCTCCATTAGCTGGATTGCCTGGCTTGCCGGAACCGGAtt                                     <  1:164892/149‑1 (MQ=255)
                                                                                                                     gTCGCCGCGAAGTAACAAAATGTATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCAc                                <  1:856259/149‑1 (MQ=255)
                                                                                                                          cgcgAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCg                           <  1:833470/149‑1 (MQ=255)
                                                                                                                                gTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCAttt                     >  2:265572/1‑149 (MQ=255)
                                                                                                                                  tgaaaaaTGGATGGTGCATTTGCACCCTCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCCCCAtttca                   >  1:166366/4‑147 (MQ=255)
                                                                                                                                  aaCAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTaa                   <  1:460297/149‑1 (MQ=255)
                                                                                                                                       aaTGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCAt              >  2:420228/1‑149 (MQ=255)
                                                                                                                                        atggatggTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGAATTTCACCATCGCCATTTAACGCATc             >  2:543842/1‑149 (MQ=255)
                                                                                                                                         tggatggTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCTTTTAACGCATCa            >  1:857852/1‑149 (MQ=255)
                                                                                                                                            atggTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGACCCGGATTTTCCCCATCGCCATTTAACGCATCACGa         >  2:786124/1‑149 (MQ=255)
                                                                                                                                             tggTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCATCACGaa        <  1:627974/149‑1 (MQ=255)
                                                                                                                                              tgggctaaTGCACCTTTCATTTTTCATGCAATTCACAAATTCTCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTTGATTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCATCACGAAt       <  1:727527/143‑1 (MQ=255)
                                                                                                                                               gTGCATTTGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCATCACGAATa      >  1:147854/1‑149 (MQ=255)
                                                                                                                                                  cAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCAATACCTGGATTGCCTGGCTTGACGGAACCGGATTTTCCCCATCGCAATTTAACGCATCACGAATAGcc   >  1:200909/1‑149 (MQ=255)
                                                                                                                                                   aaaTGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCATCACGAATAGCCg  <  2:403035/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCATCACGAATAGCCG  >  NC_000913/1703103‑1703398

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑