Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 3,390,485 Δ1 bp pseudogene (451/1997 nt) PA14_RS15455 → alkaline phosphatase D family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084633,390,4830G.99.3% 2202.9 / 3.2 453pseudogene (449/1997 nt)PA14_RS15455alkaline phosphatase D family protein
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base . (215/233);  total (216/237)
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 469 total aligned reads displayed.

CGCCGACGGCGGGGAGCCCGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGGCGGCCGATGGGCCTGCTGCAGGTAGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGCCATCCATGGGCTGATCCGCC  >  NC_008463/3390334‑3390628
                                                                                                                                                     |                                                                                                                                                 
cgccgACGGCGGGGAGCCCGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATgg                                                                                                                                                  <  1:10990860/150‑1 (MQ=255)
 gccgACGGCGGGGAGCCCGCGGGAATCGCCGACTGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCCGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATggc                                                                                                                                                 >  2:3318252/1‑149 (MQ=255)
 gccgACGGCGGGGAGCCCGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATggc                                                                                                                                                 >  1:9066473/1‑149 (MQ=255)
     aCGGCGGGGAGCCCGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCctgc                                                                                                                                            >  2:2315166/1‑150 (MQ=255)
        gcggGGAGCCCGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCtgctgc                                                                                                                                         >  1:4614973/1‑150 (MQ=255)
         cggGGAGCCCGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCa                                                                                                                                        >  1:2055002/1‑150 (MQ=255)
           gggAGCCCGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCCCGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAgg                                                                                                                                      >  1:4372029/1‑150 (MQ=255)
            ggAGCCCGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGt                                                                                                                                     <  1:6981200/150‑1 (MQ=255)
                 ccGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGAcc                                                                                                                                >  1:9477098/1‑150 (MQ=255)
                 ccGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGAcc                                                                                                                                >  1:432513/1‑150 (MQ=255)
                  cgcgGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCg                                                                                                                               >  1:10672348/1‑150 (MQ=255)
                    cgGGAATCGCCGAATGGGCGCCGCATCTGCTCGACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGgc                                                                                                                             <  2:11038878/150‑1 (MQ=255)
                     gggAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGgct                                                                                                                            <  1:3812542/150‑1 (MQ=255)
                        aaTCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGgctgct                                                                                                                         <  1:11103784/150‑1 (MQ=255)
                        aaTCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGgctgct                                                                                                                         <  1:7531805/150‑1 (MQ=255)
                          tCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTgg                                                                                                                       >  1:8533762/1‑150 (MQ=255)
                           cGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTgga                                                                                                                      <  1:5192877/150‑1 (MQ=255)
                           cGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCAGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTgga                                                                                                                      <  1:6840622/150‑1 (MQ=255)
                             ccGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTggagg                                                                                                                    >  2:8668348/1‑150 (MQ=255)
                               gAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGaa                                                                                                                  >  2:3611356/1‑150 (MQ=255)
                                aaTGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGaaa                                                                                                                 >  2:8786087/1‑150 (MQ=255)
                                aaTGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGaaa                                                                                                                 <  2:2041483/150‑1 (MQ=255)
                                  tGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAAc                                                                                                               <  1:7749931/150‑1 (MQ=255)
                                  tGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAAc                                                                                                               <  2:9227308/150‑1 (MQ=255)
                                   gggCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAAcc                                                                                                              <  2:2742216/150‑1 (MQ=255)
                                   gggCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAAcc                                                                                                              >  1:10294038/1‑150 (MQ=255)
                                      cgccgcATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACctgc                                                                                                           <  1:2606226/150‑1 (MQ=255)
                                       gccgcATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCAACACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACctgct                                                                                                          >  2:6378072/1‑150 (MQ=255)
                                          gcATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCa                                                                                                       <  1:4705338/150‑1 (MQ=255)
                                          gcATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCa                                                                                                       <  2:7904455/150‑1 (MQ=255)
                                           cATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAg                                                                                                      >  1:3067494/1‑150 (MQ=255)
                                             tCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGc                                                                                                    <  2:9679528/150‑1 (MQ=255)
                                             tCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGc                                                                                                    >  2:4230326/1‑150 (MQ=255)
                                             tCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGc                                                                                                    <  1:712939/150‑1 (MQ=255)
                                             tCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGc                                                                                                    >  1:4113346/1‑150 (MQ=255)
                                             tCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGc                                                                                                    >  1:1453193/1‑150 (MQ=255)
                                                    tACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAgc                                                                                             <  1:5577163/150‑1 (MQ=255)
                                                     aCCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAgcc                                                                                            >  1:974921/1‑150 (MQ=255)
                                                      ccccGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAgccg                                                                                           >  2:6650445/1‑150 (MQ=255)
                                                       cccGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAgccgc                                                                                          <  1:6298639/150‑1 (MQ=255)
                                                       cccGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAgccgc                                                                                          >  2:10583740/1‑150 (MQ=255)
                                                       cccGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAgccgc                                                                                          <  2:8484277/150‑1 (MQ=255)
                                                        ccGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAgccgcc                                                                                         >  1:8293339/1‑150 (MQ=255)
                                                          ggTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCg                                                                                       >  2:544625/1‑150 (MQ=255)
                                                           gTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGGTGCTCCACGGCTCCTGCCGCAAGCCGCAGCACCCGGCGGCCGATGGCCTGCTGCAGGTCGCCCGGCTGCGGGAGGACATCCTGCTGCAGGCCGAAAGCCGCCCgg                                                                                      >  2:4308938/1‑150 (MQ=255)
                                                            tCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGc                                                                                     >  2:1186377/1‑150 (MQ=255)
                                                              gCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGt                                                                                   <  1:6451950/150‑1 (MQ=255)
                                                              gCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGt                                                                                   <  2:5638041/150‑1 (MQ=255)
                                                               ccAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGtt                                                                                  >  1:899916/1‑150 (MQ=255)
                                                                 aGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGc                                                                                >  1:3637787/1‑150 (MQ=255)
                                                                 aGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGc                                                                                >  2:4847510/1‑150 (MQ=255)
                                                                 aGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGCCCTGCTGCAGGTCGACCGGCTGCTGGAGGCAAACCTGCTGCAGGCCGAAAGCCGCCCGGAGTTGc                                                                                >  2:9022472/1‑150 (MQ=255)
                                                                    gcgcCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCtga                                                                             >  1:5371878/1‑150 (MQ=255)
                                                                    gcgcCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCtga                                                                             <  2:3905796/150‑1 (MQ=255)
                                                                    gcgcCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCtga                                                                             <  2:823699/150‑1 (MQ=255)
                                                                      gcCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCtgatg                                                                           <  2:4128269/150‑1 (MQ=255)
                                                                      gcCAGCCTGGTGCACACGCCACGTCTCGACCAGTTGCTCCACGGCTCCAGCCGCAAGCCGCATCACCCGGCGGCCGAGGGCCTGCTGCAGGTCGCCCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGCAAGCCGCCCGGCGTTGCtgatg                                                                           <  2:3138932/150‑1 (MQ=255)
                                                                       ccAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCtgatga                                                                          >  2:301738/1‑150 (MQ=255)
                                                                          gCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCtgatgatga                                                                       <  1:5957634/150‑1 (MQ=255)
                                                                          gCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCtgatgatga                                                                       <  1:739860/150‑1 (MQ=255)
                                                                           ccaggTGCTCAAGCCACGTCTCGACCAGTGGCTCCCCGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAACACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAg                                                                      <  2:7794247/147‑1 (MQ=255)
                                                                             tGGTGCTCACGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCCAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAgcg                                                                    <  2:9452536/150‑1 (MQ=255)
                                                                                 gCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGa                                                                <  1:8773176/150‑1 (MQ=255)
                                                                                     caGCCACGTCTCGACCAGTTGCTCCCCGGCTCCTGCCGCCAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCGGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAg                                                            <  1:4764212/149‑1 (MQ=255)
                                                                                     aaGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAATCCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAg                                                            >  2:6934377/1‑150 (MQ=255)
                                                                                        ccACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTc                                                         <  2:4963185/150‑1 (MQ=255)
                                                                                             tctcGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGc                                                    <  1:2258607/150‑1 (MQ=255)
                                                                                             tctcGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCGGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTGGCGGATGATGAGCGGCGACCCGGTCTACGc                                                    >  2:10237877/1‑150 (MQ=255)
                                                                                              ctcGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGcc                                                   <  1:10041305/150‑1 (MQ=255)
                                                                                                cGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCcga                                                 <  1:8681671/150‑1 (MQ=255)
                                                                                                cGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGGTGCTGATGATGAGCGGCGACCAGGTCTACGCcga                                                 >  1:5494716/1‑150 (MQ=255)
                                                                                                 gACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCcgac                                                <  2:2643509/150‑1 (MQ=255)
                                                                                                  aCCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCcgacg                                               >  2:10875669/1‑150 (MQ=255)
                                                                                                   ccAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCcgacga                                              <  1:996526/150‑1 (MQ=255)
                                                                                                    cAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGCCCAGGTCTACGCCgacgac                                             >  2:5992343/1‑150 (MQ=255)
                                                                                                      gTTGCTCCACGGCTCCTGCCGCACGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGt                                           <  1:3862451/150‑1 (MQ=255)
                                                                                                         gCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGc                                        <  2:2816507/150‑1 (MQ=255)
                                                                                                         gCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGc                                        <  2:1761509/150‑1 (MQ=255)
                                                                                                         gCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGc                                        <  1:1228520/150‑1 (MQ=255)
                                                                                                          cTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGcc                                       <  1:4983654/150‑1 (MQ=255)
                                                                                                            ccACGGCTCCTGCCGCCAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCgg                                     <  1:6743451/150‑1 (MQ=255)
                                                                                                             cACGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCggg                                    <  1:722834/150‑1 (MQ=255)
                                                                                                              aCGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGc                                   >  1:10763419/1‑150 (MQ=255)
                                                                                                              aCGGCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGGCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGGTGCTGCTGATGAGCGGCGTCGAGGTCTACGCCGTCGACGTCGCCGGGc                                   >  1:4817299/1‑150 (MQ=255)
                                                                                                                 gCTCCTGCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGCCGTCGCCGGGCCGa                                >  2:7573226/1‑150 (MQ=255)
                                                                                                                       gCCGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCtgg                          <  2:7201054/150‑1 (MQ=255)
                                                                                                                         cGCAAGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCtggtg                        <  2:3517077/150‑1 (MQ=255)
                                                                                                                             aGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGcca                    <  2:6177914/150‑1 (MQ=255)
                                                                                                                             aGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGcca                    <  2:1516081/150‑1 (MQ=255)
                                                                                                                             aGCCGCATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGcca                    >  1:10495804/1‑150 (MQ=255)
                                                                                                                                  cATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGccatccat               >  1:194137/1‑150 (MQ=255)
                                                                                                                                  cATCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCGGGTGGccatccat               >  1:2719798/1‑150 (MQ=255)
                                                                                                                                   aTCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGCCATCCATg              >  2:1325351/1‑150 (MQ=255)
                                                                                                                                    tCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGCCATCCATgg             >  1:6139049/1‑150 (MQ=255)
                                                                                                                                    tCACCCGGCGGCCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGAAGACGTCGCCGGGCCGATGCTGGTGGCCATCCATgg             >  1:5739311/1‑150 (MQ=255)
                                                                                                                                        cccggcggcCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGCCATCCCTGGGCt          >  2:7123130/4‑150 (MQ=255)
                                                                                                                                               gcCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGCCATCCATGGGCTGATccgc   <  1:9343327/150‑1 (MQ=255)
                                                                                                                                               gcCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGCCATCCATGGGCTGATccgc   >  2:6584125/1‑150 (MQ=255)
                                                                                                                                               gcCGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGCCATCCATGGGCTGATccgc   >  2:4571591/1‑150 (MQ=255)
                                                                                                                                                ccGATGGCCTGCTGCAGGTCGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGCCATCCATGGGCTGATccgcc  <  1:8098889/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
CGCCGACGGCGGGGAGCCCGCGGGAATCGCCGAATGGGCGCCGCATCTGCTCTACCCCGGTCGCCAGCGCGCCAGCCTGGTGCTCAAGCCACGTCTCGACCAGTTGCTCCACGGCTCCTGCCGCAAGCCGCATCACCCGGGCGGCCGATGGGCCTGCTGCAGGTAGACCGGCTGCTGGAGGAAAACCTGCTGCAGGCCGAAAGCCGCCCGGCGTTGCTGATGATGAGCGGCGACCAGGTCTACGCCGACGACGTCGCCGGGCCGATGCTGGTGGCCATCCATGGGCTGATCCGCC  >  NC_008463/3390334‑3390628

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: