Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,122,382 A→G 37.2% *320R (TGA→CGA)  ftsN ← essential cell division protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,122,3820AG37.2% 84.5 / 108.2 110*320R (TGA→CGA) ftsNessential cell division protein
Reads supporting (aligned to +/- strand):  ref base A (39/30);  new base G (20/21);  total (59/51)
Fisher's exact test for biased strand distribution p-value = 5.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.29e-01

Only 100 of 112 total aligned reads displayed.

GTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTCTCTTTGCCTTTCACCGGGCCAATGA  >  NC_000913/4122283‑4122480
                                                                                                   |                                                                                                  
gTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCa                                                                                                    >  2:1605120/1‑100 (MQ=255)
gTCACGAGCTGAACCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCa                                                                                                    <  1:1135438/100‑1 (MQ=255)
  cACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAAc                                                                                                  >  1:707090/1‑100 (MQ=255)
  cACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAAc                                                                                                  <  2:376181/100‑1 (MQ=255)
  cACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAAc                                                                                                  >  1:1823180/1‑100 (MQ=255)
         tGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCcggc                                                                                           <  1:2074688/100‑1 (MQ=255)
          gACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGGTTCAACCCCcggcg                                                                                          >  1:444434/1‑100 (MQ=255)
            ccccTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGTGTTTCGACCCCcggcg                                                                                          >  1:318208/1‑98 (MQ=255)
             cccTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCGACCCCCGgcggcg                                                                                       <  1:79392/100‑1 (MQ=255)
             cccTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGgcggcg                                                                                       <  2:2208663/100‑1 (MQ=255)
              ccTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCGACCACCGGCGgaga                                                                                      >  1:198335/1‑100 (MQ=255)
               cTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCGACCCCCGGCGGCGAg                                                                                     >  1:246705/1‑100 (MQ=255)
               cTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAg                                                                                     >  2:1352860/1‑100 (MQ=255)
                ttGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAg                                                                                     <  1:1979666/99‑1 (MQ=255)
                 tGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGcc                                                                                   <  2:2242480/100‑1 (MQ=255)
                    ttACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAAt                                                                               <  1:1620440/100‑1 (MQ=255)
                     tACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCGACCCCCGGCGGCGAGCCGAAt                                                                               <  1:782636/100‑1 (MQ=255)
                      aCGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCGACCCCCGGCGGCGAGCCGAATg                                                                              <  1:151911/100‑1 (MQ=255)
                      aCGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGc                                                                             >  1:580375/1‑100 (MQ=255)
                       cGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCGACCCCCGGCGGCGAGCCGAATGc                                                                             <  2:1284905/100‑1 (MQ=255)
                       cGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCGACCCCCGGCGGCGAGCCGAATGc                                                                             >  1:1283173/1‑100 (MQ=255)
                           tACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGtt                                                                         <  2:217928/100‑1 (MQ=255)
                           tACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGtt                                                                         <  2:217784/100‑1 (MQ=255)
                            aCAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCGACCCCCGGCGGCGAGCCGAATGCAGttt                                                                        <  2:1127790/100‑1 (MQ=255)
                            aCAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCGACCCCCGGCGGCGAGCCGAATGCAGttt                                                                        >  2:1168742/1‑100 (MQ=255)
                             cAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGttt                                                                        >  1:1523147/1‑99 (MQ=255)
                             cAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTg                                                                       >  2:511578/1‑100 (MQ=255)
                             cAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTg                                                                       >  2:261018/1‑100 (MQ=255)
                             cAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTg                                                                       >  1:1535454/1‑100 (MQ=255)
                              agagTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGt                                                                      <  2:520790/100‑1 (MQ=255)
                               gagTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTa                                                                     <  1:979669/100‑1 (MQ=255)
                               gagTACAAACCCCGTACAAAAGTACGGGGCATAACGCAATTATAGATGGGGGGGATTTTTAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTa                                                                     >  2:836667/1‑100 (MQ=255)
                                agTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTAt                                                                    >  1:161651/1‑100 (MQ=255)
                                 gTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACACCCGGCGGCGAGCCGAATGCAGTTTGTATg                                                                   >  2:1312899/1‑100 (MQ=255)
                                  tACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGa                                                                  >  1:854614/1‑100 (MQ=255)
                                    cAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGAcc                                                                <  2:804464/100‑1 (MQ=255)
                                    cAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGAcc                                                                >  1:116850/1‑100 (MQ=255)
                                    cAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGAcc                                                                >  1:378331/1‑100 (MQ=255)
                                    cAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGAcc                                                                >  2:144903/1‑100 (MQ=255)
                                         cgcGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGcc                                                            <  2:854556/96‑1 (MQ=255)
                                            gTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCtt                                                       <  2:689683/100‑1 (MQ=255)
                                              aCAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCtt                                                       >  2:1321698/1‑99 (MQ=255)
                                              aCAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTc                                                      <  1:1997844/100‑1 (MQ=255)
                                              aCAAAAGTACAGGGCATAATTAAATTATAGATAGCGAGGAATTAGAGGGATCCAACCCCCGGCCGCGAGCCGCAGGCAGTTTGTATGACACGCCctct                                                        >  2:1279348/1‑98 (MQ=255)
                                               cAAAAGTACGGGGCATAATGTAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCa                                                     >  2:329569/1‑100 (MQ=255)
                                               cAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCGACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCa                                                     <  1:2090276/100‑1 (MQ=255)
                                               cAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATGTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCa                                                     >  1:2071757/1‑100 (MQ=255)
                                                 aaaGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCGACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCTGCCATCTTCAAc                                                   <  1:1473989/100‑1 (MQ=255)
                                                 aaaGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAAc                                                   <  1:772980/100‑1 (MQ=255)
                                                 aaaGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAAc                                                   >  2:124929/1‑100 (MQ=255)
                                                 aaaGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGAAGTTTGTATGACCCGCCATCTTCaa                                                    >  1:1509649/1‑99 (MQ=255)
                                                     tACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGAt                                               <  1:1719239/100‑1 (MQ=255)
                                                     tACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGAt                                               <  1:45323/100‑1 (MQ=255)
                                                     tACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGAt                                               <  1:398230/100‑1 (MQ=255)
                                                        ggggCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCGACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGAtt                                              <  1:1240607/98‑1 (MQ=255)
                                                        ggggCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCGACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGa                                            <  2:554298/100‑1 (MQ=255)
                                                        ggggCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCGACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGa                                            >  1:1981969/1‑100 (MQ=255)
                                                        ggggCATAAGGCAATTATAGATGGGGGGGATTGTGATGGTTTCAACCCACGGCGGCGAGCCGAATGAAGATTGTATGACCCGCCATCTTCAACCGAATg                                             >  2:1587718/1‑99 (MQ=255)
                                                         gggCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCGACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAg                                           >  2:1662078/1‑100 (MQ=255)
                                                         gggCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCTCGGCGGCGAGCCTAATGCAGTGTGTATGACCCGCCATCTTCAACCGATTGAg                                           >  2:101264/1‑100 (MQ=255)
                                                         gggCATAATGCAATTATAGATGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGCTTGTATGACCCGCCATCTTCAACCGATTGAgg                                          >  2:406732/1‑100 (MQ=255)
                                                          ggCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCGACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAgg                                          <  1:710887/100‑1 (MQ=255)
                                                             aTAATGCAATTATAGATGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGc                                      <  2:597206/100‑1 (MQ=255)
                                                              tAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCGACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGc                                      <  2:1785269/100‑1 (MQ=255)
                                                              tAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCGACCCCCGGCGGCGAGCCGAATGCAGTGTGTATGACCCGCCATCTTCAACCGATTGAGGGTGc                                      <  2:1048025/100‑1 (MQ=255)
                                                              tAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGc                                      <  2:854614/100‑1 (MQ=255)
                                                              tAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGc                                      >  1:889487/1‑100 (MQ=255)
                                                              tAATGCAATTATAGATGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCt                                     >  1:239854/1‑100 (MQ=255)
                                                              tAACGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGACGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGTGTGc                                      >  2:1761302/1‑100 (MQ=255)
                                                                aTGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCt                                     >  2:469177/1‑99 (MQ=255)
                                                                 tGCAATTATAGATGGGGGGGATTTTGAGGGTTTCGACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGt                                   >  2:1036032/1‑100 (MQ=255)
                                                                    aaTTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTg                                <  2:1221888/100‑1 (MQ=255)
                                                                      ttATAGATGGGGGGGATTTTGAGGTTTTCGACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCg                              <  2:682826/100‑1 (MQ=255)
                                                                       tataGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCTGCCCTCTTCAACCTTTTGAGGGTGCTGTCTGCGt                             <  2:2075164/100‑1 (MQ=255)
                                                                        ataGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGtt                            >  1:1895884/1‑100 (MQ=255)
                                                                          aGATGGGGGGGATTTTGAGGGTTTCGACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTtct                          >  2:636744/1‑100 (MQ=255)
                                                                           gATGGGGGGGATTTTGAGGGTTTCGACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTtctc                         >  1:1416428/1‑100 (MQ=255)
                                                                           gATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTtctc                         <  2:116850/100‑1 (MQ=255)
                                                                            aTGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTctct                        <  1:1713416/100‑1 (MQ=255)
                                                                             tGGGGGGGATTTTGAGGGTTTCGACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTCTCtt                       >  1:1563064/1‑100 (MQ=255)
                                                                               ggggggATTTTGAGGGTTTCGACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTCTCTTTg                     >  2:721630/1‑100 (MQ=255)
                                                                                gggggATTTTGAGGGTTTCGACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTCTCTTTGc                    <  1:1594799/100‑1 (MQ=255)
                                                                                gggggATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTCTCTTTGc                    >  2:1005720/1‑100 (MQ=255)
                                                                                gggggATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTCTCTTTGc                    >  1:1271897/1‑100 (MQ=255)
                                                                                 ggggATTTTGAGGGTTTCGACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTCTCTTTGcc                   <  1:1895649/100‑1 (MQ=255)
                                                                                 ggggATTTTGAGGGTTTCGACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTCTCTTTGcc                   <  2:1537716/100‑1 (MQ=255)
                                                                                    gATTTTGAGGGTTTCGACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTCTCTTTGCCttt                <  1:960182/100‑1 (MQ=255)
                                                                                      ttttGAGGGTTTCGACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTCTCTTTGCCTTTCa              >  2:155811/1‑100 (MQ=255)
                                                                                       tttGAGGGTTTCGACCCCCAGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATAGAGGGTGCTGTCTGCGTTCTCTTTGCCTTTc               >  2:493837/1‑98 (MQ=255)
                                                                                        ttGAGGGTTTCGACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTCTCTTTGCCTTTCAcc            <  1:401999/100‑1 (MQ=255)
                                                                                        ttGAGGGTTTCGACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTCTCTTTGCCTTTCAcc            >  2:839087/1‑100 (MQ=255)
                                                                                        ttGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTCTCTTTGCCTTTCAcc            <  1:511577/100‑1 (MQ=255)
                                                                                          gAGGGTTTCGACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTCTCTTTGCCTTTCACCgg          <  1:725700/100‑1 (MQ=255)
                                                                                          gAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTCTCTTTGCCTTTCACCgg          <  1:2050897/100‑1 (MQ=255)
                                                                                          gAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTCTCTTTGCCTTTCACCgg          >  1:1124500/1‑100 (MQ=255)
                                                                                            gggTTTCGACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTCTCTTTGCCTTTCACCGGGc        >  1:464910/1‑100 (MQ=255)
                                                                                            gggTTTCGACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTCTCTTTGCCTTTCACCGGGc        <  2:1989568/100‑1 (MQ=255)
                                                                                             ggTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTCTCTTTGCCTTTCACCGGGcc       >  1:328691/1‑100 (MQ=255)
                                                                                             ggTTTCAACCCCCGGCGGCGAGCCGAAAGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTCTCTTTTCCTTTAACCGGGcc       >  2:1863289/1‑100 (MQ=255)
                                                                                                  cAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTCTCTTTGCCTTTCACCGGGCCAATGa  <  1:1352866/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
GTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTCTCTTTGCCTTTCACCGGGCCAATGA  >  NC_000913/4122283‑4122480

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