Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,467,120 A→G 46.8% intergenic (‑129/+40) bfd ← / ← chiA bacterioferritin‑associated ferredoxin/periplasmic endochitinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,467,1200AG46.8% 12.7 / 122.9 96intergenic (‑129/+40)bfd/chiAbacterioferritin‑associated ferredoxin/periplasmic endochitinase
Reads supporting (aligned to +/- strand):  ref base A (25/26);  new base G (28/17);  total (53/43)
Fisher's exact test for biased strand distribution p-value = 2.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.39e-01

Only 100 of 101 total aligned reads displayed.

AAATGAAAATAGTTCTTATTTCAATTACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTACTTCGGTGCCTGTTCTTT  >  NC_000913/3467021‑3467219
                                                                                                   |                                                                                                   
aaatGATAATAGTACGTATTTAAAATACTGAAAAAATCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGa                                                                                                     <  1:1796246/100‑1 (MQ=255)
aaatGAAAATAGTTCTTATTTCAATTACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGa                                                                                                     <  2:920460/100‑1 (MQ=255)
 aatGAAAATAGTTCTTATTTCAATTACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGAc                                                                                                    >  2:1189330/1‑100 (MQ=255)
 aatGAAAATAGTTCTTATTTCAATTACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGAc                                                                                                    >  2:2575618/1‑100 (MQ=255)
 aatGAAAATAGTTCTTATTTCAATTACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGAc                                                                                                    <  1:450464/100‑1 (MQ=255)
 aatGAAAATAGTTCTTATTTCAATTACGGAAAGACTCTTTATTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGAc                                                                                                    <  1:380108/100‑1 (MQ=255)
 aatGAAAATAGTACTAATTTCAATTACGGAAAGACTCTTTTTTAAGCAAATAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGa                                                                                                     <  1:1796247/99‑1 (MQ=255)
   tGAAAATAGTTCTTATTTCAATTACGGAAAGACTCTTTTTTAAGCAAAGACTAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCaa                                                                                                  >  1:2409485/1‑100 (MQ=255)
   tGAAAATAGTTCTTATTTCAATTACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACaa                                                                                                  <  1:2050460/100‑1 (MQ=255)
   tGAAAATAGTTCTTATTTCAATTACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACaa                                                                                                  >  2:434380/1‑100 (MQ=255)
    gAAAATAGTTCTTATTTCAATTACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCaaa                                                                                                 <  2:1627336/100‑1 (MQ=255)
     aaaaTAGTTCTTATTTCAATTACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCaaaa                                                                                                >  2:1734437/1‑100 (MQ=255)
     aaaaTAGTTCTTATTTCAATTACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCaaaa                                                                                                >  1:417259/1‑100 (MQ=255)
        aTAGTTCTGATTTCTATTACGAAAAGACACTTTTTTAAGCAAGGAAGAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACaaa                                                                                                 <  1:1562199/96‑1 (MQ=255)
         tAGTTCTTATTTCAATTACGGAAAGACGCTTTTTTAAGCAAAGAATGGAAATATGAATAAAAATAATGGGCTACCAAATGTAGACAATTAGCAAAAAATg                                                                                            >  1:1432500/1‑100 (MQ=255)
          aGTTCTTATTTCAATTACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATgc                                                                                           <  2:1002935/100‑1 (MQ=255)
           gTTCTTATTTCAATTACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATgcg                                                                                          >  2:1612627/1‑100 (MQ=255)
            ttCTTATTTCAATTACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATAGGCTACCAAAAGTAGCCCATTGACAAAAAATg                                                                                            >  1:493424/1‑97 (MQ=255)
             tCTTATTTCAATTACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATgcggc                                                                                        <  2:1764472/100‑1 (MQ=255)
             tCTTATTTCAATTACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATgcggc                                                                                        <  1:66584/100‑1 (MQ=255)
              cTTATTTCAATTACCGAAAGACTCTTTTTTAAGCAAAGAATACAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATgcggcg                                                                                       <  2:545999/100‑1 (MQ=255)
                tATTTCAATTACGGAAAGACTCTTTTTTAAGCAAAAAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGAAAATAAAAAAAAAATgcggcg                                                                                       <  1:171341/98‑1 (MQ=255)
                    tCAATTACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTg                                                                                 >  2:1968190/1‑100 (MQ=255)
                     cAATTACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTgg                                                                                <  1:1010907/100‑1 (MQ=255)
                     cAATTACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTgg                                                                                <  2:1201155/100‑1 (MQ=255)
                     cAATTACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTgg                                                                                >  2:121506/1‑100 (MQ=255)
                      aaTTACGGAAAGAATCTTTTTTAAGCAAACAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAAAAAAAAAAAAAAAAATGCGGCGATACTgg                                                                                <  1:1310529/99‑1 (MQ=255)
                       aTTACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGaa                                                                              >  2:35501/1‑100 (MQ=255)
                       aTTACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGaa                                                                              <  2:1073055/100‑1 (MQ=255)
                       aTTACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGaa                                                                              <  1:2506349/100‑1 (MQ=255)
                       aTTACGGAAAGAATCTGTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCTAAAGTAGAAAAAAAAAAAAAAATGCGGCGATACTGGaa                                                                              <  1:2134933/100‑1 (MQ=255)
                         tACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAgg                                                                            <  2:1587130/100‑1 (MQ=255)
                               aaGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGc                                                                      <  2:565989/100‑1 (MQ=255)
                               aaGACTCTTTTTTAAGCAAAGAATAGAAATATAAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGc                                                                      <  2:695172/100‑1 (MQ=255)
                                aGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGcc                                                                     <  1:97591/100‑1 (MQ=255)
                                aGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGcc                                                                     <  2:1312624/100‑1 (MQ=255)
                                aGACTCTTTTTTAAGCAAAGAATAGAAATATAAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGcc                                                                     <  2:871251/100‑1 (MQ=255)
                                  aCTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCaa                                                                   <  1:1274106/100‑1 (MQ=255)
                                   ctctTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTCCCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAac                                                                  >  1:1440968/1‑100 (MQ=255)
                                   ctctTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAac                                                                  >  2:337416/1‑100 (MQ=255)
                                   ctctTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAac                                                                  >  1:1351462/1‑100 (MQ=255)
                                   ctctTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAac                                                                  >  2:830480/1‑100 (MQ=255)
                                   ctctTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAac                                                                  >  1:2398337/1‑100 (MQ=255)
                                   atctTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTTGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAac                                                                  >  2:1818836/2‑100 (MQ=255)
                                    tctTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAaca                                                                 >  2:446679/1‑100 (MQ=255)
                                     ctTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAacac                                                                <  2:191585/100‑1 (MQ=255)
                                     ctTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGATATCGCCAaaac                                                                >  1:2521943/1‑100 (MQ=255)
                                      ttttttAAGCAAAGAATAGAAATATGAATAAAAATGATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAAcaca                                                               >  1:1793527/1‑100 (MQ=255)
                                      ntttttAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCAATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAAcaca                                                               >  1:789279/2‑100 (MQ=255)
                                         tttAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACAttt                                                            >  2:1586193/1‑100 (MQ=255)
                                          ttAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTa                                                           >  1:1638330/1‑100 (MQ=255)
                                              gCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTAtt                                                         >  1:1253005/1‑98 (MQ=255)
                                              gCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGc                                                       >  2:2544086/1‑100 (MQ=255)
                                              gCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGc                                                       >  1:1483257/1‑100 (MQ=255)
                                              gCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGc                                                       <  2:1190468/100‑1 (MQ=255)
                                                aaaGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCtt                                                     >  2:1549653/1‑100 (MQ=255)
                                                 aaGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTa                                                    <  2:614740/100‑1 (MQ=255)
                                                  aGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAg                                                   <  1:1283064/100‑1 (MQ=255)
                                                  aGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAg                                                   <  1:2581741/100‑1 (MQ=255)
                                                   gAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGt                                                  <  1:2377459/100‑1 (MQ=255)
                                                   gAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGt                                                  <  2:1900429/100‑1 (MQ=255)
                                                           atatGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGTGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGgcg                                          <  2:787528/100‑1 (MQ=255)
                                                            tatGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGgcgc                                         >  2:1666923/1‑100 (MQ=255)
                                                               gAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGaa                                      <  2:692458/100‑1 (MQ=255)
                                                               gAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGaa                                      >  1:81543/1‑100 (MQ=255)
                                                               gAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGaa                                      >  1:2420441/1‑100 (MQ=255)
                                                                aaTAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCACCAACACATTTATTGCTTAGTAAACGGCGCGAAg                                     >  1:278209/1‑100 (MQ=255)
                                                                aaTAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAg                                     >  2:191140/1‑100 (MQ=255)
                                                                aaTAAAAATAATGGGCTACCAAAAGTAGAACATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAg                                     <  2:1246987/100‑1 (MQ=255)
                                                                  tAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGAt                                   >  2:860966/1‑100 (MQ=255)
                                                                   aaagatgaggggTTACCAAAAGTGGACGATTGACAAAAAATGCGGCGATAGTGGAAGGTATCGCCAACACATTTATTTTTTAGTAGACGGCGCGAAGa                                    <  1:2141092/89‑1 (MQ=255)
                                                                    aaaaTAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGCtt                                  >  2:2344805/1‑99 (MQ=255)
                                                                     aaaTAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTc                                <  1:1959753/100‑1 (MQ=255)
                                                                       ataatGGGCTACCAAAAGTAGCCCATTTGCAAAAAAGGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGc                              <  1:1372895/100‑1 (MQ=255)
                                                                       ataatGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGc                              <  2:2523868/100‑1 (MQ=255)
                                                                       ataatGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGc                              <  2:718258/100‑1 (MQ=255)
                                                                             ggCTACCAAAAGTAGCTCATTGACAAAGAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTTCTTAGTAAACGGCGCGAAGATTTCGCTAAATg                        <  2:2456700/100‑1 (MQ=255)
                                                                             ggCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGCTATCGCCAACACATTTATTGCTTACTAAACGACGCGAAGATTTCGCTAAATg                        >  1:2346320/1‑100 (MQ=255)
                                                                              gCTACCAAAAGTAGCACAAAGACAAAAAATGCGACGATACTGGAAGGTATCGCCAACACATTTACTGCTTAATAAACGGCGCGAAGATTTCGCTAAATg                        >  1:746096/1‑99 (MQ=255)
                                                                               cTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCg                      >  1:975363/1‑100 (MQ=255)
                                                                               cTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCg                      <  1:598620/100‑1 (MQ=255)
                                                                               cTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCg                      >  1:2427249/1‑100 (MQ=255)
                                                                                tACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGt                     >  1:669989/1‑100 (MQ=255)
                                                                                 acttttggTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTa                    >  1:553428/8‑100 (MQ=255)
                                                                                 aCCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTa                    >  1:2037702/1‑100 (MQ=255)
                                                                                 aCCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTa                    >  2:1000608/1‑100 (MQ=255)
                                                                                 aCCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTa                    >  1:1444147/1‑100 (MQ=255)
                                                                                 aCCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTa                    >  2:2043410/1‑100 (MQ=255)
                                                                                  ccAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTAc                   >  2:741750/1‑100 (MQ=255)
                                                                                   cAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTACt                  <  1:1722230/100‑1 (MQ=255)
                                                                                    aaaaGTAGCCCATTGACAAAGAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCTTACtt                 >  1:1588185/1‑100 (MQ=255)
                                                                                     aaaGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTACTTc                <  2:1985454/100‑1 (MQ=255)
                                                                                        gTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGCAAACGGCGCGAAGATTTCGCTAAATGCGTACTTCGGt             <  1:1385603/100‑1 (MQ=255)
                                                                                           gCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTACTTCGGTGcc          >  2:1005158/1‑100 (MQ=255)
                                                                                             ccATTGGCAAAGAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTACTTCGGTGCCTg        >  1:1312426/1‑100 (MQ=255)
                                                                                             ccATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTACTTCGGTGCCTg        >  2:469430/1‑100 (MQ=255)
                                                                                             ccATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTACTTCGGTGCCTg        >  2:657683/1‑100 (MQ=255)
                                                                                               aTTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTACTTCGGTGCCTGtt      >  2:2471508/1‑100 (MQ=255)
                                                                                                   aCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTACTTCGGTGCCTGTTCttt  <  1:2557947/100‑1 (MQ=255)
                                                                                                   aCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTACTTCGGTGCCTGTTCttt  <  2:1046532/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
AAATGAAAATAGTTCTTATTTCAATTACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTACTTCGGTGCCTGTTCTTT  >  NC_000913/3467021‑3467219

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: