Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 1,128,165 N→G intergenic (+26/‑35) 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,1650NG99.8% 1598.2 / ‑7.3 423intergenic (+26/‑35)PA14_RS05280/PA14_RS05285DUF1302 domain‑containing protein/DUF1329 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base N (0/0);  major base G (208/214);  minor base T (0/1);  total (208/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 = 9.96e-01

Only 100 of 431 total aligned reads displayed.

GGTGGACGCGGAGTACCAGAACACCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTNNNANCNGGNCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAGGCGGCCAAGCTCGGCACCAGCCTGA  >  NC_008463/1128016‑1128314
                                                                                                                                                     |                                                                                                                                                     
ggTGGACGCGGAGTACCAGACCACCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGc                                                                                                                                                      <  1:5659590/151‑1 (MQ=255)
 gTGGACGCGGAGTACCAGAACCCCTACACCGCGAGCCTCTCCTACACCAACTTCCTCGACGGCAAGTACACCCCCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGcc                                                                                                                                                     <  1:8706998/151‑1 (MQ=255)
  tGGACGCGGAGTACCAGAACACCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCg                                                                                                                                                    <  1:10686557/151‑1 (MQ=255)
  tGGACGCGGAGTACCAGAACACCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCg                                                                                                                                                    <  1:4895425/151‑1 (MQ=255)
  tGGACGCGGAGTACCAGAACACCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCg                                                                                                                                                    <  2:930558/151‑1 (MQ=255)
   ggACGCGGAGTACCAGAACACCTACACCGCGAGCCTCTCCTACACCAACTTCGTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGGGCGAACCGGGCCgg                                                                                                                                                   >  1:4756590/1‑151 (MQ=255)
        cgGAGTACCAGAACCCCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGc                                                                                                                                              <  1:8640161/151‑1 (MQ=255)
        cgGAGTACCAGAACACCTACCCCGCGAGCCTCTCCTACACCACCTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGc                                                                                                                                              <  2:6347054/151‑1 (MQ=255)
         ggAGTACCAGAACACCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCa                                                                                                                                             <  2:2808001/151‑1 (MQ=255)
          taGTACCATAACACCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAg                                                                                                                                            <  1:1981647/150‑1 (MQ=255)
            gTACCAGAACACCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGa                                                                                                                                          >  2:1364052/1‑151 (MQ=255)
              aCCAGAACACCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCGTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACt                                                                                                                                        <  2:4815771/151‑1 (MQ=255)
                cAGAACACCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGTTGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTa                                                                                                                                      >  1:461932/1‑151 (MQ=255)
                  gAACACCTACACCGCGCGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAAc                                                                                                                                    <  1:10600315/151‑1 (MQ=255)
                  gAACACCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATTAACTTCTGATGAACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAAc                                                                                                                                    <  1:381815/151‑1 (MQ=255)
                  gAACACCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAAc                                                                                                                                    >  2:379708/1‑151 (MQ=255)
                   cacacCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACg                                                                                                                                   <  1:5448096/150‑1 (MQ=255)
                    acacCTACACCGCGCGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATTACCTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGa                                                                                                                                  <  1:8077604/151‑1 (MQ=255)
                    acacCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGa                                                                                                                                  >  2:9724691/1‑151 (MQ=255)
                     ccccTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAg                                                                                                                                 <  2:4062800/149‑1 (MQ=255)
                     cacCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAg                                                                                                                                 <  2:56019/151‑1 (MQ=255)
                      acCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAgg                                                                                                                                >  2:228416/1‑151 (MQ=255)
                      acCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAgg                                                                                                                                >  1:7990352/1‑151 (MQ=255)
                       ccTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGAACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGga                                                                                                                               <  1:10085010/151‑1 (MQ=255)
                         tACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGgacg                                                                                                                             <  2:10370520/151‑1 (MQ=255)
                           cacCGCGAGCCTCTCCTACACCAACTTCTTCGACGTCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGgacgac                                                                                                                           <  2:1428728/151‑1 (MQ=255)
                           cacCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGgacgac                                                                                                                           >  2:7412531/1‑151 (MQ=255)
                           cacCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTAGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGgacgac                                                                                                                           >  2:4469243/1‑151 (MQ=255)
                              cgcgAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAg                                                                                                                        <  1:2179072/151‑1 (MQ=255)
                               gcgACCCCCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGGGGACGATCGCGCCTTCGGCGCGCACAGCTACGGCATGAACTCCTGAAGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACATCACGAGGACGACGAGc                                                                                                                       <  2:6710494/151‑1 (MQ=255)
                                cgAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCg                                                                                                                      <  2:8016793/151‑1 (MQ=255)
                                cgAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCg                                                                                                                      >  2:9439481/1‑151 (MQ=255)
                                cgAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACATCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGTGCg                                                                                                                      >  1:5137883/1‑151 (MQ=255)
                                  aGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCGTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACCAGCGaa                                                                                                                    >  1:801127/1‑151 (MQ=255)
                                   gCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAAt                                                                                                                   >  2:6896064/1‑151 (MQ=255)
                                      tctcCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCCGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGaa                                                                                                                >  1:6687182/1‑151 (MQ=255)
                                      tctcCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGaa                                                                                                                <  1:8369834/151‑1 (MQ=255)
                                        tcCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCCGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGaaaa                                                                                                              >  1:1346667/1‑151 (MQ=255)
                                         ccTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCATCGGCAAGAACTTCTGATGTACGCCGGCCCGGTTCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAaac                                                                                                             >  2:9066549/1‑151 (MQ=255)
                                             cacCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACg                                                                                                         >  2:212856/1‑151 (MQ=255)
                                               ccAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTACTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGaa                                                                                                       <  2:3319634/151‑1 (MQ=255)
                                                cAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAg                                                                                                      >  2:8573473/1‑151 (MQ=255)
                                                cAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAg                                                                                                      <  1:10166708/151‑1 (MQ=255)
                                                  aCTTCTTCGACGGCACGTACCCCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCGGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGAt                                                                                                    <  2:7511004/151‑1 (MQ=255)
                                                    ttcttcGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGAttc                                                                                                  >  2:600016/1‑151 (MQ=255)
                                                    ttcttcGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACATAACGAGGACGACGAGCGAATGAAAACAACGAAAAtta                                                                                                  >  1:1239471/1‑150 (MQ=255)
                                                       ttcGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGAttcttc                                                                                               >  2:6178527/1‑151 (MQ=255)
                                                       ttcGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACATCTGATTTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAATAttcttc                                                                                               >  1:1514606/1‑151 (MQ=255)
                                                        tcGACGGCAAGTACACCACCGTGGACGATCGCGACTACGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATtcttct                                                                                              <  2:8272291/151‑1 (MQ=255)
                                                         cGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTg                                                                                             <  1:5591780/151‑1 (MQ=255)
                                                          gACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGc                                                                                            >  2:8399781/1‑151 (MQ=255)
                                                           aCGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCCGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTATGCa                                                                                           >  1:8863479/1‑151 (MQ=255)
                                                             ggCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGACCCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATa                                                                                         <  1:7188798/151‑1 (MQ=255)
                                                               cAAGTACCCCACCGTGGACGATCGCGACTTCGTCGCGCACAGCTTCGGCCTGAACGTCTGATGTACGCCGGCCCGGTGCGACCCGGGCCGGCGAGCAGGACTTACCGAGGACGACGAGCGAATGAAAACAACGCAGATTCTTCTGCATAcc                                                                                       <  2:2558042/151‑1 (MQ=255)
                                                                 aGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTGTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTGTGCATAACgg                                                                                     >  1:7557560/1‑151 (MQ=255)
                                                                  gTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGgt                                                                                    <  1:580704/151‑1 (MQ=255)
                                                                     caccacCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTa                                                                                 >  2:7936455/1‑151 (MQ=255)
                                                                        caccGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGTCAACAACGAAGATTCTTCTGCATACCGGTGTACTg                                                                              >  2:5223779/1‑151 (MQ=255)
                                                                          ccGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGc                                                                            <  1:9806182/151‑1 (MQ=255)
                                                                            gTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGTACGTCTGATGTAAGCCGGCCCGGTGCGAACCGGGCCGGCGCGCAGGACTTAACGAGGACGACGAGCGACTGAAAACAACGAAGTTTCTTCTGCATACCGGTGTACTGGccc                                                                          >  2:3454362/1‑151 (MQ=255)
                                                                             tGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATAACGGTGTACTGGCCCt                                                                         >  1:6095965/1‑151 (MQ=255)
                                                                               gACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGa                                                                       <  2:5907448/151‑1 (MQ=255)
                                                                                 cGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGc                                                                     >  2:3155230/1‑151 (MQ=255)
                                                                                  gATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGcc                                                                    >  2:1296370/1‑151 (MQ=255)
                                                                                  gATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGcc                                                                    >  1:5942980/1‑151 (MQ=255)
                                                                                   ctCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGCCGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCct                                                                   <  2:7124309/150‑1 (MQ=255)
                                                                                   aTCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACACCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCct                                                                   <  1:7260624/151‑1 (MQ=255)
                                                                                    tCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCctg                                                                  <  1:777014/151‑1 (MQ=255)
                                                                                      gcgACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCctgct                                                                <  2:4243680/151‑1 (MQ=255)
                                                                                          cTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGCCGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGcc                                                            >  1:2718035/1‑151 (MQ=255)
                                                                                                 gcgcTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAgg                                                     >  1:7441798/1‑151 (MQ=255)
                                                                                                     tCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGAt                                                 >  1:4452653/1‑151 (MQ=255)
                                                                                                      cAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATg                                                <  1:648529/151‑1 (MQ=255)
                                                                                                       aGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATgg                                               >  2:1069246/1‑151 (MQ=255)
                                                                                                       aGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATgg                                               >  1:9446616/1‑151 (MQ=255)
                                                                                                        gCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATggc                                              >  1:2304462/1‑151 (MQ=255)
                                                                                                         cTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATggcg                                             >  2:921752/1‑151 (MQ=255)
                                                                                                            cGGCATGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGgcggcg                                          <  1:9928409/151‑1 (MQ=255)
                                                                                                                aTGAACTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTAt                                      >  1:7010908/1‑151 (MQ=255)
                                                                                                                   aaCTTCTGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATccg                                   <  1:6879279/151‑1 (MQ=255)
                                                                                                                         tGATGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACATAACGAGGACGACGAGCGAATGAAAACAACGAATATTCTTCTTCATACCGGTGTACTTGCCCTGATCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCcgacg                             >  1:9052788/1‑151 (MQ=255)
                                                                                                                            tGTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAgg                          <  1:3708985/151‑1 (MQ=255)
                                                                                                                             gTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAggc                         <  1:3618313/151‑1 (MQ=255)
                                                                                                                             gTACGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTGAACGCGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGAGGGCGGCGGTATCCGCCGACGAggc                         <  2:1006108/151‑1 (MQ=255)
                                                                                                                               aCGCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAggcgg                       >  2:6547384/1‑151 (MQ=255)
                                                                                                                                 gCCGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAGGCGGcc                     >  2:5335675/1‑151 (MQ=255)
                                                                                                                                   cGGCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAGGCGGCCaa                   <  1:3126809/151‑1 (MQ=255)
                                                                                                                                    ggCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAGGCGGCCAAg                  <  1:6249817/151‑1 (MQ=255)
                                                                                                                                    ggCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAGGCGGCCAAg                  >  2:4658651/1‑151 (MQ=255)
                                                                                                                                    ggCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAGGCGGCCAAg                  <  1:2577793/151‑1 (MQ=255)
                                                                                                                                     gCCCGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAGGCGGCCAAGc                 <  2:8095700/151‑1 (MQ=255)
                                                                                                                                       ccGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAGGCGGCCAAGCTc               <  1:4019834/151‑1 (MQ=255)
                                                                                                                                        cGGTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAGGCGGCCAAGCTCg              >  1:4261144/1‑151 (MQ=255)
                                                                                                                                         ggTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAGGCGGCCAAGCTCgg             >  2:2445906/1‑151 (MQ=255)
                                                                                                                                         ggTGCGAACCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAGGCGGCCAAGCTCgg             >  2:5092759/1‑151 (MQ=255)
                                                                                                                                             cgaaCCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAGGCGGCCAAGCTCGGCAcc         >  2:1771315/3‑151 (MQ=255)
                                                                                                                                               aaCCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAATACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAGGCGGCCAAGCTCGGCACCAg       <  2:2003823/151‑1 (MQ=255)
                                                                                                                                               aaCCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAGGCGGCCAAGCTCGGCACCAg       >  2:8860419/1‑151 (MQ=255)
                                                                                                                                               aaCCGGGCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAGGCGGCCAAGCTCGGCACCAg       <  1:8476824/151‑1 (MQ=255)
                                                                                                                                                  cgggCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATggcggc                                           >  2:1562534/2‑112 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGTGGACGCGGAGTACCAGAACACCTACACCGCGAGCCTCTCCTACACCAACTTCTTCGACGGCAAGTACACCACCGTGGACGATCGCGACTTCGTCGCGCTCAGCTTCGGCATGAACTTCTGATGTACGCCGGCCCGGTNNNANCNGGNCCGGCGAGCAGGACTTAACGAGGACGACGAGCGAATGAAAACAACGAAGATTCTTCTGCATACCGGTGTACTGGCCCTGAGCCTGCTGGCCACCCAGGTGATGGCGGCGGTATCCGCCGACGAGGCGGCCAAGCTCGGCACCAGCCTGA  >  NC_008463/1128016‑1128314

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: ‑