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% 454.9 / NA 282intergenic (‑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 (126/155);  minor base T (0/1);  total (126/156)
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
*NC_0084634,957,5500CT99.7% 460.0 / NA 282intergenic (‑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);  major base T (126/155);  minor base G (0/1);  total (126/156)
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

Only 100 of 282 total aligned reads displayed.

TTGCGCCGACCTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGATGAAGCTTCGTTTCGAGCCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGGAGATAGAAACTG  >  NC_008463/4957454‑4957694
                                                                                               ||                                                                                                                                                
ttGCGCCGACCTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACt                                                                                            <  2:9160770/151‑1 (MQ=2)
 tGCGCCGACCTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTc                                                                                           <  1:7940740/151‑1 (MQ=9)
   cgcCGACCTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCt                                                                                         <  1:9073146/151‑1 (MQ=9)
    gcCGACCTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTg                                                                                        >  2:9274517/1‑151 (MQ=9)
     ccGACCTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACTACCGGCAGGAGCGCTTCACAACTCCTgg                                                                                       >  2:8912942/1‑151 (MQ=11)
       gACCTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGACCAACCTGCCGGAGCGCTTCACCACTCCTGGga                                                                                     >  2:4068236/1‑151 (MQ=9)
        aCCTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGgag                                                                                    <  1:6140418/151‑1 (MQ=9)
        aCCTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGgag                                                                                    >  2:2360313/1‑151 (MQ=9)
        aCCTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGgag                                                                                    <  1:1189138/151‑1 (MQ=9)
        aCCTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGgag                                                                                    <  1:1186701/151‑1 (MQ=9)
          cTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGgagag                                                                                  >  2:2447473/1‑151 (MQ=9)
          cTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGgagag                                                                                  <  1:4571251/151‑1 (MQ=9)
           ttCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGgagaga                                                                                 <  2:10623058/151‑1 (MQ=9)
           ttCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGgagaga                                                                                 >  1:9628335/1‑151 (MQ=9)
            tCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGaa                                                                                >  2:1102066/1‑151 (MQ=9)
               aaCCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATcc                                                                             <  1:1137723/151‑1 (MQ=9)
               aaCCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATcc                                                                             <  1:9852280/151‑1 (MQ=9)
                aCCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCg                                                                            <  1:5368079/151‑1 (MQ=9)
                 ccATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGACCAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTGTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGt                                                                           <  2:9518409/151‑1 (MQ=11)
                  cATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTa                                                                          >  1:8380616/1‑151 (MQ=9)
                  cATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTa                                                                          <  1:10645755/151‑1 (MQ=9)
                   aTCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTaa                                                                         >  2:5131940/1‑151 (MQ=9)
                     cTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAAt                                                                       <  2:1251722/151‑1 (MQ=9)
                     cTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAAt                                                                       >  1:5579120/1‑151 (MQ=9)
                      ttCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATc                                                                      <  1:1378255/151‑1 (MQ=9)
                      ttCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATc                                                                      <  2:2503421/151‑1 (MQ=9)
                      ttCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATc                                                                      >  2:6798201/1‑151 (MQ=9)
                      ttCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATc                                                                      >  2:5665433/1‑151 (MQ=9)
                       tCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATcc                                                                     >  2:10339512/1‑151 (MQ=9)
                        ccGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCt                                                                    >  2:5735895/1‑151 (MQ=9)
                        ccGATCAGCGGGAGGCGCAATCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCt                                                                    >  2:9352850/1‑151 (MQ=9)
                         cGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCtt                                                                   >  2:6057651/1‑151 (MQ=9)
                          gATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTg                                                                  <  2:7202112/151‑1 (MQ=9)
                           aTCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGa                                                                 <  1:2670457/151‑1 (MQ=9)
                           aTCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGa                                                                 >  1:5153250/1‑151 (MQ=9)
                           aTCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGa                                                                 >  1:394932/1‑151 (MQ=9)
                            tCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGAt                                                                <  1:10774923/151‑1 (MQ=9)
                             cAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGAtt                                                               <  2:7838229/151‑1 (MQ=9)
                             cAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGAtt                                                               >  2:1534977/1‑151 (MQ=9)
                             cAGCGGGAGGCGCAGTCTACAGCGTTCAAAGTCGCTGTCAAGCACCTCGGTGAAGCTTCGGTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGAtt                                                               <  2:6582474/151‑1 (MQ=11)
                              aGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGAttt                                                              <  1:8047156/151‑1 (MQ=9)
                               gCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGAtttt                                                             <  1:7861973/151‑1 (MQ=9)
                               gCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGAtttt                                                             <  2:1750917/151‑1 (MQ=9)
                               gCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGAtttt                                                             <  2:8183634/151‑1 (MQ=9)
                                  ggAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCaa                                                          <  1:8775143/151‑1 (MQ=9)
                                   gAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAg                                                         >  2:2103285/1‑151 (MQ=9)
                                      gcgcATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCAt                                                      <  2:379871/151‑1 (MQ=9)
                                       cgcATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCAtt                                                     >  2:7772583/1‑151 (MQ=9)
                                        gcATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCAttt                                                    >  1:998928/1‑151 (MQ=9)
                                          aTTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTAc                                                  >  2:6941677/1‑151 (MQ=9)
                                          aTTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTAc                                                  >  1:2203818/1‑151 (MQ=9)
                                            tCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTActc                                                >  2:5503585/1‑151 (MQ=9)
                                            tCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCACTTTCTGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTActc                                                >  2:3428826/1‑151 (MQ=255)
                                               aCAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTtc                                             <  1:9179425/151‑1 (MQ=9)
                                               aCAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTtc                                             <  1:5385631/151‑1 (MQ=9)
                                                 aGCGTTCAAAGTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTtcgt                                           <  2:7515827/151‑1 (MQ=9)
                                                  gCGTTCAAATTCGCTGTCAAGCTCCTCGGTGAATCTTCGTTTCGAATCCCTTTAGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGTGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTtcgtc                                          <  1:10003147/151‑1 (MQ=11)
                                                  gCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTtcgtc                                          <  2:1334411/151‑1 (MQ=9)
                                                  gCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTtcgtc                                          <  1:1733038/151‑1 (MQ=9)
                                                   cGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTtcgtcg                                         <  2:4521506/151‑1 (MQ=9)
                                                    gTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTtcgtcgt                                        <  2:4374959/151‑1 (MQ=9)
                                                    gTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTtcgtcgt                                        <  1:9379866/151‑1 (MQ=9)
                                                    gTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTtcgtcgt                                        <  1:864203/151‑1 (MQ=9)
                                                     ttCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTtcgtcgtc                                       <  1:2344162/151‑1 (MQ=9)
                                                     ttCAAATTCGCTGTCAAGCACCTCGGTGAAGCTACGTTTCGAATCACTTTCGGGACCGACCCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTtcgtcgtc                                       >  2:1945718/1‑151 (MQ=11)
                                                       cAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGc                                     >  2:3643067/1‑151 (MQ=9)
                                                          aTTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCgg                                  >  2:6210788/1‑151 (MQ=9)
                                                           ttCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCggg                                 <  1:1922471/151‑1 (MQ=9)
                                                              gCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAgt                              <  1:9877445/151‑1 (MQ=9)
                                                              gCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAgt                              >  1:2060423/1‑151 (MQ=9)
                                                               cTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAgtg                             <  2:3206453/151‑1 (MQ=9)
                                                               cTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAgtg                             <  1:4442185/151‑1 (MQ=9)
                                                                  tCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAgtggtg                          <  1:8941936/151‑1 (MQ=9)
                                                                   cAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTgc                         >  2:10712222/1‑151 (MQ=9)
                                                                    aaGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTgcg                        >  2:4912869/1‑151 (MQ=9)
                                                                       cACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCAt                     >  1:1131264/1‑151 (MQ=9)
                                                                            cggtGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAg                <  2:7318461/148‑1 (MQ=9)
                                                                             ggtGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAgg               <  1:328303/149‑1 (MQ=9)
                                                                             ggtGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAgg               <  1:4681304/149‑1 (MQ=9)
                                                                             ggtGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCAGTTACTCGTCGTCGTCGCCGGGAGTGGTGCGCATTATCgg               >  1:6012672/3‑151 (MQ=11)
                                                                               tGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGga             <  1:10372101/151‑1 (MQ=14)
                                                                               tGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGga             <  1:3853444/151‑1 (MQ=14)
                                                                               tGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGga             >  1:7468051/1‑151 (MQ=14)
                                                                               tGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCCGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGga             <  1:7273263/151‑1 (MQ=9)
                                                                                gAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGgag            >  1:6299815/1‑151 (MQ=14)
                                                                                gAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGgag            >  2:10794054/1‑151 (MQ=14)
                                                                                gAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGgag            <  2:8785659/151‑1 (MQ=14)
                                                                                gAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGgag            >  2:8073718/1‑151 (MQ=14)
                                                                                 aaGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGGCGGGAGTGGTGCGCAGTATCGGgaga           >  2:2747296/1‑151 (MQ=255)
                                                                                 aaGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGgaga           <  1:3073010/151‑1 (MQ=14)
                                                                                  aGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGGAGAt          <  1:4157446/151‑1 (MQ=14)
                                                                                    cTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGGAGATAg        >  1:7663765/1‑151 (MQ=14)
                                                                                    cTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGGAGATAg        <  1:2931094/151‑1 (MQ=14)
                                                                                     ttCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGGAGATAGa       <  2:8526806/151‑1 (MQ=14)
                                                                                       cGTTTCGAATCCCTTTCGGGCCCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCGACACCACTCCTGGGAGAGAATCCGTAAATCCTTGCTTTTCAAGCATTTACACTTCGTCGTCGCCGGGAGTGGTGCGCATTCTAGGGAGATAGaaa     <  2:8706662/151‑1 (MQ=11)
                                                                                       cGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGGAGATAGaaa     <  1:5795276/151‑1 (MQ=14)
                                                                                        gTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGGAGATAGAAAc    >  1:6666140/1‑151 (MQ=14)
                                                                                        gTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGGAGATAGAAAc    >  1:1663552/1‑151 (MQ=14)
                                                                                         tttCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGGAGATAGAAACt   <  1:3694434/151‑1 (MQ=14)
                                                                                          ttCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGGAGAGAGAAACTg  >  2:3829571/1‑151 (MQ=9)
                                                                                               ||                                                                                                                                                
TTGCGCCGACCTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGATGAAGCTTCGTTTCGAGCCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGGAGATAGAAACTG  >  NC_008463/4957454‑4957694

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: