Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,006,113 C→A 5.4% A54A (GCC→GCA)  zapC → FtsZ stabilizer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,006,1130CA 5.4% 252.6 / 3.2 86A54A (GCC→GCA) zapCFtsZ stabilizer
Reads supporting (aligned to +/- strand):  ref base C (35/46);  new base A (2/3);  total (37/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

CCAGACGATAACTGGCGTTGGTATTACGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAGCCGAACTGGTGCTGAATGCACTGGTAGCAATTCGCTATCTCAAACCGCAAATGCCGAAAAGCTGGCATTTTGTTTCCCATGGTG  >  NC_000913/1005964‑1006258
                                                                                                                                                     |                                                                                                                                                 
ccAGACGATAACTGGCGTTGGTATTACGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGAcgcc                                                                                                                                                   >  2:133450/1‑150 (MQ=255)
    aCGATAACAGGCGTTGGTATTACGATGAAGAGCACGATCGAATGAAGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATCCCTTTTCTCCCGCATGCTTTTGCGTTGACGACGCCgcgc                                                                                                                                               >  2:658873/1‑150 (MQ=255)
     cGATAACTGGCGTTGGTATTACGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGct                                                                                                                                              >  1:403443/1‑150 (MQ=255)
     cGATAACTGGCGTTGGTATTACGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGct                                                                                                                                              <  1:1474026/150‑1 (MQ=255)
     cGATAACTGGCGTTGGTATTACGATGAAGAGCACGATCGAATGATGCTCGATTTAGCAAATGGTATGCTATTTCGATCACGTTTTGCGCGAAAGATGTTGACTACAGATGAATTTTCACCCGCAGGCTTTTGCGTTGACTACCCcgcgat                                                                                                                                              >  1:701329/1‑148 (MQ=255)
       ttAACTGGCGTTGGTATTACGATGAAGAGCACGATCGAATGATGGTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGAGCAAGTTGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGctct                                                                                                                                            <  2:1386557/149‑1 (MQ=255)
          aCTGGCGTTGGTATTACGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTAtt                                                                                                                                         >  1:1021670/1‑150 (MQ=255)
          aCTGGCGTTGGTATTACGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTAtt                                                                                                                                         >  2:319697/1‑150 (MQ=255)
          aCTGGCGTTGGTATTACGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTAtt                                                                                                                                         >  2:1694299/1‑150 (MQ=255)
            tGGCGTTGGTATTACGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTtc                                                                                                                                       >  1:1583862/1‑150 (MQ=255)
             ggCGTTGGTATTACGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTtct                                                                                                                                      >  1:995214/1‑150 (MQ=255)
               cGTTGGTATTACGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCTCTCACGTTTTGCGCGCAAGATGTTGAATCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTctct                                                                                                                                    >  1:1083132/1‑150 (MQ=255)
               cGTTGGTATTACGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACCCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTctct                                                                                                                                    >  2:1332645/1‑150 (MQ=255)
                  tGGTATTACGCTGAAGAGCACGATCGAAAGATGCTCGATTTTGCCAATGGTCTGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCGCTATTTCTCtttt                                                                                                                                 <  2:1129678/150‑1 (MQ=255)
                  tGGTATTACGATGAAGAGCACGCTCGAATGATGCGCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCtttt                                                                                                                                 <  2:328541/150‑1 (MQ=255)
                  tGGTATTACGATGAAGAGCACGATCGAATGCTGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCtttt                                                                                                                                 <  1:1693497/150‑1 (MQ=255)
                  tGGTATTACGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCtttt                                                                                                                                 <  2:467407/150‑1 (MQ=255)
                   ggTATTACGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTg                                                                                                                                <  2:1262131/150‑1 (MQ=255)
                   ggTATTACGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTg                                                                                                                                <  1:1581051/150‑1 (MQ=255)
                      tttACGAGGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTgaag                                                                                                                             <  2:675869/149‑1 (MQ=255)
                      aTTACGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTgaag                                                                                                                             <  1:664151/150‑1 (MQ=255)
                          cGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGaaaa                                                                                                                         <  1:136035/150‑1 (MQ=255)
                            aTGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTga                                                                                                                               <  2:1045717/142‑1 (MQ=255)
                            aTGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTga                                                                                                                               >  1:1045717/1‑142 (MQ=255)
                               aaGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGcc                                                                                                                    <  2:1335704/150‑1 (MQ=255)
                               aaGAGCACGATCGAATGATGCTCGA‑TTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCcg                                                                                                                   >  2:663995/1‑150 (MQ=255)
                                agagCACGATCGAATGATTATAGATTTAGACAATTGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGAAGCCTTTTCTCCCGCAGGCTTTTGCGTTTACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCcg                                                                                                                   <  2:701329/150‑1 (MQ=255)
                                agagCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCcg                                                                                                                   <  2:1526788/150‑1 (MQ=255)
                                 gagCACGATCGAATGATGCTCGATTAAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCcgc                                                                                                                  <  2:1582510/150‑1 (MQ=255)
                                      cGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGAttt                                                                                                             >  2:796902/1‑150 (MQ=255)
                                       gATTGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGAtttt                                                                                                            >  2:747638/1‑150 (MQ=255)
                                        aTCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTa                                                                                                           >  1:622839/1‑150 (MQ=255)
                                          cGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAAt                                                                                                         <  1:744087/150‑1 (MQ=255)
                                           gAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAAtt                                                                                                        <  2:1338481/150‑1 (MQ=255)
                                             atgatgCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTa                                                                                                      <  2:501527/150‑1 (MQ=255)
                                              tgatgCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGGTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTAt                                                                                                     <  2:185517/150‑1 (MQ=255)
                                                      gATTTAGCCAATGGTATGCTATATCGCTCACGATTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGaa                                                                                             <  2:737612/150‑1 (MQ=255)
                                                           aGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACaaaa                                                                                        <  1:796902/150‑1 (MQ=255)
                                                            gCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGCCGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTTAAGAACaaaaa                                                                                       >  1:337890/1‑150 (MQ=255)
                                                            gCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACaaaaa                                                                                       >  2:604220/1‑150 (MQ=255)
                                                            gCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACaaaaa                                                                                       >  1:1394761/1‑150 (MQ=255)
                                                             ccAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAg                                                                                      <  2:995214/150‑1 (MQ=255)
                                                              cAATGGTATGCTATTTCGCTCCCGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGATTTTGCGTTGACGACGCCTCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAGc                                                                                     <  1:129638/150‑1 (MQ=255)
                                                              cAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAGc                                                                                     >  1:339114/1‑150 (MQ=255)
                                                              cAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAGc                                                                                     >  2:1005460/1‑150 (MQ=255)
                                                              cAATGGTATGCTATTTCGCTCACGTATTGCGCGCAAGATGTTGACTCCAGATGCC‑TTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTATTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAGcc                                                                                    >  1:1330165/1‑150 (MQ=255)
                                                                   gTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAGCCGAAc                                                                                >  1:1051477/1‑150 (MQ=255)
                                                                   gTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAGCCGAAc                                                                                >  2:1015501/1‑150 (MQ=255)
                                                                     ttGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAGCCGAACTg                                                                              <  1:319697/149‑1 (MQ=255)
                                                                     aTGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAGCCGAACTg                                                                              <  1:1694299/150‑1 (MQ=255)
                                                                       gCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAGCCGAACTGGt                                                                            <  2:974949/150‑1 (MQ=255)
                                                                       gCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAGCCGAACTGGt                                                                            >  2:994189/1‑150 (MQ=255)
                                                                         tATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAGCCGAACTGGTGc                                                                          >  1:839232/1‑150 (MQ=255)
                                                                         tATTCCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGAAGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAGCCGAACTGGTGc                                                                          <  2:900017/150‑1 (MQ=255)
                                                                              cGCTCACGTTTTGCGCGCAAGCTGTTGACTCCAGATGCCTTTTCTCCCGCCGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAGCCGAACTGGTGCTGAAt                                                                     >  2:222748/1‑150 (MQ=255)
                                                                              cGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAGCCGAACTGGTGCTGAAt                                                                     <  1:663995/150‑1 (MQ=255)
                                                                                     gTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAGCCGAACTGGTGCTGAATGCACTgg                                                              >  2:1275432/1‑150 (MQ=255)
                                                                                     gTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAGCCGAACTGGTGCTGAATGCACTgg                                                              <  2:965746/150‑1 (MQ=255)
                                                                                         tGCGCGCAAGATGTTGACCCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAGCCAAACTGGTGCTGAATGCACTGGTAGc                                                          <  1:1332645/150‑1 (MQ=255)
                                                                                               cAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAGCCGAACTGGTGCTGAATGCACTGGTAGCAATTCg                                                    <  1:1005460/150‑1 (MQ=255)
                                                                                               cAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCAGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAGCCGAACTGGTGCTGAATGCACTGGTAGCAATTCg                                                    >  2:1189874/1‑150 (MQ=255)
                                                                                               cAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCAGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAGCCGAACTGGTGCTGAATGCACTGGTAGCAATTCg                                                    >  2:1136849/1‑150 (MQ=255)
                                                                                                 aGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAGCCGAACTGGTGCTGAATGCACTGGTAGCAATTCGCt                                                  <  2:1583862/150‑1 (MQ=255)
                                                                                                 aGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAGCCGAACTGGTGCTGAATCCACTGGTAGCAATTAGCt                                                  >  2:1207533/1‑150 (MQ=255)
                                                                                                  gATGTTGACTCCAGATGCCATTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAGCCGAACTGGTGCTGAATGCACTGGTAGCAATTCGCTa                                                 <  2:724604/150‑1 (MQ=255)
                                                                                                            ccAGATGCCTTTTCTCCCGCAGGCTTTTGTGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAGCCGAACTGGTGCTGAATGCACTGGTAGCAATTCGCTATCTCAAACCg                                       >  2:1473407/1‑150 (MQ=255)
                                                                                                                  gCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAGCCGAACTGGTGCTGAATGCACTGGTAGCAATTCGCTATCTCAAACCGCAAATg                                 <  1:133450/150‑1 (MQ=255)
                                                                                                                  gCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAGCCGAACTGGTGCTGAATGCACTGGTAGCAATTCGCTATCTCAAACCGCAAATg                                 <  2:839232/150‑1 (MQ=255)
                                                                                                                    cTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAGCCGAACTGGTGCTGAATGCACTGGTAGCAATTCGCTATCTCAAACCGCAAATGcc                               >  2:753273/1‑150 (MQ=255)
                                                                                                                     ttttCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAGCCGAACTGGTGCTGAATGCACTGGTAGCAATTCGCTATCTCAAACCGCAAATGCCg                              >  2:870502/1‑150 (MQ=255)
                                                                                                                     ttttCTCCCGCAGGCTTTTGCGGTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAGCCGAACTGGTGCTGAATGCACTGGTAGCAATTCGCTATCTCAAACCGCAAATGCCg                              <  2:1083132/150‑1 (MQ=255)
                                                                                                                      tttCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAGCCGAACTGGTGCTGAATGCACTGGTAGCAATTCGCTATCTCAAACCGCAAATGCCGa                             <  1:578386/150‑1 (MQ=255)
                                                                                                                          tcCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGACGCGATTTTAATTTATCTAAAGAACAAAAAGCCTAACTGGTGCTGAATGCACTGGTAGCAATTCGCTATCTCAAACCGCAAATGCCGAAAAg                         <  2:1184598/150‑1 (MQ=255)
                                                                                                                             cGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAGCCGAACTGGTGCTGAATGCACTGGTAGCAATTCGCTATCTCAAACCGCAAATGCCGAAAAGCTg                      >  1:1556149/1‑150 (MQ=255)
                                                                                                                               cAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAGCCGAACTGGTGCTGAATGCACTGGTAGCAATTCGCTATCTCAAACCGCAAATGCCGAAAAGCTGGc                    <  1:1481255/150‑1 (MQ=255)
                                                                                                                                      ttGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAGCCGAACTGGTGCTGAATGCACTGGTAGCAATTCGCTATCTCAAACCGCAAATGCCGAAAAGCTGGCATTTTGt             <  2:1394761/150‑1 (MQ=255)
                                                                                                                                        gCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAGCCGAACTGGTGCTGAATGCACTGGTAGCAATTCGCTATCTCAAACCGCAAATGCCGAAAAGCTGGCATTTTGttt           <  1:1460100/150‑1 (MQ=255)
                                                                                                                                        gCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAGCCGAACTGGTGCTGAATGCACTGGTAGCAATTCGCTATCTCAAACCGCAAATGCCGAAAAGCTGGCATTTTGttt           <  2:622839/150‑1 (MQ=255)
                                                                                                                                         ctttGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGGGCCGCGATTTTAATTTATCTAAAGAACAAAAAGCCGAACTGGTGCTGAATGCACTGGTAGCAATTCGCTATCTCAAACCGCAAATGCCGAAAAGCTGGCATTTTGTTTc          <  1:1416058/148‑1 (MQ=255)
                                                                                                                                         cGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAGCCGAACTGGTGCTGAATGCACTGGTAGCAATTCGCTATCTCAAACCGCAAATGCCGAAAAGCTGGCATTTTGTTTc          >  2:936048/1‑150 (MQ=255)
                                                                                                                                          gTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAGCCGAACTGGTGCTGAATGCACTGGTAGCAATTCGCTATCTCAAACCGCAAATGCCGAAAAGCTGGCATTTTGTTTcc         <  2:1051477/150‑1 (MQ=255)
                                                                                                                                            tGACGTCGCAGAGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACATAAAGCCGAACTGGTGCTGAATGCACTGGTAGCAATTCGCTATCTCAAACCGCAAATGCCGAAAAGGTGGCATTTTGTTTccct       <  2:1330165/150‑2 (MQ=255)
                                                                                                                                            tGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAGCCGAACTGGTGCTGAATGCACTGGTAGCAATTCGCTATCTCAAACCGCAAATGCCGAAAAGCTGGCATTTTGTTTCCCa       <  1:658873/150‑1 (MQ=255)
                                                                                                                                             gacgacGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAGCCGAACTGGTGCTGAATGCACTGGTAGCAATTCGCTATCTAAAACCGCAAATGCCGAAAAGCTGGCATTTTGTTTCCCAt      >  1:1426260/1‑150 (MQ=255)
                                                                                                                                              acgacgCAGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAGCCGAACTGGTGCTGAATGCACTGGTAGCAATTCGCTATCTCAAACCGCAAATGCCGAAAAGCTGGCATTTTGTTTCCCATg     <  1:1136849/150‑1 (MQ=255)
                                                                                                                                              acgacgCAGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAGCCGAACTGGTGCTGAATGCACTGGTAGCAATTCGCTATCTCAAACCGCAAATGCCGAAAAGCTGGCATTTTGTTTCCCATg     <  1:1189874/150‑1 (MQ=255)
                                                                                                                                               cgacgCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAGCCGAACTGGTGCTGAATGCACTGGTAGCAATTCGCTATCTCAAACCGCAAATGCCGAAAAGCTGGCATTTTGTTTCCCATgg    >  1:211016/1‑150 (MQ=255)
                                                                                                                                                 acgCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAGCCGAACTGGTGCTGAATGCACTGGTAGCAATTCGCTATCTCAAACCGCAAATGCCGAAAAGCTGGCATTTTGTTTCCCATGGTg  >  1:983826/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
CCAGACGATAACTGGCGTTGGTATTACGATGAAGAGCACGATCGAATGATGCTCGATTTAGCCAATGGTATGCTATTTCGCTCACGTTTTGCGCGCAAGATGTTGACTCCAGATGCCTTTTCTCCCGCAGGCTTTTGCGTTGACGACGCCGCGCTCTATTTCTCTTTTGAAGAAAAGTGCCGCGATTTTAATTTATCTAAAGAACAAAAAGCCGAACTGGTGCTGAATGCACTGGTAGCAATTCGCTATCTCAAACCGCAAATGCCGAAAAGCTGGCATTTTGTTTCCCATGGTG  >  NC_000913/1005964‑1006258

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