Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,560,985 C→A 11.8% 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,9850CA11.8% 149.5 / 2.6 81intergenic (+87/‑383)yffL/yffMCPZ‑55 prophage; uncharacterized protein/CPZ‑55 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (32/36);  new base A (9/3);  total (41/40)
Fisher's exact test for biased strand distribution p-value = 1.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGGCAGAAAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAA  >  NC_000913/2560839‑2561134
                                                                                                                                                  |                                                                                                                                                     
tCGGCAGAAAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCaaa                                                                                                                                                    >  1:1780112/1‑150 (MQ=255)
 cGGCAGAAAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCaaaa                                                                                                                                                   <  1:434749/150‑1 (MQ=255)
  ggCAGAAAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAg                                                                                                                                                  <  2:1243879/150‑1 (MQ=255)
   gCAGAAAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTGTGTGTTGTGTCaaaatt                                                                                                                                                 >  2:1594724/1‑148 (MQ=255)
       aaaGTGAGTGCAACGTATACCAGTTATTATTCAAGCCGTACTACATGAAATTACGGAATCATTCGGGCAATGTTCCCGCCTTTTCTTTGCCTTACTTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGcc                                                                                                                                             <  1:593031/150‑1 (MQ=255)
        aaGGGAGTGCAACGAATACCAGTTTTTATTCAAGCCGTACAACATGTAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGACTTAATTCAATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCg                                                                                                                                            >  1:1347748/1‑150 (MQ=255)
            gAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGTATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTaa                                                                                                                                        <  2:826853/150‑1 (MQ=255)
               tGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATACCGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGATGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGt                                                                                                                                     <  2:1509580/150‑1 (MQ=255)
               tGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGt                                                                                                                                     <  2:727206/150‑1 (MQ=255)
                gCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTg                                                                                                                                    <  1:1582323/150‑1 (MQ=255)
                gCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTACATGAATTCAGGTGGGTCAAAAGTTGCCGTTATTGGAGGGTAATAATTTGCCGTTAAggat                                                                                                                                    >  1:181664/1‑148 (MQ=255)
                gCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTATCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTTGAGGGTCAAAAGTTGCAGTTAAGGTg                                                                                                                                    >  1:374293/1‑150 (MQ=255)
                 cAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAATTTGACTTTAAGGTgg                                                                                                                                   >  1:444446/1‑150 (MQ=255)
                          ccAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAATGTTGTGTCAAAATTTGCCTTGTAGGTTGGTCAATAGt                                                                                                                          >  1:92418/1‑150 (MQ=255)
                          ccAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGAGGGCAATGTTCCCGCCTTTTCTATGCCTTAATTACATTAATTCAGGTGGGTCAAAAGTTTCAGTTAGTGGATGGTCAAAAGTTGAATTTAAGGAGAGTCAAATGt                                                                                                                          >  1:531652/1‑150 (MQ=255)
                           cAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCTTGAATTCAGGTGGGTCAAAAGTTGCCGTTATTGCTGGTTAAAAAGTTGCCTTAAAGGAGTGGCAAAAGtt                                                                                                                         >  2:996466/1‑150 (MQ=255)
                           cAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGtt                                                                                                                         >  1:537580/1‑150 (MQ=255)
                           cAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGATTTTGGTGGGTAAAAATTTCCCTTTAAGGGGGGGAAAAAGtt                                                                                                                         >  1:760742/1‑150 (MQ=255)
                            aGTTATTATTCAAGCAGTACAACATGAAATAACTGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCTGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTg                                                                                                                        <  1:723703/150‑1 (MQ=255)
                               tattatTCAAGCCGTACAACATGAAATAACGGAATCAGGCGTGCAATGTTCCCGCCTTTGCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCg                                                                                                                     <  1:1282419/150‑1 (MQ=255)
                                  tattCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAATTTGCCGTTAAGGTGTGTAAAATGTTTCAGTTa                                                                                                                  >  1:625113/1‑150 (MQ=255)
                                  tattCAAGCCGTACAACATGAAATAACGCAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTa                                                                                                                  <  1:1217781/150‑1 (MQ=255)
                                   attCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTATTGGTGGGTAAAAATTTGACTTTTAGGAGTGTCAAAATTTgaattaaa                                                                                                                 >  2:1219485/1‑143 (MQ=255)
                                      cAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGTTGGGTAAAAATTTGCCGTTAAGGAGTCTAAAAAGATGAAGTTAAGGt                                                                                                              >  2:461370/1‑150 (MQ=255)
                                      cAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGGCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGt                                                                                                              <  1:1024923/150‑1 (MQ=255)
                                       atgCCGTACAACCTGAAATACCGGACTCAGCCGGGCAATTTTCCCGCCTTTTCTTTGACTAAATTCCATTAATTCAGTGGGGTCAAAAGTTGCCGTTTTGGGTGGGTCAACATTTTCCGTTAAGGTGGGGCAAAAGTTGCCGTTAAGGTg                                                                                                             <  2:1347180/148‑1 (MQ=255)
                                       aaGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTg                                                                                                             <  2:1999431/150‑1 (MQ=255)
                                        aGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTgg                                                                                                            <  1:1624027/150‑1 (MQ=255)
                                        aGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGGCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTgg                                                                                                            <  1:1155427/150‑1 (MQ=255)
                                        aGCAGTATATCATGAAATCACGGTATTAGGGGTGCATTGATCCGGCCTTTTGTTTTCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTCGCCGTTAGTGTTGGGTGAATAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTgg                                                                                                            <  1:1461779/150‑1 (MQ=25)
                                         gCCGTACAACATGAAATAACGGAATCAGGCGGGCATTGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTTGTGGTTGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTggg                                                                                                           <  1:1479276/150‑1 (MQ=255)
                                           cGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTc                                                                                                         >  2:1765882/1‑150 (MQ=255)
                                               aaacaTGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTATGCCTTAATTCCATGAATTCAGGTGGGGCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGGTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGGCaaaa                                                                                                     <  1:1578026/149‑1 (MQ=255)
                                                 acaTGAAATAACGGAATCAGGCGGGCTATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATACAGGTGGGTCAAAAGTTGCCGTTTGAGGTGGGTCAAAAGTTGACTTTAAGTTGTGTCAAAATTTTCAGTTAAGGAGGGTGTAAAGt                                                                                                   >  1:270933/1‑150 (MQ=37)
                                                  caTGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTTGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCTAAAGTTGCGGTTTAGGTGGGTCAAAAGtt                                                                                                  >  1:1913498/1‑150 (MQ=255)
                                                     gAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGTGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCAGTTAAGGTCGGACTAAAGtttaa                                                                                               >  2:391097/1‑147 (MQ=255)
                                                     gAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCTTTAAGGTGGGTCAATTGTTGCCGTGTAGGTGGGGCAAATGTTGcg                                                                                               >  1:1506481/1‑149 (MQ=255)
                                                        aTAAAGGAATCTGCCTTGCAATCTTCCCGCCTTATCTTTTCCATAATTCCATGAATTCAGGTGTGTCAATCGCTGGCGTTAGTGGTGGGTCAAAAGTTGCCGTAAAGTTGGGTCAAAAGTTGCCGTTAAGGTGGGGCAAAAGTTGCCGAg                                                                                            <  2:1231375/150‑1 (MQ=21)
                                                           aCGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCAGGTTAGGTGGGGCAAAAGTTGCCGAGTaa                                                                                         >  2:1690451/1‑150 (MQ=255)
                                                            cGGAATCAGGCGGGCAATTTTCCCGCCTTTTATTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGCTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAg                                                                                        <  1:1691111/150‑1 (MQ=255)
                                                               aaTCAGGCGGGCAATGTTACCGCCTTTTCTTTGCCTTAAATCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTACTGGTGGGTAAAAAGTTGCCTTTAAGGTGGGTCAAAAGTTGCAGTTAAGGGGGGTATAATGTTTCCGAGTAAGAAt                                                                                     >  2:775758/1‑150 (MQ=37)
                                                                   cgtctggCAAATTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCCGTTGGGTCAAAAGTTGCCGTTAGTGGGGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTcta                                                                                 <  2:1674613/145‑1 (MQ=37)
                                                                   aGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAGTTGCCGAGTAAGAATTctac                                                                                >  1:696531/1‑150 (MQ=255)
                                                                    ggCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTATTGGTGGGTAAAAAGTTGCCTTTAAGGTGGGTCAAAAGTTGCAGTTTAGGTGGGTCTAAAGTTTACGAGTAAGAATTctaa                                                                                >  1:546606/1‑149 (MQ=37)
                                                                                 cccGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTAAAAATTTGCCTTTTAGGAGGCTCAAAAGTTGAAGAAAAGGGGTCTCAAAATTTGACGAGAAAGGATTCTACAATTAGCAAttgt                                                                   >  2:66729/1‑148 (MQ=19)
                                                                                  ccGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGAGTCAAAAGTTGCCGTTTAGGTGGGTCAAAATTTGCCGAGTAACCATTCTAATATTAAGAATtata                                                                  >  1:1070881/1‑150 (MQ=37)
                                                                                   cGCCTTTTGTTTTCCTTTATTACATGAATTCAGTGTTGTCAAAAGTTGCCTTTATTGGTGGGTAATTTGTTGCCTTTAAGGTGGGTCAAATGTTGCCTTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAg                                                                 <  2:595591/150‑1 (MQ=25)
                                                                                   cGCCTTTTCTTTTCCTTAATTCCATGAATTCAGGTGGGGCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGGCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAg                                                                 <  2:1572488/150‑1 (MQ=255)
                                                                                    gCCTTTTCTTTTCCTTTATTCATTGACTTCAGGTGGGTCTATAGTTGGCTTTACTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAATAGTTGCCGTTAAGGTGTGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAgg                                                                <  2:1572314/150‑1 (MQ=25)
                                                                                    gCCTTTTCTTTGCCTTTATTTCATGAATTCTTGTGGGTAAAAAGTTGGCTTTATTGGTGGGTAAATAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAgg                                                                <  2:1662278/150‑1 (MQ=37)
                                                                                    gCCTTTTCTTTGCCTTAAATCCATGAATTCAGGTGGGACAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAgg                                                                <  1:773226/150‑1 (MQ=255)
                                                                                    gCCGTTTCTTTCCCTTCAATCAATGAATTCAGGGGGGTCAAAAGTTGGCTTTAGTGGTGGGTCAAAAGTTGCAGTTAAGGGGGGTAAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAgg                                                                <  2:1691523/150‑1 (MQ=37)
                                                                                      cTTTTCTATTCCTTTAATTAATCAATTCAGTTTGGTATAAATATGCGTTTATTGGTGGGTCAAAAGTTGCATTTTAGTTGGGTAAAAAGTTGCATTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCg                                                              <  2:486360/150‑1 (MQ=16)
                                                                                          tCTTTGCCTTTATTAAATCAATTCATGTTGGTAAAAAGTTTCGCTTTGTGGTGGGTCAAAAGTTGCCGTTACGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACgg                                                          <  2:677999/150‑1 (MQ=25)
                                                                                             ttGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAATTTGCCGTTAAGGAGGGTCAAAATTTGACGTTAAGCTGGGGCAAAATTTGACTAGAAAGAATTATAAAATAAAGAATAATACGCGACGTGat                                                       >  2:1580863/1‑149 (MQ=21)
                                                                                                 cTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAATTTGCCGTTAAGGTGGGTCAAAAGTTGCAGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGtttt                                                   >  1:1981254/1‑150 (MQ=255)
                                                                                                  ttAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTAAAAAGTTGCATTTAAGGTGTGTCAAAAGTTGACGTTAAGGTGGCGGAAAATTTGCCTAGAAAAGATTGTAACATTAAGAATTCTAGGCGGAGGTAGttttt                                                  >  2:850194/1‑150 (MQ=21)
                                                                                                     attaaaTGAATTCAGGTGGGTAAAAAGTTGGCGTTATTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAg                                               <  2:1608710/145‑1 (MQ=255)
                                                                                                        ccATGAATTCAGGGGGGTCAAAAGTTGCCGTTAGTGGAGGGTCAAAAGATGCCTTAAAGCTTGGTAAAAATTTGCCGTTAAGGTGGGTCAAAAGTAGCCGAGTAAGAATTCTACAATTAAGAATAATAGGCGACGTGATTTATTTAGccg                                            >  1:360695/1‑149 (MQ=255)
                                                                                                         cATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGTGTCAAAATTTGCCGTTAAGGTGGGTCAAAATTTGAAGTAAAGGTGGGGCAAAATTTGAAGAGTAACAAGTATAAAATAAACAATAATAAGCGTCGGGATTTTTTTAGACta                                           >  2:865628/1‑149 (MQ=16)
                                                                                                         cATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCATTTAAGGTGTGTCAAAAGTAGACGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGACta                                           >  1:1732248/1‑149 (MQ=37)
                                                                                                              aTTCAGGTGGGTCAAAAGTTGCCGTTAGTGCTGGGTCAAAAGTTGCCGTGAAGGTGTGTCAAATGTTTACGTTAAAGTGTGTCAAAAGTTGACTAGTAAAAATTATAAGATTAAGAATTATAGTCGACGGTATTTTTTTAGACTAAcggt                                      >  2:1364516/1‑150 (MQ=19)
                                                                                                                   ggTGGGAAAAAAGTTGGCGTTAGGGGTGGGTCAAAACTTGCCGTTAAGGTGGGTCAAATGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCggtcggtc                                 <  1:1078287/150‑1 (MQ=38)
                                                                                                                    gTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCAGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGATTTTTTTAGCCTACCGGTCGGTCa                                >  2:718650/1‑150 (MQ=255)
                                                                                                                      gggTCAAAAGTTGCCGTTAGTGGTGGGTCAAAACTTGCCGTTAAGGTGGGGCAAGAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATa                              <  1:1920665/150‑1 (MQ=255)
                                                                                                                         tCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATagat                           >  2:1890865/1‑149 (MQ=255)
                                                                                                                           aaaaGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACTGGAGTTTTTTAGCCTCCAGGTCGGTCATagagag                         >  1:1668099/1‑150 (MQ=255)
                                                                                                                               gTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTAAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTTCCGAGTAAGAATTCTACTATTAACAATTCTAGGCGACGGTACTTTTTTAGCCTCACGGTCGGTCATAGAGAGTAGc                     >  1:1564788/1‑150 (MQ=37)
                                                                                                                                tagCCGTTAGCGGTGGGTCAAAAGCTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGcc                    <  1:480882/148‑1 (MQ=255)
                                                                                                                                gtGCCGTTCCTGTTGGGTCAAAAGTTGCCTTTAAGGTGTTTCAAAAGTTGGCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTATGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGcc                    <  2:539163/149‑1 (MQ=38)
                                                                                                                                 tGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCt                   >  2:1588250/1‑150 (MQ=255)
                                                                                                                                 tGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAACAATTCTACTATTAATAATTCTAGGCGACGGGAGTTTTTAAGCCTACCGGTCGGTAATAGAGAGTCGCCt                   >  2:1909662/1‑150 (MQ=37)
                                                                                                                                   ccGTTAGTGGTGGGTCAAAAGGAGCCGTTAAGGGGGGTCAAAAGTTGCCGTTAAGGGGTGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTag                 <  1:1831079/150‑1 (MQ=39)
                                                                                                                                     gTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTagag               <  1:483824/150‑1 (MQ=255)
                                                                                                                                        atgggtggGTAATTAGTTTCCTTTAAGTTGGGTCAATAGTTGCCTTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTa            <  2:546606/147‑1 (MQ=255)
                                                                                                                                             tggGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACag       >  1:1475236/1‑150 (MQ=255)
                                                                                                                                             ggggTCAAAAGCAGCCGTTAAGGTGGGTCAAAAGATGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACag       <  1:1399100/149‑1 (MQ=255)
                                                                                                                                              gggTCAAAACTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACaga      <  1:429082/150‑1 (MQ=255)
                                                                                                                                               ggTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACagaa     >  1:1293760/1‑150 (MQ=255)
                                                                                                                                               ggTCAAAACTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACagaa     <  1:850195/150‑1 (MQ=255)
                                                                                                                                                gTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACagaag    >  1:968061/1‑150 (MQ=255)
                                                                                                                                                  cAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAgaagaa  >  1:1256762/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TCGGCAGAAAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAA  >  NC_000913/2560839‑2561134

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