Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 292,302 C→G 7.0% intergenic (+71/+20) yagJ → / ← yagK CP4‑6 prophage; uncharacterized protein;Phage or Prophage Related/CP4‑6 prophage; conserved protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913292,3020CG 7.0% 411.7 / 17.3 129intergenic (+71/+20)yagJ/yagKCP4‑6 prophage; uncharacterized protein;Phage or Prophage Related/CP4‑6 prophage; conserved protein
Reads supporting (aligned to +/- strand):  ref base C (60/60);  new base G (4/5);  total (64/65)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.61e-01

Only 100 of 131 total aligned reads displayed.

ATGTGCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTG  >  NC_000913/292155‑292451
                                                                                                                                                    |                                                                                                                                                     
ttgtgCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAAAGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATAGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATCAGGCCTTCTCACTg                                                                                                                                                     <  1:1042766/149‑1 (MQ=255)
aTGTGCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTTGAAACAATACAAATAAGGCCTTCTCACt                                                                                                                                                      <  2:394648/150‑1 (MQ=255)
aTGTGCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTg                                                                                                                                                     <  2:1631783/150‑1 (MQ=255)
aTGTGCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTg                                                                                                                                                     >  1:680804/1‑150 (MQ=255)
  gtgCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATACGGCCTTCGCACTgag                                                                                                                                                   >  2:887140/1‑150 (MQ=255)
  gtgCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTgag                                                                                                                                                   >  1:1073098/1‑150 (MQ=255)
  gtgCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTgag                                                                                                                                                   <  1:686427/150‑1 (MQ=255)
    gCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGaa                                                                                                                                                 >  1:622839/1‑150 (MQ=255)
    gCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGaa                                                                                                                                                 >  1:455466/1‑150 (MQ=255)
    gCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAAAGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTTCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACATTACAAATAAGGCCTTCTCACTGAGaa                                                                                                                                                 <  1:1209575/150‑1 (MQ=255)
     cAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAg                                                                                                                                                <  1:1683921/150‑1 (MQ=255)
     cAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGACTTCTCACTGAGAAg                                                                                                                                                >  2:583941/1‑150 (MQ=255)
       ggATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGc                                                                                                                                              >  2:1329277/1‑150 (MQ=255)
        gATGCTGTTTCTACCGCTTAAAATGCAGAAACGATATTTAGTAGCGGACATGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCAGTGAGAAGGcc                                                                                                                                             <  2:90376/150‑1 (MQ=255)
         aTGCTGTTGCTACCGCTGAAAATGCATAAACTTTATTTAGTAGAGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCAGTGAGAAGGcct                                                                                                                                            <  1:926056/150‑2 (MQ=255)
         aTGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCAGTGAGAAGGcct                                                                                                                                            <  1:1326204/150‑2 (MQ=255)
         aTGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCa                                                                                                                                            >  1:653126/1‑150 (MQ=255)
          tGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACTAAACTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCAGTGAGAAGGcctt                                                                                                                                           <  2:941240/150‑3 (MQ=255)
           gCTGTTGCTACCGCTGAAAATGCAGAAACTTTATTTAGTAGCGTACTTGATGATGAATATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGATACAATACAAATAAGGCCTTCTCACTGAGAAGGCCAtt                                                                                                                                          <  1:254152/150‑1 (MQ=255)
           gCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCAtt                                                                                                                                          <  2:847178/150‑1 (MQ=255)
             tGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTTTCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTAc                                                                                                                                        <  1:555403/150‑1 (MQ=255)
                tGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGa                                                                                                                                     <  1:171770/150‑1 (MQ=255)
                 gCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGTACTTGCTGATGACTATCAGTAATCAATAGAAAATTGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGATAAGGCCATTACCGAc                                                                                                                                    <  2:864294/150‑1 (MQ=255)
                   tACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACtt                                                                                                                                  >  2:799403/1‑150 (MQ=255)
                   aaCCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGAAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACtt                                                                                                                                  <  2:846805/149‑1 (MQ=255)
                     ccGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTAc                                                                                                                                >  2:819212/1‑150 (MQ=255)
                       gCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAg                                                                                                                              >  2:1541671/1‑150 (MQ=255)
                       gCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCAATTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAg                                                                                                                              >  1:931699/1‑150 (MQ=255)
                        cTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGt                                                                                                                             <  1:1642196/150‑1 (MQ=255)
                         tGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGtt                                                                                                                            <  1:571250/150‑1 (MQ=255)
                           aaaaTGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCg                                                                                                                          >  2:1621158/1‑150 (MQ=255)
                           aaaaTGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACAAAACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGATCg                                                                                                                          <  2:1067298/150‑1 (MQ=255)
                                 cAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTgg                                                                                                                    >  1:1541843/1‑150 (MQ=255)
                                    aaaCGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCAGTg                                                                                                                                                     >  2:419740/1‑114 (MQ=255)
                                         ttATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTATTAA‑TTTTTTGAAACAATACAAATTATGCCTTCTCACTGAGTAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGcc                                                                                                            <  2:1013605/150‑1 (MQ=255)
                                             ttAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAg                                                                                                        <  2:511450/150‑1 (MQ=255)
                                                  agCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTAACAGTTCGATTTGGCTACAGCCGAAGTTACg                                                                                                   >  1:1047287/1‑150 (MQ=255)
                                                   gCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACgg                                                                                                  <  2:373594/150‑1 (MQ=255)
                                                      gACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAGACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGtct                                                                                               <  1:1296252/150‑1 (MQ=255)
                                                       aCTTGCTGCTGACTATCAGTAATCTATAGACAATTGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACACTTCGATTTGGCTACAGACGAAGTTACGGtctc                                                                                              >  1:688905/1‑150 (MQ=255)
                                                        cTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGTCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTcc                                                                                             >  2:146296/1‑150 (MQ=255)
                                                        cTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTcc                                                                                             <  1:411645/150‑1 (MQ=255)
                                                        cTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTAAAGTTCGATTTGGCTACATCCGAAGTTACTGTCTcc                                                                                             >  2:1286722/1‑150 (MQ=255)
                                                          tGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCtt                                                                                           >  2:883575/1‑150 (MQ=255)
                                                           gCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTc                                                                                          <  1:750551/150‑1 (MQ=255)
                                                            cTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCa                                                                                         >  2:1295762/1‑150 (MQ=255)
                                                              gatgaCTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACAATACTAATAAGTCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCAcc                                                                                       <  2:719294/150‑1 (MQ=255)
                                                              gatgaCTATAAGTAATCTATAGAAATTTGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCAcc                                                                                       <  2:1718436/150‑1 (MQ=255)
                                                                 gaCTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCaaa                                                                                    >  2:146726/1‑150 (MQ=255)
                                                                  aCTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCTAac                                                                                   >  2:750178/1‑150 (MQ=255)
                                                                    tATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAccac                                                                                 >  2:519669/1‑150 (MQ=255)
                                                                       cAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACttt                                                                              <  2:712587/150‑1 (MQ=255)
                                                                        aGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTg                                                                             <  2:1433607/150‑1 (MQ=255)
                                                                         gTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGc                                                                            <  2:1143975/150‑1 (MQ=255)
                                                                                ataGAAAATTGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAAc                                                                     >  2:338963/1‑150 (MQ=255)
                                                                                 taGAAAATTGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACt                                                                    >  1:36645/1‑150 (MQ=255)
                                                                                   taaaaTTGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACttt                                                                  <  1:1252629/149‑1 (MQ=255)
                                                                                   gAAAATTGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACttt                                                                  >  2:1282167/1‑150 (MQ=255)
                                                                                    aaaaTTGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTg                                                                 <  1:12541/150‑1 (MQ=255)
                                                                                    aaaaTTGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTg                                                                 <  1:176230/150‑1 (MQ=255)
                                                                                    aaaaTTGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACAAAACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTg                                                                 <  2:319095/150‑1 (MQ=255)
                                                                                     aaaTTGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTgg                                                                <  1:1486324/150‑1 (MQ=255)
                                                                                     aaaTTGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTgg                                                                <  1:1526582/150‑1 (MQ=255)
                                                                                           cGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGt                                                          >  1:1370101/1‑150 (MQ=255)
                                                                                                cGAAACTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTcc                                                     <  2:1709704/150‑1 (MQ=255)
                                                                                                  aaaCTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAg                                                   <  1:312071/150‑1 (MQ=255)
                                                                                                    aCTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTc                                                 <  2:1151374/150‑1 (MQ=255)
                                                                                                       gCATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCggg                                              <  2:382514/150‑1 (MQ=255)
                                                                                                        cATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCTAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGt                                             <  1:750178/150‑1 (MQ=255)
                                                                                                        cATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGt                                             >  2:1523695/1‑150 (MQ=255)
                                                                                                        cATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGt                                             <  1:764358/150‑1 (MQ=255)
                                                                                                         aTTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGtt                                            <  1:1528489/150‑1 (MQ=255)
                                                                                                               gTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAgc                                      >  2:857155/1‑150 (MQ=255)
                                                                                                               gTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAgc                                      >  2:376867/1‑150 (MQ=255)
                                                                                                               gTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAgc                                      >  2:977453/1‑150 (MQ=255)
                                                                                                                tAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAgcg                                     >  1:1277957/1‑150 (MQ=255)
                                                                                                                 aa‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCTGGTTAACAgcgc                                    >  1:615023/1‑150 (MQ=255)
                                                                                                                 aa‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAgcgc                                    >  1:1098905/1‑150 (MQ=255)
                                                                                                                     tttttGAAACAATACAAATAAGGCATTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGTTCTCATTAACCAAACACTTTGCTTTCCACTTTGGTCAGGTAGTCAAGGCGGGTTAACAGCGCCTg                                 >  2:31864/1‑150 (MQ=255)
                                                                                                                         tGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAg                             <  1:1300416/150‑1 (MQ=255)
                                                                                                                         tGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAg                             >  2:1468476/1‑150 (MQ=255)
                                                                                                                               aaTACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTATCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTctgc                       >  2:1454251/1‑150 (MQ=255)
                                                                                                                               aaTACAAATAAGGCCTTCTCACTGAGAAGGCCATTAAAGACTTACAGTTAGATATCGCAAAATCCGAAAATACTATATACTTCACAAAACACTTTGCTTTCAAATATGGTAAGGAAGTCAAGTCGTGTTAACAGCGCCTGGTAGttatta                       >  2:330764/1‑146 (MQ=255)
                                                                                                                                aTACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCACACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTctgct                      >  1:1509152/1‑150 (MQ=255)
                                                                                                                                aTACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTctgct                      <  2:1025951/150‑1 (MQ=255)
                                                                                                                                    cggtAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCGTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGaa                  >  2:1152749/4‑150 (MQ=255)
                                                                                                                                     aaTAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAg                 >  1:909506/1‑150 (MQ=255)
                                                                                                                                     aaTAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAg                 >  1:445198/1‑150 (MQ=255)
                                                                                                                                         tggCCTTCTCAGTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAg             >  1:1199030/2‑150 (MQ=255)
                                                                                                                                         tggCCTTCTCAGTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAg             >  1:1390400/2‑150 (MQ=255)
                                                                                                                                            ccTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTAc          >  1:1015955/1‑150 (MQ=255)
                                                                                                                                            ccTTCTCACTGAGAACGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGAAGTCCAGTCGGGATAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTAc          >  1:357864/1‑150 (MQ=255)
                                                                                                                                               tctcACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTAt       >  1:1681097/1‑150 (MQ=255)
                                                                                                                                                  cACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTgg    >  1:583856/1‑150 (MQ=255)
                                                                                                                                                  cACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTgg    <  2:622839/150‑1 (MQ=255)
                                                                                                                                                  cACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTgg    <  2:421189/150‑1 (MQ=255)
                                                                                                                                                   aCTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAATCACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTgtt   >  1:805768/1‑148 (MQ=255)
                                                                                                                                                   aCTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGgt   >  1:424997/1‑150 (MQ=255)
                                                                                                                                                   aCTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGgt   <  2:1710873/150‑1 (MQ=255)
                                                                                                                                                   aCTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACATTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGgt   >  2:554264/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ATGTGCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAA‑TTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTG  >  NC_000913/292155‑292451

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