Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,560,996 G→T 21.6% 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,9960GT21.6% 112.8 / 7.7 70intergenic (+98/‑372)yffL/yffMCPZ‑55 prophage; uncharacterized protein/CPZ‑55 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (24/28);  new base T (7/10);  total (31/39)
Fisher's exact test for biased strand distribution p-value = 7.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGC  >  NC_000913/2560849‑2561143
                                                                                                                                                   |                                                                                                                                                   
gggAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGTAATAACGGAATCAGGCGGGCAATGTTGCGGCCTTTTCTTTGACTTATTTCCATGAATTCAGGTGGGTCAAAAGTTGCCCTTTGTGGTGGNTCAAAAGNTGCCttt                                                                                                                                                   >  1:1668906/1‑150 (MQ=255)
gggAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTAAAAATTTGACGtt                                                                                                                                                   >  1:1677976/1‑150 (MQ=255)
 ggAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGGGGGTCAAAAGTTGCCGTTa                                                                                                                                                  <  2:1131941/150‑1 (MQ=255)
 ggAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGATCCCGCCGTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTa                                                                                                                                                  <  2:66225/150‑1 (MQ=255)
 ggAGTGCAAAGAATACCAGTTATTAATCAAGCCGTACAACATGAACTAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGACTTCAGGTGGGTCAAAAGTTGCCGTTAGTGGGGGGTCAAAAGTTGCCGTTa                                                                                                                                                  <  2:1178756/150‑1 (MQ=255)
    gTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAgg                                                                                                                                               <  2:776319/150‑1 (MQ=255)
      gCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTATTGTTTGGTAAAAATTTGCCTTTaatgag                                                                                                                                             >  2:269156/1‑146 (MQ=255)
        aaCGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGGCAAAAGTTGCCGGTAGTGGTGGGTCAAAAGGTGCCGTTAAgggggg                                                                                                                                           <  2:1973263/150‑1 (MQ=255)
        aaCGAATACCAGTTATAATTCAAGCCGTACAACATGAAATAACGGAATCACCCGGGCAATGTTCCCGCCTTTTGTTTGCCTTAATTCCATGACTTCAGGTGGGTCAAAAGTTGCCGTTAGTGTTGGGTAAAAAGTTGCCGTTAAGGTggg                                                                                                                                           <  2:306783/150‑1 (MQ=255)
          cGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTAGACGTTAAGgtgtgta                                                                                                                                         >  2:1022898/1‑149 (MQ=255)
              tACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCaaaa                                                                                                                                     <  2:2070679/150‑1 (MQ=255)
              tACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCaaaa                                                                                                                                     <  1:1882480/150‑1 (MQ=255)
                   gTTATTATTCAAGCCGTACAACATGAAATAACGGAATCCGGCGTTCAATTTTACCTCCTTTTCTTTGCCTAAATTCCATGAATTCAGGTGGGTCAAAATTTGACTTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGc                                                                                                                                <  2:2176801/150‑1 (MQ=255)
                     tattatTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCg                                                                                                                              <  2:1041847/150‑1 (MQ=255)
                       ttattCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGtt                                                                                                                            <  1:780142/150‑1 (MQ=255)
                            cAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTTCTGGTGGGTCAAATGTTGCCGTATAGGAGGGTCAAAAGTTTACGGTAAggg                                                                                                                       >  2:708450/1‑149 (MQ=255)
                            cAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCAGCCTTTTCTTTGACTTAATACCATGAATTCACGTGGGTCAAAAGTTGCCGTAAGTGGTGGGTCAAAAGTTGCCTATAAGTTGGGTCAAAAGTTGCAGTTAAGGt                                                                                                                       >  1:381228/1‑150 (MQ=255)
                                 cGTACAACATGAAATAACGGAAACAGGCGTGCAATTTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGGGGGTCAAAAGATGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTc                                                                                                                  <  1:177380/150‑1 (MQ=38)
                                     caacaTGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAATTTGCCGTGAAGGTGGGTCaaat                                                                                                              >  2:1916669/1‑149 (MQ=255)
                                          tGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGACGTTAAGGTGTCTCAAAAGTTTACGATAACGTGTGTCaaaattaga                                                                                                         >  2:380102/1‑145 (MQ=255)
                                           gAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGcc                                                                                                        <  1:750791/150‑1 (MQ=255)
                                           gAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAATGGTGGGTAAAAATTTGCCTTTAAGGTGTGTCAAAATTTGAATTTAAGGTGTGTCTAAATTTgaa                                                                                                        >  2:1840783/1‑148 (MQ=255)
                                            aaaTAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTATTGGTGGGTCAAAAGTTGCATTTAAGGAGAGTCAAAATTTGAAGTTAACGTGGCGCAAAATTTgaag                                                                                                       >  2:2176023/1‑147 (MQ=255)
                                              ttAACTGAATCAGGCGGGCAATGTTGCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAg                                                                                                     <  1:1540043/149‑1 (MQ=255)
                                              aTAACGGAATCAGGCGGGCAATGTTCCAGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCTTTAAGGAGTGTCAAAAGTTACCGTTAAGGTGGCTCAAAAGTTgaaaag                                                                                                     >  2:1492644/1‑145 (MQ=38)
                                               tAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGTGTCAAAAGTTGCCGTTAAGGTGGGTCAAAATTTGACGAGt                                                                                                    >  2:1909441/1‑150 (MQ=255)
                                                  cGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGTTGGGTCAAAAGTTGCCGTTTAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAATTTGCCGAGTAAg                                                                                                 >  2:795137/1‑150 (MQ=255)
                                                  cGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAATGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAg                                                                                                 <  1:1723439/150‑1 (MQ=255)
                                                       acAGTCGTGCACTCTTCTCGCCTTGATTTTGCCTTAATTACATGAAGTCAGGATGGTACAAAGTTGCCCTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGGCCAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTc                                                                                            <  1:294268/149‑1 (MQ=21)
                                                         aGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGGCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAATGTTGCCGAGTAAGAATTcta                                                                                          >  1:1750503/1‑150 (MQ=255)
                                                          gtcGGGCAATGATCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGTTGGGTCAAAAGTAGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTctac                                                                                         <  1:1835102/148‑1 (MQ=255)
                                                          ggCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAATGGTGGGTAAAAATTTGCCTTTAAGGAGCGTCAAAATTTGAAGAAAAGGTGGCGCAAAATTTGCCGCGAAAGAATTgtaa                                                                                         >  2:1618148/1‑149 (MQ=25)
                                                              ggCAATCATCCCGCCTTTTCTTTTCCATAATTCCATGAATTCCGGTGGGTCAAAAGTAGGCGTTAGTGGGGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTAtt                                                                                     <  1:2159610/150‑1 (MQ=37)
                                                                 aaTGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTATTGGTGGGTCAAAATTTGCCGTTAAGGAGGGTCAAAATTTGAAGTTAAGGTGGGTCAAAATTTGAAGAGTAAGAATTATACTATTaaa                                                                                  >  2:1850021/1‑149 (MQ=37)
                                                                    gTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGGCGTTATGGTGGGTCTAAAGTTGCCGAGTAAGAATTCTACTTTTAAGAAt                                                                               >  1:1157992/1‑150 (MQ=255)
                                                                      tcaagacTTTTTTTTTCCATAAATAAATGAATTCATGTGGGTGAATAGTTGGCGTTTATTGTGGGTAATAAGTTGCCTTTAAGTTGGGTAAAAAGTTTCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTc                                                                             <  1:1295425/144‑1 (MQ=19)
                                                                        ccGTCTTTTATTTTCCTTAATTACATGAATTCATTGTTCTTAAAAGTTGCCGTTAGTGGGGGGTATAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTa                                                                           <  2:1994331/150‑1 (MQ=25)
                                                                         cGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTTAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAg                                                                          >  1:503314/1‑150 (MQ=255)
                                                                                 cTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGTTGGGTCAAAATTAGCCTTAAAGGAGTGTCAAAATTTGACGTAAAGGTGGGTCAAAAGTTGACTCGTAAGAATTGTAAAATAAATAATTATATGCGACggc                                                                  >  2:1995224/1‑149 (MQ=255)
                                                                                    tGCCTTAAATCCATGAATACCGTCTGGTCAAAATTTCCCTTTCGTGGGGGGTCAAAAGCCGCCGTTTAGGTGGCTAAAAAGTTGCCCTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGt                                                               <  1:2126743/150‑1 (MQ=25)
                                                                                               cATGAATTCAGGTGGGTCAAAAGTTGCCGATAGTGGAGGGTAAAAAGTTGCCGTTAAGGAGTGTGAAAATTTGAAGAAAAGCGGGGGGAAAATTTGACTAGTAAGAATTAAACAAGTAAGAATACTAGGCGACGGGATTTTTTTAgagtt                                                    >  2:1939208/1‑146 (MQ=255)
                                                                                                 ggAATTCCGGTGGGTCAAAAGTCGGCGTTAGTGGTGGCTCAATAGTAGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTccc                                                  <  1:1150031/149‑1 (MQ=38)
                                                                                                       ctggggggTAAAAAGGAGCCCTTTGTGGTGGGTCAAATGTAGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCcggtcgg                                            <  1:909449/145‑1 (MQ=37)
                                                                                                            gggTCAAAAGTCGCCGTTATTGGGGGGTAAAAAGCTGCCTTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATa                                       <  1:846233/150‑1 (MQ=38)
                                                                                                              gTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTATAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGAAGAGTAAGAATTCTACTATTAAGAATTCTATGCTACGGGAGATTTTTATCCTAAAGGTCGGTCATaga                                     >  2:1428873/1‑150 (MQ=25)
                                                                                                                 aaaaGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGTTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATagagag                                  <  2:2010035/150‑1 (MQ=255)
                                                                                                                      ttGCCGTTAGTGGGGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGcc                             <  1:1022898/150‑1 (MQ=255)
                                                                                                                       tggcGTTATTGTTGGGTCAAAAGTTGCCTTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCt                            <  1:1994398/147‑1 (MQ=255)
                                                                                                                       tGCCTTTAGTGGTGGGTAAAAACATGCCGTTAAGTTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCt                            <  2:1663097/150‑1 (MQ=38)
                                                                                                                        gCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTa                           >  1:1721888/1‑150 (MQ=255)
                                                                                                                           gTTAGTGGGGGCGCAAAACATGCCCTTAAGTTCGGTCAAAAGTTGCCGTTAAGGGGGGGCAAAAGTTGCCGAGTAAGAATTCTACAATTAAGACTTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTagag                        <  1:69311/150‑1 (MQ=255)
                                                                                                                             tAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGTTCGGTCATAGAGAGTCGCCTAGAGtt                      >  2:1855948/1‑150 (MQ=255)
                                                                                                                              atttgtggGTCAAAACTTGCCTTTTAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTa                     <  2:1126872/146‑1 (MQ=255)
                                                                                                                              aGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGTGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTa                     <  1:995642/150‑1 (MQ=255)
                                                                                                                                 gtggggTCAATACTAGCCGTTACGGTGGTGAAATCGTCGCCGTTCAGTTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCAAGGCGACGGGATTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAAc                  <  2:1844096/147‑1 (MQ=255)
                                                                                                                                 ggtggGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGATTTTTTTAGCCTCCCGGTCGGTCTTAGAGAGTCGCCTAGAGTTATAc                  >  1:325089/1‑150 (MQ=255)
                                                                                                                                 ggtggGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGCAGTATTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAAc                  >  2:1310903/1‑150 (MQ=255)
                                                                                                                                    gggTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAAGTCGACTAGTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGGCGGTCATAGAGAGTCGCCTAGAGTTAAACcga               >  1:93008/1‑150 (MQ=38)
                                                                                                                                     ggTCAAAAGTCGCCTTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACagaa              <  1:874426/150‑1 (MQ=255)
                                                                                                                                      gTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACagaag             >  1:126809/1‑150 (MQ=255)
                                                                                                                                        cAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAgaagaa           >  2:488788/1‑150 (MQ=255)
                                                                                                                                        cAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAgaagaa           >  1:1652055/1‑150 (MQ=255)
                                                                                                                                          tttgTTTCCTTTTAGTTGGGTCAATTGTTGCCTTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAatt         <  2:159384/147‑1 (MQ=255)
                                                                                                                                           atgTTGCCTTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAatta        <  2:1239824/148‑1 (MQ=255)
                                                                                                                                           aaGGTGCCTTTAAGGTGGTTCAAAAGTAGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAATAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAatta        <  1:1065583/150‑1 (MQ=39)
                                                                                                                                            atttGCCTTTAAGTTGGGTCAAAAGTTGCCTTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAattat       <  2:1650567/148‑1 (MQ=255)
                                                                                                                                               tGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTc    >  1:574521/1‑150 (MQ=255)
                                                                                                                                                gacGTTAAGTTGGGTAAAAAGTTGCCTTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTcg   <  1:1388866/148‑1 (MQ=255)
                                                                                                                                                gCCTTTACGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTcg   <  1:251607/150‑1 (MQ=255)
                                                                                                                                                 ccGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTcgc  >  1:1578287/1‑150 (MQ=255)
                                                                                                                                                 ccGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTcgc  >  1:1162617/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGC  >  NC_000913/2560849‑2561143

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: