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.5% 1554.3 / NA 421intergenic (+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 (205/214);  minor base A (1/1);  total (206/215)
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 423 total aligned reads displayed.

GGGGGTGGACGCGGAGTACCAGAACACCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTNNNANCNGGNCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAGGCGGCCAAGCTCGGCACCAGCC  >  NC_008463/1128013‑1128311
                                                                                                                                                     |                                                                                                                                                     
gggggTGGACGCGGAGTACCAGAACACCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCg                                                                                                                                                      <  1:11238035/151‑1 (MQ=255)
gggggTGGACGCGGAGTACCAGAACACCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCg                                                                                                                                                      <  1:3865976/151‑1 (MQ=255)
  gggTGGACGCGGAGTACCAGAACACCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGACCCggg                                                                                                                                                    >  2:4076951/1‑150 (MQ=255)
   ggTGGACGCGGAGTACCAGAACACCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGc                                                                                                                                                   >  2:7065713/1‑151 (MQ=255)
   ggTGGACGCGGAGTACCAGAACACCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGc                                                                                                                                                   <  1:2246410/151‑1 (MQ=255)
      ggACGCGGAGTACCAGAACACCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCgg                                                                                                                                                <  1:6501801/151‑1 (MQ=255)
        aCGCGGAGTACCAGAACACCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCg                                                                                                                                              <  2:7731274/151‑1 (MQ=255)
            ggAGTACCAGAACACCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGCCTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGTTGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCa                                                                                                                                          >  2:10486455/1‑151 (MQ=255)
                 aCCAGAACCCCTACACCGCGAGCCACTCCTACACCAACTTCTACGACGGCAAGTACACCCCCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACt                                                                                                                                     <  1:10335941/151‑1 (MQ=255)
                   cAGAACACCTACCCCGCGAGCCTCTCCTACACCACCTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTa                                                                                                                                   <  1:6714611/151‑1 (MQ=255)
                         acCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAgg                                                                                                                             <  2:112579/151‑1 (MQ=255)
                         acCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAgg                                                                                                                             <  1:8214413/151‑1 (MQ=255)
                         acCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAgg                                                                                                                             >  2:9304685/1‑151 (MQ=255)
                         acCTACACCGCGAGCCTCTCCTACACCAACTACGTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAgg                                                                                                                             <  2:7655434/151‑1 (MQ=255)
                          ccTACCCCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGga                                                                                                                            <  2:2478410/151‑1 (MQ=255)
                          ccTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGga                                                                                                                            >  2:389214/1‑151 (MQ=255)
                          ccTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGga                                                                                                                            >  2:3579311/1‑151 (MQ=255)
                           cTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGgac                                                                                                                           <  1:10020512/151‑1 (MQ=255)
                            tACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGTTGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGgacg                                                                                                                          >  2:3761162/1‑151 (MQ=255)
                            tACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGgacg                                                                                                                          <  1:4224830/151‑1 (MQ=255)
                            cacacCGCGAGCCTCTCCTACACCCACGCCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGGGCGAACCGGGCCGGCGAGCAGGACTTAACGAGgacg                                                                                                                          <  2:8813013/150‑1 (MQ=255)
                             acacCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACATCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGgacga                                                                                                                         >  2:8617004/1‑151 (MQ=255)
                                 cgcgAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAg                                                                                                                     >  2:1647272/1‑151 (MQ=255)
                                  gcgAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGc                                                                                                                    >  2:661456/1‑151 (MQ=255)
                                  gcgAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGc                                                                                                                    <  1:1252004/151‑1 (MQ=255)
                                  gcgAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACATCGTCGCGCTCAGCTTCGGCCTGAAATTCTGATGTACGCCGGCCCGGGGGGCACCCGGCCGGCGAGCAGGACATAACGAGGACGACGAGc                                                                                                                    >  1:54104/1‑151 (MQ=255)
                                     aGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCCTGAACTTCTGATGTACGCCGGCCCCGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGaa                                                                                                                 >  2:4740030/1‑151 (MQ=255)
                                        ctctCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCGTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGa                                                                                                              >  1:10726820/1‑151 (MQ=255)
                                        ctctCCTACACCAACGTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGa                                                                                                              <  2:11023233/151‑1 (MQ=255)
                                            ccTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGGACTTCTGATGTACGCCGGCCCCGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAAAGAAaac                                                                                                          >  1:9711143/1‑151 (MQ=255)
                                            ccTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAaac                                                                                                          <  1:8012548/151‑1 (MQ=255)
                                                  ccAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGaa                                                                                                    >  2:6436190/1‑151 (MQ=255)
                                                    aaCTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGGATGTACTTCTGATGTAAGCCGGGCCGGTGCGAACCGGGCCGGCGAGCAGGTCTTAACGAGGACGACGTGCGAATGTAAACAaagaagg                                                                                                  >  1:189609/1‑150 (MQ=255)
                                                    aaCTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCAGGAACTTCTGATGGACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACCACGAAGa                                                                                                  >  2:8114029/1‑151 (MQ=255)
                                                     aCTTCGACGACGGCAATTACACCACCGTGGACTATCGCTACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCGGGTGCTAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGAt                                                                                                 <  1:2978804/151‑1 (MQ=255)
                                                         cttcGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGAttctt                                                                                             <  2:2868015/151‑1 (MQ=255)
                                                             gACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGc                                                                                         <  1:697135/151‑1 (MQ=255)
                                                              aCGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATTTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAAAGAAAACAACGAAGATTCTTCTGCa                                                                                        >  1:5238203/1‑151 (MQ=255)
                                                               cGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCAt                                                                                       <  1:9806605/151‑1 (MQ=255)
                                                               cGGCAAGTACACCACCGTGGACGATCGCGACATCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCAt                                                                                       <  1:3982842/151‑1 (MQ=255)
                                                                ggCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATa                                                                                      >  2:14163/1‑151 (MQ=255)
                                                                ggCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATa                                                                                      >  2:428565/1‑151 (MQ=255)
                                                                 gCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATAc                                                                                     <  2:6687986/151‑1 (MQ=255)
                                                                  cAAGTTCACCACCGTGGACGATCGCGACTTCGTCGCGCGCAGCTTCGGCATGAACTTAAGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAAAGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATAcc                                                                                    >  2:7542311/1‑151 (MQ=255)
                                                                      tACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGgtg                                                                                <  1:9143217/151‑1 (MQ=255)
                                                                      tACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGgtg                                                                                <  1:3513507/151‑1 (MQ=255)
                                                                       acacCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGgtgt                                                                               <  1:6316323/151‑1 (MQ=255)
                                                                       acacCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATAACGgtgt                                                                               >  2:5578566/1‑151 (MQ=255)
                                                                        caccacCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGACCTTCTGCTGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTa                                                                              <  1:759793/151‑1 (MQ=255)
                                                                        caccacCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGGACGCCGGCCCGGTGCGACCCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTa                                                                              >  1:10897336/1‑151 (MQ=255)
                                                                         accaccGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTAc                                                                             >  2:7923583/1‑151 (MQ=255)
                                                                         accaccGTGGACGATCGCGACTTAGTCGCGCTCAGCTTCGGCATTAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGACTGAAAACAACGAAGATTCGTCTGCATACCGGTGTag                                                                             >  2:10891933/1‑150 (MQ=255)
                                                                           caccGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTg                                                                           <  1:1396625/151‑1 (MQ=255)
                                                                             ccGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGTAAACAACGAAGATTCTTCAGCATACCGGTGTACTGGc                                                                         >  2:11146504/1‑151 (MQ=255)
                                                                             ccGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGc                                                                         <  1:10592412/151‑1 (MQ=255)
                                                                              cGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGcc                                                                        <  1:55536/151‑1 (MQ=255)
                                                                                tGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCt                                                                      <  2:9851039/151‑1 (MQ=255)
                                                                                     gATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGcc                                                                 <  1:7427227/151‑1 (MQ=255)
                                                                                       tCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAaacaac                                                                                                       >  2:9163149/1‑111 (MQ=255)
                                                                                       tCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCctg                                                               >  1:5804467/1‑151 (MQ=255)
                                                                                        cgcgACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCctgc                                                              >  2:9013709/1‑151 (MQ=255)
                                                                                          cgACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCctgctg                                                            >  2:8417595/1‑151 (MQ=255)
                                                                                          cgACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCctgctg                                                            <  2:511046/151‑1 (MQ=255)
                                                                                           gACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGACCCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTgg                                                           >  2:3384556/1‑151 (MQ=255)
                                                                                            aCTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGGGTACTGGCCCTGAGCCTGCTGGc                                                          <  2:5959310/151‑1 (MQ=255)
                                                                                                cgtcgCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCAcc                                                      <  2:1899643/151‑1 (MQ=255)
                                                                                                 gtcgCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCAccc                                                     >  2:4945386/1‑151 (MQ=255)
                                                                                                   cgcgCTCAGCTTAGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAg                                                   >  1:4746862/1‑151 (MQ=255)
                                                                                                       cTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGa                                               >  1:3312036/1‑151 (MQ=255)
                                                                                                       cTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGa                                               <  2:6911521/151‑1 (MQ=255)
                                                                                                         cAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATg                                             >  2:8228847/1‑151 (MQ=255)
                                                                                                          aGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATgg                                            >  1:9909762/1‑151 (MQ=255)
                                                                                                            cTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATggcg                                          >  2:20625/1‑151 (MQ=255)
                                                                                                              tCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATggcggc                                        <  1:2146830/151‑1 (MQ=255)
                                                                                                                ggCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGcggcgg                                      >  1:238619/1‑151 (MQ=255)
                                                                                                                 gCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGt                                     <  1:11089705/151‑1 (MQ=255)
                                                                                                                 gCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGt                                     <  2:9616131/151‑1 (MQ=255)
                                                                                                                   aTGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTAt                                   >  1:7351358/1‑151 (MQ=255)
                                                                                                                   aTGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTAt                                   <  1:6116958/151‑1 (MQ=255)
                                                                                                                     gAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATcc                                 <  1:8419398/151‑1 (MQ=255)
                                                                                                                           cTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCcgac                           <  1:9529123/151‑1 (MQ=255)
                                                                                                                             gATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAACCAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCcgacga                         <  1:2541661/151‑1 (MQ=255)
                                                                                                                             gATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCcgacga                         <  1:785053/151‑1 (MQ=255)
                                                                                                                             gATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCcgacga                         <  1:5062554/151‑1 (MQ=255)
                                                                                                                              acgaaCGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGCCGACGAGCGTATGAAACCAACGAAGATTATTCTTCGTACCGGTGTACTGGCCCTGAGCCTGCTGGCCTACCAGGTGATGGCGGCGGTATCCGCCGACGAg                        <  1:637465/147‑1 (MQ=255)
                                                                                                                               tGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGAGTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTCTCCGCCGACGAgg                       >  2:5760939/1‑151 (MQ=255)
                                                                                                                                   cGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAggcggc                   <  2:2599373/151‑1 (MQ=255)
                                                                                                                                    gCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTACCGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAGGCGGcc                  <  2:7188208/151‑1 (MQ=255)
                                                                                                                                    gCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAATCAACGAAGATTCTTCTCCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAGGCGGcc                  >  1:8828892/1‑151 (MQ=255)
                                                                                                                                       ggCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACCACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAGGCGGCCAAg               >  2:2329354/1‑151 (MQ=255)
                                                                                                                                         cccGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAGGCGGCCAAGCt             >  2:1003220/1‑151 (MQ=255)
                                                                                                                                          ccGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAGGCGGCCAAGCTc            <  1:3761162/151‑1 (MQ=255)
                                                                                                                                           cGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAGGCGGCCAAGCTCg           >  1:7195616/1‑151 (MQ=255)
                                                                                                                                            ggTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAGGCGGCCAAGCTCgg          >  1:9376158/1‑151 (MQ=255)
                                                                                                                                            ggTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAGGCGGCCAAGCTCgg          >  1:8750056/1‑151 (MQ=255)
                                                                                                                                                 gaaCCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAGGCGGCCAAGCTCGGCACCa     <  2:7320099/150‑1 (MQ=255)
                                                                                                                                                 gaaCCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAGGCGGCCAAGCTCGGCACCa     <  1:9291861/150‑1 (MQ=255)
                                                                                                                                                 gaaCAGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATCACGGTGTACTGGCCCTGAGCCTGCTGGCAACCCAGGTGATGGCGGCGGTATCCGCCGACGAGGCGGCCAAGCTCGGCACCa     <  1:1714081/150‑1 (MQ=255)
                                                                                                                                                 gaCCCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGACTGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAGGCGGCCAAGCTCGGCACCa     >  2:5393459/2‑151 (MQ=255)
                                                                                                                                                  aaCCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCCCCCAGGTGATGGCGGCGGTATCCGCCGACGAGGCGGCCAAGCTCGGCACCAg    >  2:2804548/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: