Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 4,957,549 2 bp→AT intergenic (‑371/‑300) PA14_RS22700 ← / → PA14_RS22705 16S ribosomal RNA/exonuclease SbcCD subunit D C‑terminal domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084634,957,5490GA99.7% 584.7 / NA 285intergenic (‑371/‑301)PA14_RS22700/PA14_RS2270516S ribosomal RNA/exonuclease SbcCD subunit D C‑terminal domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (118/165);  minor base T (1/0);  total (119/166)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0084634,957,5500CT100.0% 595.2 / NA 285intergenic (‑372/‑300)PA14_RS22700/PA14_RS2270516S ribosomal RNA/exonuclease SbcCD subunit D C‑terminal domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (119/166);  total (119/166)

Only 100 of 285 total aligned reads displayed.

TGCGCCGACCTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGATGAAGCTTCGTTTCGAGCCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGGAGATAGAAACT  >  NC_008463/4957455‑4957693
                                                                                              ||                                                                                                                                               
tGCGCCGACCTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACt                                                                                           <  1:8359645/150‑1 (MQ=2)
  cgcCGACCTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTcc                                                                                         <  1:3340194/150‑1 (MQ=12)
   gcCGACCTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCt                                                                                        <  1:5536288/150‑1 (MQ=12)
    ccGACCTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTg                                                                                       <  1:8639678/150‑1 (MQ=12)
         cTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGgaga                                                                                  >  2:9140114/1‑150 (MQ=12)
         cTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGgaga                                                                                  >  2:4858780/1‑150 (MQ=12)
           tCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGgagaga                                                                                <  2:7495570/150‑1 (MQ=12)
           tCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGgagaga                                                                                <  2:4783325/150‑1 (MQ=12)
             aaaCCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAAt                                                                              <  1:3098313/150‑1 (MQ=12)
              aaCCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATc                                                                             >  2:2806716/1‑150 (MQ=12)
               aCCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATcc                                                                            <  2:8980641/150‑1 (MQ=12)
               aCCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATcc                                                                            >  1:1589740/1‑150 (MQ=12)
               aCCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATcc                                                                            <  2:1084221/150‑1 (MQ=12)
                ccATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGAACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCg                                                                           >  1:8830921/1‑150 (MQ=255)
                 cATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGt                                                                          >  2:1901988/1‑150 (MQ=12)
                  aTCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTa                                                                         <  1:4559793/150‑1 (MQ=12)
                  aTCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTa                                                                         <  1:10142550/150‑1 (MQ=12)
                  aTCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTa                                                                         >  1:7485511/1‑150 (MQ=12)
                  aTCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTa                                                                         >  2:8256586/1‑150 (MQ=12)
                       ccGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATcc                                                                    <  1:7326149/150‑1 (MQ=12)
                        cGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCt                                                                   <  2:8811964/150‑1 (MQ=12)
                         gATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCtt                                                                  <  2:6536815/150‑1 (MQ=12)
                         gATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCtt                                                                  >  2:6691138/1‑150 (MQ=12)
                         gATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCtt                                                                  >  1:7823363/1‑150 (MQ=12)
                         gATCAGCGGGAGGCGCATGCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAACTCCtt                                                                  >  2:10632032/1‑150 (MQ=11)
                          aTCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTg                                                                 >  2:2936435/1‑150 (MQ=12)
                          aTCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTg                                                                 >  1:11007025/1‑150 (MQ=12)
                             aGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGAtt                                                              <  2:1589740/150‑1 (MQ=12)
                             aGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGAtt                                                              >  2:1520584/1‑150 (MQ=12)
                             aGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGAtt                                                              <  1:7079453/150‑1 (MQ=12)
                             aGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGAtt                                                              <  1:6879604/150‑1 (MQ=12)
                              gCGGGAGGCTCATTCTACCGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGAttt                                                             <  2:9912262/150‑1 (MQ=11)
                              gCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGAttt                                                             >  1:2488914/1‑150 (MQ=12)
                              gCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGAttt                                                             <  2:1702077/150‑1 (MQ=12)
                               cGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGAtttt                                                            <  1:10326824/150‑1 (MQ=12)
                               cGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGAtttt                                                            <  1:6193524/150‑1 (MQ=12)
                                gggAGGCGCATTCTACAGCGTTCAAGTTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTc                                                           >  2:10267375/1‑150 (MQ=255)
                                gggAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTc                                                           >  1:10493481/1‑150 (MQ=12)
                                gggAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTc                                                           >  1:1291579/1‑150 (MQ=12)
                                gggAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGATTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTc                                                           <  2:10375235/150‑1 (MQ=255)
                                 ggAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCCCCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCa                                                          >  2:4626401/1‑150 (MQ=255)
                                  gAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCaa                                                         <  2:4917904/150‑1 (MQ=12)
                                     gcgcATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCa                                                      <  1:7915974/150‑1 (MQ=12)
                                     gcgcATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCa                                                      >  1:3380456/1‑150 (MQ=12)
                                       gccttCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACTACCTGCCGGAGCGCGTCACCACTCCTGGGAGAGAATCCGTGACTCCTTGATTTTCAAGCagt                                                    >  2:4433383/4‑148 (MQ=11)
                                        cATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGACGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCAttt                                                   >  2:5250731/1‑150 (MQ=255)
                                         aTTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTa                                                  <  1:2121613/150‑1 (MQ=12)
                                          ttCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTAc                                                 <  1:3509229/150‑1 (MQ=12)
                                            cTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTActc                                               <  2:1701624/150‑1 (MQ=12)
                                            cTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTActc                                               <  1:1892802/150‑1 (MQ=12)
                                            cTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTActc                                               <  1:2975911/150‑1 (MQ=12)
                                              aCAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCtt                                             <  1:10189095/150‑1 (MQ=12)
                                               cAGCGTTCAAAGTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTtc                                            <  2:5778314/150‑1 (MQ=255)
                                                aGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTtcg                                           <  1:5075070/150‑1 (MQ=12)
                                                aGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGATTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTtcg                                           <  2:10412210/150‑1 (MQ=255)
                                                 gCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTtcgt                                          >  1:232268/1‑150 (MQ=12)
                                                   gTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTtcgtcg                                        <  1:7456877/150‑1 (MQ=12)
                                                    ttCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTtcgtcgt                                       >  1:650619/1‑150 (MQ=12)
                                                    ttCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTtcgtcgt                                       <  1:2760579/150‑1 (MQ=12)
                                                    ttCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTtcgtcgt                                       <  1:9255980/150‑1 (MQ=12)
                                                        aaTTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGcc                                   <  1:9530223/150‑1 (MQ=12)
                                                        aaTTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGcc                                   <  1:9528567/150‑1 (MQ=12)
                                                        aaTTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCCGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGcc                                   >  1:10118357/1‑150 (MQ=255)
                                                         aTTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCg                                  <  1:1762043/150‑1 (MQ=12)
                                                         aTTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCg                                  >  2:9337973/1‑150 (MQ=12)
                                                         aTTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCg                                  <  1:3773367/150‑1 (MQ=12)
                                                         aTTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCg                                  <  1:1128320/150‑1 (MQ=12)
                                                         aTTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCg                                  >  2:5896276/1‑150 (MQ=12)
                                                          ttCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCgg                                 <  2:3670832/150‑1 (MQ=12)
                                                             gCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAg                              <  2:8667558/150‑1 (MQ=12)
                                                               tGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAgtg                            <  1:638554/150‑1 (MQ=12)
                                                                gTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAgtgg                           <  2:10067132/150‑1 (MQ=12)
                                                                gTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAgtgg                           <  1:9892527/150‑1 (MQ=12)
                                                                 tCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAgtggt                          <  2:4644734/150‑1 (MQ=12)
                                                                 tCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAgtggt                          >  2:3482327/1‑150 (MQ=12)
                                                                 tCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCAGTTACTCTTCGTCGTCGCCGGGAgtggt                          >  2:10203408/1‑150 (MQ=255)
                                                                       aCCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCAt                    <  1:7189211/150‑1 (MQ=12)
                                                                        ccTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCAtt                   >  1:9964335/1‑150 (MQ=12)
                                                                        ccTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCAtg                   >  2:8439039/1‑149 (MQ=12)
                                                                         cTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATta                  <  2:2051315/150‑1 (MQ=12)
                                                                          tCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATtat                 <  2:3553940/150‑1 (MQ=12)
                                                                           cggtGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATtata                >  2:78642/4‑150 (MQ=12)
                                                                            ggtGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAg               >  2:8008747/3‑150 (MQ=12)
                                                                             gtGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAgg              >  2:6130445/2‑150 (MQ=17)
                                                                              tGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGACTCCGTAAATCCTTGATGTTCAAGCAGTTACTCTTCGTCGGCGCCGGGAGTGGTGCGCAGTTTAggg             >  2:7012921/1‑150 (MQ=11)
                                                                              tGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAggg             <  2:5437049/150‑1 (MQ=17)
                                                                               gAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGga            <  2:1441258/150‑1 (MQ=17)
                                                                                 aGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGgaga          <  1:4974334/150‑1 (MQ=17)
                                                                                 aGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGgaga          >  2:3233799/1‑150 (MQ=17)
                                                                                 aGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGgaga          <  2:4204833/150‑1 (MQ=17)
                                                                                  gCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGGAGAt         <  1:6596947/150‑1 (MQ=17)
                                                                                  gCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGGAGAt         >  1:6001744/1‑150 (MQ=17)
                                                                                   cTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGGAGATa        <  2:3072540/150‑1 (MQ=17)
                                                                                   cTTCGTTTCGAATCCCTTTCGGGACCGAACCGCGGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGCTTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTCTAGGGAGATa        <  2:5159826/150‑1 (MQ=11)
                                                                                    ttCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGGAGATAg       <  1:5615772/150‑1 (MQ=17)
                                                                                     tCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGGAGATAGa      >  1:690262/1‑150 (MQ=17)
                                                                                     tCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGGAGATAGa      <  1:4298547/150‑1 (MQ=17)
                                                                                     tCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGGAGATAGa      >  2:2410002/1‑150 (MQ=17)
                                                                                         ttCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGGAGATAGAAACt  >  1:3980273/1‑150 (MQ=17)
                                                                                         ttCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGGAGATAGAAACt  >  2:2920026/1‑150 (MQ=17)
                                                                                              ||                                                                                                                                               
TGCGCCGACCTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGATGAAGCTTCGTTTCGAGCCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGGAGATAGAAACT  >  NC_008463/4957455‑4957693

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