Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,560,996 G→T 16.7% intergenic (+98/‑372) 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,9960GT16.7% 147.8 / 2.9 82intergenic (+98/‑372)yffL/yffMCPZ‑55 prophage; uncharacterized protein/CPZ‑55 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (36/29);  new base T (8/8);  total (44/38)
Fisher's exact test for biased strand distribution p-value = 7.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCG  >  NC_000913/2560847‑2561145
                                                                                                                                                     |                                                                                                                                                     
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)
                                                                                                                                        gttaataGTTGCCTTTTAGGTGGGTCAAAAGTTGCCTTTCAGTGGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACagaag               <  2:90304/144‑1 (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)
                                                                                                                                           aaaaGTTGCCGTTATGGTGGGTCAAAAGTTGCCTTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAat            <  2:1502711/150‑1 (MQ=255)
                                                                                                                                             aaGATGCCTTTAAGGTGGGTCAAAAGGTGCCCTAAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAatta          <  2:1373973/150‑1 (MQ=255)
                                                                                                                                               gTTGCCGTTAAGGTGGGTAAAAAGTTTCCTTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAattatt        <  1:758255/150‑1 (MQ=255)
                                                                                                                                                  gCCGTTAATGTGGGTAAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTcg     <  1:1403933/150‑1 (MQ=255)
                                                                                                                                                   ccGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTcgc    >  2:745235/1‑150 (MQ=255)
                                                                                                                                                   ccGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTcgc    >  1:355121/1‑150 (MQ=255)
                                                                                                                                                     gTTTAGTTGGGTTAAATGTTTCCTTTATGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTcgccg  <  2:1144898/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCG  >  NC_000913/2560847‑2561145

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: