Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 1,128,162 N→C intergenic (+23/‑38) PA14_RS05280 → / → PA14_RS05285 DUF1302 domain‑containing protein/DUF1329 domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084631,128,1620NC99.6% 1721.6 / NA 463intergenic (+23/‑38)PA14_RS05280/PA14_RS05285DUF1302 domain‑containing protein/DUF1329 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base N (0/0);  major base C (218/242);  minor base G (0/2);  total (218/245)
Fisher's exact test for biased strand distribution p-value = 5.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 470 total aligned reads displayed.

GGGGGTGGACGCGGAGTACCAGAACACCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTNNNANCNGGNCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAGGCGGCCAAGCTCGGCACCAGCC  >  NC_008463/1128013‑1128311
                                                                                                                                                     |                                                                                                                                                     
gggggTGGACGCGGAGTACCAGAACACCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCg                                                                                                                                                      <  1:1340206/151‑1 (MQ=255)
gggggTGGACGCGGAGTACCAGAACACCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCg                                                                                                                                                      >  2:4748644/1‑151 (MQ=255)
 ggggTGGACGCGGAGTACCAGAACACCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCgg                                                                                                                                                     <  1:7207532/151‑1 (MQ=255)
 ggggTGGACGCGGAGTACCAGAACACCTACACCGCGAGCCTCACCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCGGATGTACGCCGGCCCGGTGCGAACCgg                                                                                                                                                     <  2:3482600/151‑1 (MQ=255)
  gggTGGACGCGGAGTACCAGAACACCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCggg                                                                                                                                                    >  2:9906511/1‑150 (MQ=255)
  gggTGGACGCGGAGTACCAGAACACCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCggg                                                                                                                                                    >  1:11154658/1‑150 (MQ=255)
    gTGGACGCGGAGTACCAGAACACCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGcc                                                                                                                                                  >  1:7547729/1‑151 (MQ=255)
     tGGACGCGGAGTACCAGAACACCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCg                                                                                                                                                 >  2:7838456/1‑151 (MQ=255)
       gACGCGGAGTACCAGAACACCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGc                                                                                                                                               <  1:9417649/151‑1 (MQ=255)
       gACGCGGAGTACCAGAACACCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGc                                                                                                                                               >  2:2204028/1‑151 (MQ=255)
        aCGCGGAGTACCAGAACACCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTAAGCCGGCCCGGTGCGAACCGGGCCGGCg                                                                                                                                              >  2:9544934/1‑151 (MQ=255)
               gtcccAGAACACCTACACCGCGAGCCTCTCCTACACCACCTTCTTCGACGGCAAGTACCCCCCCGTGGCCGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGa                                                                                                                                       <  2:10476964/148‑1 (MQ=255)
               gTACCAGAACACCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGa                                                                                                                                       <  2:5250804/151‑1 (MQ=255)
                 aCCAGAACACCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACCCCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACt                                                                                                                                     >  1:1195041/1‑151 (MQ=255)
                  ccAGAACACCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACtt                                                                                                                                    <  2:2828436/151‑1 (MQ=255)
                   cAGAACACCTACACCGCGCGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCCCCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTa                                                                                                                                   <  2:11277432/151‑1 (MQ=255)
                         acCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAgg                                                                                                                             <  1:6451270/151‑1 (MQ=255)
                          ccTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACATCGTCGCGCTCAGCTTCGGCATGAACTTATGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGga                                                                                                                            >  1:7869395/1‑151 (MQ=255)
                              cacCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGgacgac                                                                                                                        <  1:2404681/151‑1 (MQ=255)
                              cacCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACATCGTCGCGCTCAGCTTCGGCATTAACTTATGATGTACGCCGGCCCCGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGgacgac                                                                                                                        >  1:10810361/1‑151 (MQ=255)
                               acCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGacgacg                                                                                                                       >  2:10795146/1‑151 (MQ=255)
                                 cgcgAGCCTCTCCTACACCCACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAg                                                                                                                     <  2:8585307/151‑1 (MQ=255)
                                 cgcgAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATAGCGAATTCGTAGCGCTCAGCTTAGGCATGAAATTCTGATGTAAGCCGGCCCGGTTCGAACCGGGCCGGCGTGCAGGACTTAACTATGACGAagag                                                                                                                     >  1:5417089/1‑151 (MQ=255)
                                  gcgAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGc                                                                                                                    <  1:6754726/151‑1 (MQ=255)
                                  gcgAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTGCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGc                                                                                                                    <  2:3060343/151‑1 (MQ=255)
                                   cgAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCg                                                                                                                   <  1:6152463/151‑1 (MQ=255)
                                    gAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGa                                                                                                                  <  2:9503108/151‑1 (MQ=255)
                                    gAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGa                                                                                                                  >  1:743313/1‑151 (MQ=255)
                                    gAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTAGTCGCGCTCAGCTTCGGCATTAACTTCTGATGTACGCCGGCCCGGGGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGa                                                                                                                  >  1:3273702/1‑151 (MQ=255)
                                    gAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACGTCGTCGCGCACAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGTGGACGTCGAGCGa                                                                                                                  >  1:449100/1‑151 (MQ=255)
                                     aGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGaa                                                                                                                 <  1:10084460/151‑1 (MQ=255)
                                     aGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGaa                                                                                                                 >  2:11451956/1‑151 (MQ=255)
                                      gCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAAt                                                                                                                <  2:418399/151‑1 (MQ=255)
                                      gCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGGGCGAACCGGGCCGGCGAGCAGGACTTAAAGAGGACGACGAGCGAAt                                                                                                                >  1:572985/1‑151 (MQ=255)
                                       ccTCTCCTCCACCATCTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGAAGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATg                                                                                                               <  2:9313960/151‑1 (MQ=255)
                                       ccTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATg                                                                                                               >  2:6622409/1‑151 (MQ=255)
                                            ccTACCCCACCTACTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAaac                                                                                                          <  2:7453354/151‑1 (MQ=255)
                                            ccTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAaac                                                                                                          <  2:4294067/151‑1 (MQ=255)
                                            ccTACACCAACTTATTCGACGGCAAGTACACCACCGTGGACTATCGCGACTTCGTCGCGCTCAGCTTCGGCATGACATTCTGATGTACGCCGGCCCGTTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAaac                                                                                                          <  1:4771913/151‑1 (MQ=255)
                                               acacCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAaacaac                                                                                                       >  2:8482556/1‑151 (MQ=255)
                                               acacCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAaacaac                                                                                                       >  2:11110057/1‑151 (MQ=255)
                                               acacCAACTTCTTCGACGCCAAGTACACCCCCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTACTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAaacaac                                                                                                       <  1:9998473/151‑1 (MQ=255)
                                                cacCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACg                                                                                                      <  2:1691664/151‑1 (MQ=255)
                                                 acCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGa                                                                                                     <  2:9101977/151‑1 (MQ=255)
                                                  ccAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGaa                                                                                                    >  1:8383105/1‑151 (MQ=255)
                                                   cAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGAATTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGTGCGAATGAAAACAACGAAg                                                                                                   >  2:11368860/1‑151 (MQ=255)
                                                    aaCTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTTCGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGa                                                                                                  >  1:6627474/1‑151 (MQ=255)
                                                      cttcttCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGAtt                                                                                                <  1:10756678/151‑1 (MQ=255)
                                                          ttcGACGGCAAGTACACCACCGTGGACGATCGCGTCTTCGTCGCGCGCAGCGTAGGCATGAACTTGTGATGTACGCCGGCCCGGGGGGCACCGGGCCGGCGAGCCGGACTTTACGAGGACGACGAGCGAATGTCAACCACGAAGAttcgtc                                                                                            >  2:3803951/1‑151 (MQ=255)
                                                              aCGGCACGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCAGGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTATCGAGGACGACGAGCGAATGAACACAACGAAGAGTCTTCTGCa                                                                                        >  2:2696239/1‑151 (MQ=255)
                                                              aCGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCa                                                                                        <  1:9875594/151‑1 (MQ=255)
                                                              aCGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCa                                                                                        <  2:2472672/151‑1 (MQ=255)
                                                                ggCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATa                                                                                      <  1:5537370/151‑1 (MQ=255)
                                                                  cAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACGTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATAcc                                                                                    >  2:8219377/1‑151 (MQ=255)
                                                                    aGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTTATGTACGCAGGCCCGGTGCGAACCGGGCCGGCGAGCATGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCgg                                                                                  <  2:2311500/151‑1 (MQ=255)
                                                                    aGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCgg                                                                                  >  2:6767312/1‑151 (MQ=255)
                                                                         accaccGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTAc                                                                             >  2:3661768/1‑151 (MQ=255)
                                                                            accGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTgg                                                                          <  2:738042/151‑1 (MQ=255)
                                                                             ccGTGGACGAGCGCGACTTCGTCGCGCTCAGCATCGGCCAGAACTTAATATGTACGCCCGCCCGGTGCGAACCGGGCCGGCGAGCAGGACATAAAGAGGACGACGAGCGAATGAAAACAACGAATATTCTTCTTCATACCGGTGTACAGGc                                                                         >  1:11022643/1‑151 (MQ=255)
                                                                              cGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTCACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGcc                                                                        <  2:6059172/151‑1 (MQ=255)
                                                                               gTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGccc                                                                       <  1:852657/151‑1 (MQ=255)
                                                                                  gACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGa                                                                    <  2:4602400/151‑1 (MQ=255)
                                                                                    cGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGc                                                                  >  1:4886719/1‑151 (MQ=255)
                                                                                     gATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGcc                                                                 >  1:10525832/1‑151 (MQ=255)
                                                                                     gATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGcc                                                                 >  1:6929949/1‑151 (MQ=255)
                                                                                      aTCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCct                                                                <  2:6556741/151‑1 (MQ=255)
                                                                                          cgACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCctgctg                                                            >  1:1944580/1‑151 (MQ=255)
                                                                                           gACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTgg                                                           >  1:775691/1‑151 (MQ=255)
                                                                                              ttCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCa                                                        <  1:4050600/151‑1 (MQ=255)
                                                                                               tcggcgCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGACAACCACGAAGATTCTTCTGCATACCGGTGGACTGGCCCTGAGCCTGCTGGCCAc                                                       <  2:8050048/151‑1 (MQ=255)
                                                                                                     cgcTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGt                                                 >  2:5701988/1‑151 (MQ=255)
                                                                                                         cAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATg                                             <  1:4523678/151‑1 (MQ=255)
                                                                                                           gCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATggc                                           <  1:8118319/151‑1 (MQ=255)
                                                                                                           gCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATggc                                           >  2:8855328/1‑151 (MQ=255)
                                                                                                              tCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATggcggc                                        >  1:7363830/1‑151 (MQ=255)
                                                                                                              acGGCATGAACTTCTCATGGACGCCGGCCCGGTGCGACCCGGGCCGGCGCGCAGGACTTACCGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATggcggc                                        <  1:898582/150‑1 (MQ=255)
                                                                                                               cGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGACTGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGgcggcg                                       >  1:4332202/1‑151 (MQ=255)
                                                                                                               cGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGgcggcg                                       <  1:869514/151‑1 (MQ=255)
                                                                                                               cGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGgcggcg                                       >  2:5617472/1‑151 (MQ=255)
                                                                                                                 gCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGt                                     <  2:1500330/151‑1 (MQ=255)
                                                                                                                 gCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGt                                     >  1:2974855/1‑151 (MQ=255)
                                                                                                                     gAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATcc                                 >  2:10226046/1‑151 (MQ=255)
                                                                                                                     gAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGGACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATcc                                 <  1:518890/151‑1 (MQ=255)
                                                                                                                      aaCTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATccg                                >  2:6311461/1‑151 (MQ=255)
                                                                                                                       aCTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATccgc                               <  2:4944578/151‑1 (MQ=255)
                                                                                                                        cTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATccgcc                              <  2:7168324/151‑1 (MQ=255)
                                                                                                                           cTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCcgac                           <  2:10547074/151‑1 (MQ=255)
                                                                                                                              aTGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAg                        <  1:9272701/151‑1 (MQ=255)
                                                                                                                                    gCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAGGCGGcc                  >  2:10912536/1‑151 (MQ=255)
                                                                                                                                    gCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAGGCGGcc                  <  1:8898174/151‑1 (MQ=255)
                                                                                                                                       ggCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGGCGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAGGCGGCCAAg               >  1:2269651/1‑151 (MQ=255)
                                                                                                                                       ggCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAGGCGGCCAAg               >  2:5517097/1‑151 (MQ=255)
                                                                                                                                          ccGGTGCGACCCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCAGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAGGCGGCCAAGcac            <  1:1635557/151‑3 (MQ=255)
                                                                                                                                            ggTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAGGCGGCCAAGCTCgg          >  2:697830/1‑151 (MQ=255)
                                                                                                                                            ggTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAGGCGGCCAAGCTCgg          >  2:694544/1‑151 (MQ=255)
                                                                                                                                              tGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAGGCGGCCAAGCTCGGCa        <  2:7923859/151‑1 (MQ=255)
                                                                                                                                              tGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAGGCGGCCAAGCTCGGCa        <  2:11167754/151‑1 (MQ=255)
                                                                                                                                                 gaaCCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAGGCGGCCACGCTCGGCACCa     >  2:1264597/2‑151 (MQ=255)
                                                                                                                                                 gaaCCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAGGCGGCCAAGCTCGGCACCa     >  2:3342202/2‑151 (MQ=255)
                                                                                                                                                    ccGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAGGCGGCCAAGCTCGGCACCAGcc  >  2:5099682/1‑151 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGGGGTGGACGCGGAGTACCAGAACACCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTNNNANCNGGNCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAGGCGGCCAAGCTCGGCACCAGCC  >  NC_008463/1128013‑1128311

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