Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,560,985 C→A 16.6% intergenic (+87/‑383) yffL → / → yffM CPZ‑55 prophage; uncharacterized protein/CPZ‑55 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,560,9850CA16.6% 149.6 / 2.6 78intergenic (+87/‑383)yffL/yffMCPZ‑55 prophage; uncharacterized protein/CPZ‑55 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (34/29);  new base A (7/8);  total (41/37)
Fisher's exact test for biased strand distribution p-value = 7.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCAGAAAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAA  >  NC_000913/2560841‑2561134
                                                                                                                                                |                                                                                                                                                     
ggCAGAAAGTGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAg                                                                                                                                                  <  1:229774/150‑1 (MQ=255)
ggCAGAAAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGGCAAAAGTTGCCGTTAGTGGTGGGTCAAAAg                                                                                                                                                  <  2:42393/150‑1 (MQ=255)
  cAGAAAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGtt                                                                                                                                                >  1:731227/1‑150 (MQ=255)
   aGAAAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTg                                                                                                                                               <  2:466852/150‑1 (MQ=255)
    gAAATGGAGTGCAACGAATACCAGTTTTTATTCAAGCCGTACAACATGAAATAACGGAAACAGGCGTTCAATTTTACCGACTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCATAAGTTGCCGTTTGTGGTGGGTCAAAAGTTGc                                                                                                                                              <  1:333173/150‑1 (MQ=255)
      aaGGGAGCGCAACGAATACCAGTTATTCTTCAAGCCGTACATCATGAAATAACGGCATCAGGCGCGCACTGGTCCCGCCTTTTCTTTGCCTTAATTCCATGTATTCAGGTGGGTCAAAAGGTGCCGTTAGTGGTGGGTCAAAAGTTGCCg                                                                                                                                            <  2:328617/150‑1 (MQ=255)
         ggAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAAAAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGGCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTa                                                                                                                                         <  1:1284012/150‑1 (MQ=255)
           aGTGCACCGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAg                                                                                                                                       <  1:1317569/150‑1 (MQ=255)
              gCAACGAATTCCAGTTATTATTCAAGCCGTACAACCTGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCATAATTCCATGAATTCAGGTGGGTCAAAACTTGCCGTTAGAGGTGGGTCAAAAGTTGCCGTTAACGTg                                                                                                                                    <  2:868308/150‑1 (MQ=255)
              gCAACGAATACCAGTTTTTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGATGCCGATAGTGTTGGGTAAAAATTTGCCGATAAGGTg                                                                                                                                    >  2:17972/1‑150 (MQ=255)
              gCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAATGGTGGGTCAAAAGTTGCCGTTAAGGTg                                                                                                                                    <  1:1237981/150‑1 (MQ=255)
               cAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTTGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAATTTGCCGTTAaggagg                                                                                                                                   >  2:970119/1‑150 (MQ=255)
               cAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTATTGGTGGGTAAAAAGTTGCAGTTTAGGTgg                                                                                                                                   >  1:664211/1‑150 (MQ=255)
               cAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTACTGTTGGGTCAAAAGATGCCTTGAaggagg                                                                                                                                   >  1:343410/1‑150 (MQ=255)
               cAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAATGGTGGGTCAAAAGTTGCCTTTAaggagg                                                                                                                                   >  1:670963/1‑150 (MQ=255)
               cAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGATAATGGTTGCTAAAAATTTGAATTTAAggagt                                                                                                                                   >  2:1403933/1‑147 (MQ=255)
               cAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCGTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCAGTTAAGGTgg                                                                                                                                   >  1:909296/1‑150 (MQ=255)
                 aCGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGGCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGt                                                                                                                                 <  2:854789/150‑1 (MQ=255)
                   gAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGGGGGTCAAAAGTTGCCGTTAAGGTGGGTCa                                                                                                                               <  1:1466059/150‑1 (MQ=255)
                   gAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAAAGGTTGGTAAAATTTAGCCGTAAAGGAGAGTCa                                                                                                                               >  2:1661480/1‑150 (MQ=255)
                    aaTACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTTATGGTGGGTCAAAATTTGACTTATAGGTGGGTaaa                                                                                                                              >  1:692727/1‑150 (MQ=255)
                       aCCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAg                                                                                                                           <  1:650672/150‑1 (MQ=255)
                         cAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGTTGGGTCAAAATATGCATTTAAGGTGAGTCAAAAttt                                                                                                                         >  2:644856/1‑150 (MQ=255)
                                 attCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGGTAGTGGTGGGTCAAATGTTGCGGTTACGGTGTGTCAAAAGTTggctttta                                                                                                                 >  1:90304/1‑143 (MQ=255)
                                   tCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTAAAAATTTGCAGTTAAGGAGGGTCAAAATTTGCCGTTaaag                                                                                                               >  2:999871/1‑148 (MQ=255)
                                   tCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAAATCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAgg                                                                                                               <  1:1105221/150‑1 (MQ=255)
                                    cAAGCCGTACAACATTAAATAAAGGAATCAGGCGTTCAATTTTCACGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTTGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGt                                                                                                              <  1:1018741/150‑1 (MQ=38)
                                    cAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTAAAAAGTTGCCGTGTAGGTGGGTCAAAATTTGCCGTTTAGGt                                                                                                              >  1:1140003/1‑150 (MQ=255)
                                     aaGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCTTTAAGGTGTGGCAAAATTTGCAGATAAGCTg                                                                                                             >  1:1240206/1‑150 (MQ=255)
                                     aaGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGAAGTTAAGGAGAGTCAAAATTTTCCGTTAACGTg                                                                                                             >  1:1144898/1‑150 (MQ=255)
                                      agacGTACAACATGAAATAACGGAAAGAGCCGGGCATTGTTCCCGCATGATCTATCCCATAATTCCATGAATTCAGGTGGGGAAACATCTGCCGTTACTGCTGGGTCAAAACTAGCCCTTCAGGTGGGGCAAAAGTTGCCGTTAAGGTgg                                                                                                            <  2:273678/147‑1 (MQ=255)
                                         cGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCAGTTTAGGTGGGtg                                                                                                         >  1:1595502/1‑149 (MQ=255)
                                         cGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGACGTTTAGGTGGTTc                                                                                                         >  1:1588034/1‑150 (MQ=255)
                                              cacaTGAAATCACGGAATCAGGCTGGCTATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTATCCGGTAGTGGTGGGTAAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAg                                                                                                    <  1:130935/149‑1 (MQ=38)
                                               acaTGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTATTGGTGGGTCAAATGTTGCAGTTAAGGAGGGTGAATAGTTGACGTAAAGGAGGCTCAAATGt                                                                                                   >  2:173117/1‑150 (MQ=255)
                                               acaTGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGTTGGGTCAAATGTTGCCTTTTAGGTGTGGAAAAAGTTGACGTATAGGTGTGTGaaattt                                                                                                   >  2:134791/1‑147 (MQ=255)
                                                     aaTAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGACGa                                                                                             >  1:1519857/1‑150 (MQ=255)
                                                     aaTAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTAAAAAGTTGCCGa                                                                                             >  1:1633407/1‑150 (MQ=255)
                                                      aTAACGTTATCAGGCGGGCAGTGATCCCGCCTTTTCTTGGCCTTAATTCCAGGAATTCCGTGGGGGAAAAAGATCCCGTTAGTGGTGTGGCTAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAg                                                                                            <  2:1683136/150‑1 (MQ=21)
                                                      aTAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAg                                                                                            <  1:1484822/150‑1 (MQ=255)
                                                               tCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTTAGGTGTGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGACGAGTAAGAAtta                                                                                   >  1:1766975/1‑149 (MQ=255)
                                                                cAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTATTGGTGGGTAAAAATTTGCCTTTAAGGTGGGTCAAAATTTGACGTATAGGTGTCGCAAAATTTGCAGAGTAACAAtcat                                                                                  >  2:624659/1‑147 (MQ=25)
                                                                 aGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGACGTTAAGGTGGGTCAAAAGTTGCCGTTTAGGTGGGTGAAAATTAGACGAGTAAGAATTcta                                                                                 >  2:880004/1‑150 (MQ=38)
                                                                   tatggCAATCATCCAGCCTTTTCTATTCCTTAATTTCATGAATTCATGTGGCTTTAAATTTGGCGTTTGTGGGGGGTCAATAGTTGCCTTTAAGTTGGGTCAAAAGTTGCATTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTctact                                                                               <  2:1477193/147‑1 (MQ=19)
                                                                     gggCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTGAAGGTGGGTCATAAGTTGCCGAGTAAGAATTCTACTAt                                                                             >  1:1375059/1‑150 (MQ=255)
                                                                         attGTTTCCGCCTTTCGTTTGCCTTAAATTAATGAATTCATGTGGGTTAAAAGTTGGCCTTTGTGGTGGGTAAATAGTTGCAGTTCAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGGCAAAAGTTGCCGAGTAAGAATTCTACTATTAAg                                                                         <  1:1180974/148‑1 (MQ=25)
                                                                                ccGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTATTGGTGGGTCAAAATTTGCCGTTAAGGTGTGTAAAAATTTGACGTTAAGGTGGGGGAAATGTTGACGAGTAAGAAGTCTACAATTATGAAGTCTa                                                                  >  1:356144/1‑150 (MQ=25)
                                                                                           tacccatAAATCCATGAATTCAGGTGGGTCAAAAGGCGCCCTTATTTGTGGGTCAAAAGATGCCTTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAg                                                       <  2:1093276/144‑1 (MQ=37)
                                                                                               catAATTTCATGAATTCTTGTTGCTAAAAATTTGCCTTTACTTGGGGGTATATAGTTGCCTTTTAGGTGGGTCAAAAGTTGCAGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGtttt                                                   <  2:762300/148‑1 (MQ=25)
                                                                                                   aatCCATGAATTCAGGGGGGTCAAAAGTTGCCGTTAGTGGTGGGTAAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAg                                               <  1:1341957/148‑1 (MQ=255)
                                                                                                   aTTACATGAATTCAGGGTGGTGAACAGTTGCCGTTTGTGGTGGGTAAAATGTTTCCGTTAAGGTGGTGCAAAAGTTGCCTTTAAGGTGGGTCAATAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACTGGAGTTTTTTAg                                               <  1:1778493/150‑1 (MQ=255)
                                                                                                    taccATGAATTCTGTGTGGTCAAAAGTTGCCGTTAGTTGTGGGTCAAATGTTTCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGGCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGc                                              <  1:1661480/148‑1 (MQ=37)
                                                                                                    atCCCTGTATTCATGCGGGCAAAAAGATCGCCTTATTGGTGGCTAAAAACCTGCCTTTAAGGTGGGTCAAAAGTTGCCTTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGc                                              <  1:1067494/149‑1 (MQ=25)
                                                                                                      ccATGAATTCAGGTGGGTCAAAAGTTGCCTTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCt                                            <  2:1603641/150‑1 (MQ=255)
                                                                                                         tGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAATGTTGCCGTTAAGGAGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAAATCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTccc                                         >  1:1668629/1‑150 (MQ=255)
                                                                                                             ttCAGGTGTGTCAAAAGTTGCCGTTAGTGGGGGGTCAAAAGATGCCGTTAAGGTGGGTCAAAAGTAGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCcggtc                                     <  1:722874/150‑1 (MQ=255)
                                                                                                              tCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTAAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCcggtcg                                    <  1:618809/150‑1 (MQ=255)
                                                                                                               cAGGTGGGGAAAAAGTAGCCGTTAGTGGTGGGTCAATAGTTGCCTTTAAGGTGGTTCAAATGTTGCCGTTTAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCcggtcgg                                   <  1:997705/150‑1 (MQ=38)
                                                                                                               cAGGGGGGTCCAAAGATGCCCTTAGGGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCcggtcgg                                   <  1:838114/150‑1 (MQ=37)
                                                                                                                 ggTGGGTCAAAAGATGCCGTTAGTGGGGGGTCACAAGTTGCCGTAAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCggtcggtc                                 <  1:1352821/150‑1 (MQ=255)
                                                                                                                     ggTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTAAAAAGTTGCAGTTAAGGTGGGTCAAAAGTTGACGAGTAAGAAGTCTACAATTAAGAAATCTAGGCGACGGGAGTTGTTTATCCTACCGGACCGACATag                             >  1:881454/1‑150 (MQ=25)
                                                                                                                      gTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATaga                            >  1:669918/1‑150 (MQ=255)
                                                                                                                      gTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGACTAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGACGGTCATaga                            >  2:402045/1‑150 (MQ=255)
                                                                                                                        cAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGAGGGTCAAAATTTGCCTAGTAAGAATTCTACTATTAAGAATTCTAGCCGAAGCGAGTTTTTTAGCCTCCCGGTCGGTCATagaga                          >  1:1150419/1‑150 (MQ=38)
                                                                                                                        cAAAAGGTGCCGTTAGGGGTGGGTCAAAAGTCGCCGTTACGGTGGGTCAAAAGTTGCCGTTAAGGTGGGGCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATagaga                          <  1:1795723/150‑1 (MQ=38)
                                                                                                                         aaaaGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGATTTTTAGCCTCCCGGTCCGTCATagagac                         >  1:585870/1‑149 (MQ=255)
                                                                                                                           aaGTTGCCTTTACTGGTGGGTCAATAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTc                       <  2:1866805/150‑1 (MQ=255)
                                                                                                                             gTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGc                     >  1:1611673/1‑150 (MQ=255)
                                                                                                                                 ccGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTag                 >  2:776069/1‑150 (MQ=255)
                                                                                                                                  cGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCATAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTATACTCCCGGTCGGTCATAGAGAGTCGCATaga                >  2:1006197/1‑150 (MQ=38)
                                                                                                                                   tttATTGGGGGGTAAAATCTTGCCTTTTAGGTGGGTCAAAAGTTGCCTTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTagag               <  2:685941/149‑1 (MQ=38)
                                                                                                                                   gttacgtgtggGTCATATCTTGCCGTTTAGGTGGGGTAATAGTTTCCGTTAAGGGGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTagag               <  2:1144925/143‑1 (MQ=255)
                                                                                                                                     tAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGtt             >  2:1443314/1‑150 (MQ=255)
                                                                                                                                      cggggtggGTCAAAACTTGCCGTTACGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATCCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTa            <  1:624659/147‑1 (MQ=255)
                                                                                                                                      aGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGGCGGTCATAGAGAGTCGCCTAGAGTTa            >  2:1023054/1‑150 (MQ=255)
                                                                                                                                       gtggtgGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTaa           >  2:1803283/1‑150 (MQ=255)
                                                                                                                                              gTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACagaag    >  1:315336/1‑150 (MQ=255)
                                                                                                                                              gTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGACTAGAGTTTAACagaat    >  2:1380704/1‑149 (MQ=255)
                                                                                                                                                cAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAgaagaa  >  1:415387/1‑150 (MQ=255)
                                                                                                                                                cAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAgaagaa  >  1:1701104/1‑150 (MQ=255)
                                                                                                                                                cAAAAGTTGCCGTTAAGGTGGCTCAAAACTTGCCGTTAAGGTGGGGCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAgaagaa  <  1:1609579/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
GGCAGAAAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAA  >  NC_000913/2560841‑2561134

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: