Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,560,985 C→A 17.0% 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,9850CA17.0% 107.3 / 3.9 62intergenic (+87/‑383)yffL/yffMCPZ‑55 prophage; uncharacterized protein/CPZ‑55 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (29/21);  new base A (7/5);  total (36/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATCGGCAGAAAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAA  >  NC_000913/2560837‑2561134
                                                                                                                                                    |                                                                                                                                                     
gATCGGCAGAAAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGTGCAATGTTCACGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGTTCa                                                                                                                                                      >  1:280524/1‑150 (MQ=255)
    ggCAGAAAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACTTGAAATAACGGAATCAGGCGGGCATGGTTCCCGCCGTTTCTTTGCCTTAATTCCATGAATTCAGGGGGGTCAAAAGTTGCCGTTATTGGTGGGTCAAAAg                                                                                                                                                  >  1:757719/1‑150 (MQ=255)
    ggCAGAAAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCaaaat                                                                                                                                                  >  1:430316/1‑149 (MQ=255)
     gCAGAAAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGt                                                                                                                                                 <  1:1714162/150‑1 (MQ=255)
     gCAGAAAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGt                                                                                                                                                 >  2:378774/1‑150 (MQ=255)
        gAAAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGATATTGATTGGTaaaattttga                                                                                                                                              >  1:32532/1‑144 (MQ=255)
          aaGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCg                                                                                                                                            <  2:845859/150‑1 (MQ=255)
           aGGTAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATGCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGt                                                                                                                                           <  1:1411986/150‑1 (MQ=255)
            gggAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGtt                                                                                                                                          <  2:1861507/150‑1 (MQ=255)
            gggAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGtt                                                                                                                                          <  1:682426/150‑1 (MQ=255)
                 tGCAACGAATACCTGTTATTATTCAAGCCGTACAACATCAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTTATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGATAGTGGAGGGTAAAAAGTTGCCGTTAAGGt                                                                                                                                     >  2:1535631/1‑150 (MQ=255)
                  gCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTg                                                                                                                                    >  1:1704698/1‑150 (MQ=255)
                    aaCGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTggg                                                                                                                                  <  2:522561/150‑1 (MQ=255)
                    aaCGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTggg                                                                                                                                  >  1:1895242/1‑150 (MQ=255)
                          tACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCaaaa                                                                                                                            <  2:1740816/150‑1 (MQ=255)
                           aCCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGACGTTAAGGTGGGTCaaaat                                                                                                                           >  1:849803/1‑149 (MQ=255)
                           aCCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGATAGTGTTGGGTCAAAATTTACCTTTAACTTTAGTCaaaat                                                                                                                           >  1:382851/1‑149 (MQ=255)
                              aGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCACGACTGTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAATTTTGAGGTAAAGGAGGGTAAAAAGTTg                                                                                                                        >  1:1044642/1‑150 (MQ=255)
                               gTTATTATTCAAGCCGTACAACATGAACTAACGGACTCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTAAGTGGGGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGc                                                                                                                       <  1:905537/150‑1 (MQ=255)
                               gTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAATTTGACGTTAAGGTGGGTCAAAAGTTga                                                                                                                       >  1:504497/1‑149 (MQ=255)
                               gTTAGTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGTGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGc                                                                                                                       <  1:1188391/150‑1 (MQ=255)
                                 taatattCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTATTCTATGCCTTAAATCCATGAATTCAGGTGGGGGAAAAGTTGCCGTTAGTGGGGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCg                                                                                                                     <  2:1445972/147‑1 (MQ=255)
                                  attattCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCATTTTCTTTGCCTTAATTCCATGAAATCAGGTGGGTCAAAAGTTGCCGTTATTGGTGGGTAAAAATTTGCAGTTAATGAGGCTCAAAAGTTgacta                                                                                                                    >  2:1443936/1‑146 (MQ=255)
                                        cAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAATTTGACTTAAAGGAGTGTAAAAAGTTGACGTTAAgct                                                                                                              >  2:497048/1‑148 (MQ=255)
                                        cAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGt                                                                                                              >  2:1902836/1‑150 (MQ=255)
                                        aaaGCCGTACAACATGAAATTACGGAATCAGGCGGGCAATGTTCCCGCATTTTCTTTGCCTTAATTCCAGGAATTCAGGTGGGTAAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGt                                                                                                              <  2:1384809/149‑1 (MQ=255)
                                         aaGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTg                                                                                                             <  1:1505610/150‑1 (MQ=255)
                                                 caacaTGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGTTGGGTAAAAAGTTGCCGTTAAGGTGGGTGAAAAGTTGCCGTTAACGTGGCTCaaaa                                                                                                     >  2:759839/1‑150 (MQ=255)
                                                   acaTGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCATAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTGaaaatt                                                                                                   >  1:61510/1‑148 (MQ=255)
                                                     attaaaTAACGGAATCAGTCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTAAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTg                                                                                                 <  1:1702067/147‑1 (MQ=255)
                                                     aTGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTTGTGGGTCAAAAGTTGCCTTTAAGGTGGGTCAAAAGTTGCCGTTAAGTTGGGTCAAAAGTTg                                                                                                 >  1:487161/1‑150 (MQ=255)
                                                            aaCGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTTGGTAAAAAGTTGCCGTTTAGGTGGCTCAAAAGTTGACGTTTACGTGGGGGTAATGTTGCCTAGTa                                                                                          >  1:1529513/1‑150 (MQ=255)
                                                                   accgtcGGGCAATTTTCCAGCCATTTCTTTGCCTTATTTCCATGTATTCAGGTGGGTCAAAAGTTGCCGTTACTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTc                                                                                   <  2:1167673/145‑1 (MQ=37)
                                                                    cAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAATTTGCCGTTAAGGTGGGTCAAAAGTTGCCGATAAGGTGGGTCAAAAGTTGACGAGTAAGAAGTct                                                                                  >  2:1598592/1‑150 (MQ=255)
                                                                    cAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGAGTCAAAAGTTGCAGTTAAGGTGGGTCAAAAGTTGACTAGTAAGAAttgt                                                                                  >  1:378417/1‑148 (MQ=255)
                                                                            cAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTTCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTaa                                                                          >  1:821277/1‑150 (MQ=255)
                                                                                 ttCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAAta                                                                     >  2:1885965/1‑149 (MQ=255)
                                                                                      gCCTTTTCTTTTCCTTAATTACATGAATTCAGGTGGGTCAAAATTTGGCGTTATTGGTGGGTAAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAgg                                                                <  2:644178/150‑1 (MQ=37)
                                                                                      gCCTTTCGTTTTCCTTAAATTCATGAATTCAGTGGGGTAAAAAGTCGCCGTTCCTGGTGGCTCAAACGCTGCCGTTACGGTGGGTCAAAAGTTGCCTTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAgg                                                                <  2:959342/150‑1 (MQ=25)
                                                                                        cTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCTAAAGTTGCCGTTAGTGGTGGGTAAAAAGTTGCCTATAAGGTGGGTCAAAAGTTGCCTTTAAGGTGTGTCAAAAGTTGCCGAGAAAGAATTCTACTATTAAGAATTCTAGGCg                                                              >  1:270729/1‑150 (MQ=255)
                                                                                          tttCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTATTGGTGGCTCAAAAGTTGCCGTTAAGGAGGGTCAAAAGTTTCCGTTAAGGTGGGTCAAAATTTGACGAGAAAAAATTCTACTATAAAGAATTATAGGCGAc                                                            >  2:1406167/1‑150 (MQ=37)
                                                                                            tCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAATGGTGGGTAAAAAGTTGCCTTGAAGGAGTGTCTAAAGTTGACGCGAAGCTGGCTCACATTTAGCGGAGTATGAATTCTAAAATTAAGAATTCTAAGCGAcat                                                          >  1:1692725/1‑148 (MQ=19)
                                                                                                  ccTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTTGTGGTGGGTAAAAAGTTGCCTTCAAGGAGGCTCAATAGTTGCCGTTAAGGTGGGTCTAATTTTGCCGAGAATGAATTATACTATTAAGAATTCTAGGCGTCGTGAGttt                                                    >  1:1688970/1‑150 (MQ=37)
                                                                                                    ttAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTAAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTCCCGAGTAAGAATACTACTATTAAGAATAATACGCGACGGCAGttttt                                                  >  1:1615316/1‑150 (MQ=37)
                                                                                                     tAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTTAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTAAAAAGTATCCGAGTAAGAATTATACTATTAAGAATTATACGCGAAGGTAGtttttt                                                 >  1:1740191/1‑150 (MQ=37)
                                                                                                     tAATTCCATGAATTCAGGGGGGTCAAAAGATGCCGTTAGTGGTGGGTCAAAACTTGCCGTTACGGTGGGTCAAAAGTTGCCGTTAAGGGGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGtttttt                                                 <  1:712521/150‑1 (MQ=38)
                                                                                                      aaatcaaTGAATTCAGGTTTGTTAAAATCTGCCTTTAGGGGTGGGTAATTAGCTGCCTTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTa                                                <  2:530340/144‑1 (MQ=25)
                                                                                                      aaatCCAGGACTTCATGCGGGGCAAAAGCTGCCGTTCGTGGTGGGTCAAAAGTGGCCGTTAAGGTGGGGCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTa                                                <  2:739845/147‑1 (MQ=37)
                                                                                                        ttCCATGAATTCAGGTGCGTAAAAACATGCACTTATTGGGGGCTCTAATGCCTCCTTTCAGGTGGGTCAAAAGTTGCCTTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGc                                              <  1:1712864/150‑1 (MQ=25)
                                                                                                         tttttgtaTTTAGTTGTGTAAAATGTTGGCTTTTTTGGTGGTTAAAAATTTTACTTTAAGTTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGcc                                             <  1:1636122/143‑1 (MQ=21)
                                                                                                          ccATGCATTCAGGGGGGTCAAAAGATGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAACGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCt                                            <  1:456671/150‑1 (MQ=255)
                                                                                                          ccATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGAGGGTCAAAAGTTGCCGTTAACGTGGGTAAAAATTTGAATAGTAAGAATTCTAATATTTTGAATAATAGTAGACGGGAGTATTTTataat                                            >  2:1781314/1‑146 (MQ=21)
                                                                                                             tGAATTCAGTTGGGTAAAAATTTGCCGTTATTGGTGGGTAAAATGTTGCCGTTTAGGTGGGTGAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTccc                                         <  2:435116/150‑1 (MQ=38)
                                                                                                                   cAGTTGTGTCAAAAGTAGGCGTTATTGGGGGGTAAAAAGTTGCCGTTTAGGTGGGTCAAAAGTTTCCGTTAAGGTGGGTCAAAAGTTGCCGAGTATGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCcggtcgg                                   <  2:1758086/150‑1 (MQ=37)
                                                                                                                     ggTGGGCCAAATGCTGGCGTTTGTGGTGGGTCAAAACTGGCCGTTAAGGTGGGTCAAAAGGTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCggtcggtc                                 <  1:1490505/150‑1 (MQ=37)
                                                                                                                       tGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTGTAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAATAATTCTAGGCGACGGGAGTTTTTTTTACTCCCGGTCGGTAAt                               >  1:803873/1‑150 (MQ=37)
                                                                                                                                    gCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTa                  >  1:597712/1‑150 (MQ=255)
                                                                                                                                      cGTTAGTGGTGGGTCAAAAGTAGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTaga                <  2:1485663/150‑1 (MQ=255)
                                                                                                                                         gaGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCCAGTAAGAATTCTACTATTAAGAATACTATGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCAAGAGtt             >  1:1336072/2‑150 (MQ=255)
                                                                                                                                             ggtggGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAAc         >  1:1210175/1‑150 (MQ=255)
                                                                                                                                               tggGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACat       >  2:684124/1‑149 (MQ=255)
                                                                                                                                                  gTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGATTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGATATAAACagaag    >  2:1211343/1‑150 (MQ=255)
                                                                                                                                                  gTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACagaag    >  2:423207/1‑150 (MQ=255)
                                                                                                                                                    cAAAAGTTGCCGGTAAGGTGGGTCAAAAGTTGCCGTTTAGGTGGGTCAAAAGTTGCCCAGTAAGAATTCTACTATTAATAATTATAGGCCACGGGAGTTGTTTAGCCTCCCGGTCGTTCATATATAATCGCCTAGAGGTTAACAgaagaa  >  2:496410/1‑150 (MQ=25)
                                                                                                                                                    |                                                                                                                                                     
GATCGGCAGAAAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAA  >  NC_000913/2560837‑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: ‑