Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 2,001,811 A→C Q66P (CAG→CCG)  pgl → 6‑phosphogluconolactonase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084632,001,8110AC99.8% 1524.7 / NA 397Q66P (CAG→CCG) pgl6‑phosphogluconolactonase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (173/223);  total (174/223)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 397 total aligned reads displayed.

GCAGCGCCGCCGAGCAGGCCCGCGGCCTGGCCGCCGAGGTCGCCGGCCGGTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGATCAGGTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTGGCGAGGCGGCGAAGGCGCGCTTCATCGGCCTCTACCAGC  >  NC_008463/2001667‑2001954
                                                                                                                                                |                                                                                                                                               
gCAGCGCCGCCGAGCAGGCCCGCGGCCTGGCCGCCGAGGTCGCCGGCCGGTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGa                                                                                                                                           >  1:7278128/1‑151 (MQ=255)
 cAGCGCCGCCGAGCAGGCCCGCGGCCTGGCCGCCGAGGTCGCCGGCCGGTTGCGCGAGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCAGCTGGAc                                                                                                                                          >  1:892147/1‑151 (MQ=255)
  aGCGCCGCCGAGCAGGCCCGCGGCCTGGCCGCCGAGGTCGCCGGCCGGTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACt                                                                                                                                         <  1:6206971/151‑1 (MQ=255)
  aGCGCCGCCGAGCAGGCCCGCGGCCTGGCCGCCGAGGTCGCCGGCCGGTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACt                                                                                                                                         >  2:2063687/1‑151 (MQ=255)
   gcgcCGCCGAGCAGGCCCGCGGCCTGGCCGCCGAGGTCGCCGGCCGGTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTg                                                                                                                                        <  2:6552617/151‑1 (MQ=255)
      ccgccgAGCAGGCCCGCGGCCTGGCCGCCGAGGTCGCCGGCCGGTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTc                                                                                                                                     <  1:2563528/151‑1 (MQ=255)
      ccgccgAGCAGGCCCGCGGCCTGCCCGCCGAGGTCGCCGGCCGGTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTc                                                                                                                                     >  2:5035792/1‑151 (MQ=255)
        gccgAGCAGGCCCGCGGCCTGGCCGCCGAGGTCGCCGGCCGGTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAGGCCTTGAGCGAGGAGCCGCTGGACTGGTcgc                                                                                                                                   >  1:3880693/1‑151 (MQ=255)
        gccgAGCAGGCCCGCGGCCTGGCCGCCGAGGTCGCCGGCCGGTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTcgc                                                                                                                                   >  2:5544672/1‑151 (MQ=255)
         ccgAGCAGGCCCGCGGCCTGGCCGCCGAGGTCGCCGGCCGGTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTcgcg                                                                                                                                  <  2:6968881/151‑1 (MQ=255)
            aGCAGGCCCGCGGCCTGGCCGCCGAGGTCGCCGGCCGGTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGAt                                                                                                                               <  1:998394/151‑1 (MQ=255)
              cAGGCCCGCGGCCTGGCCGCCGAGGTCGCCGGCCGGGTGCGCTCGGCGCGGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCa                                                                                                                             <  2:10124182/151‑1 (MQ=255)
               aGGCCCGCGGCCTGGCCGCCGAGGTCGCCGGCCGGTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCAc                                                                                                                            >  1:1091147/1‑151 (MQ=255)
                   ccGCGGCCTGGCCGCCGAGGTCGCCGGCCGGTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTc                                                                                                                        <  2:7599265/151‑1 (MQ=255)
                   ccGCGGCCTGGCCGCCGAGGTCGCCGGCCGGTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTc                                                                                                                        <  1:6743062/151‑1 (MQ=255)
                   ccGCGGCCTGGCCGCCGAGGTCGCCGGCCGGTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTc                                                                                                                        <  2:5386883/151‑1 (MQ=255)
                         ccTGGCCGCCGAGGTCGCCGGCCGGTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTg                                                                                                                  <  1:3047289/151‑1 (MQ=255)
                         ccTGGCCGCCGAGGTCGCCGGCCGGTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTg                                                                                                                  >  1:1148608/1‑151 (MQ=255)
                          cTGGCCGCCGAGGTCGCCGGCCGGTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTgg                                                                                                                 <  1:8099222/151‑1 (MQ=255)
                           tGGCCGCCGAGGTCGCCGGCCGGTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGc                                                                                                                >  1:791402/1‑151 (MQ=255)
                           tGGCCGCCGAGGTCGCCGGCCGGTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGc                                                                                                                <  2:2197760/151‑1 (MQ=255)
                           tGGCCGCCGAGGTCGCCGGCCGGTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGc                                                                                                                <  1:4569133/151‑1 (MQ=255)
                            ggCCGCCGAGGTCGCCGGCCGGTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGcc                                                                                                               <  1:5678121/151‑1 (MQ=255)
                               cgccgAGGTCGCCGGCCGGTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCcgac                                                                                                            >  2:2885304/1‑151 (MQ=255)
                                gccgAGGTCGCCGGCCGGTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCcgacg                                                                                                           <  2:4674500/151‑1 (MQ=255)
                                  cgAGGTCGCCGGCCGGTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAg                                                                                                         >  2:2699409/1‑151 (MQ=255)
                                  cgAGGTCGCCGGCCGGTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAg                                                                                                         <  2:9843540/151‑1 (MQ=255)
                                  cgAGGTCGCCGGCCGGTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAg                                                                                                         <  1:5081727/151‑1 (MQ=255)
                                    aGGTCGCCGGCCGGTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAgcg                                                                                                       <  1:9760787/151‑1 (MQ=255)
                                     ggTCGCCGGCCGGTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAgcgc                                                                                                      >  1:373657/1‑151 (MQ=255)
                                      gTCGCCGGCCGGTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAgcgcg                                                                                                     >  2:5819782/1‑150 (MQ=255)
                                           cggccggTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTg                                                                                                >  2:5665671/1‑151 (MQ=255)
                                            ggccggTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGc                                                                                               <  1:7359660/151‑1 (MQ=255)
                                            ggccggTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGc                                                                                               >  2:6932710/1‑151 (MQ=255)
                                             gccggTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGcc                                                                                              <  1:1978220/151‑1 (MQ=255)
                                                 gTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAg                                                                                          >  1:7768977/1‑151 (MQ=255)
                                                   tGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTc                                                                                        <  1:10477161/151‑1 (MQ=255)
                                                    gcgcTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCg                                                                                       <  1:10058301/151‑1 (MQ=255)
                                                    gcgcTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCGGGGTGCCGGAGTCg                                                                                       >  2:2525512/1‑151 (MQ=255)
                                                     cgcTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGc                                                                                      >  2:415757/1‑151 (MQ=255)
                                                     cgcTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGc                                                                                      <  2:3427101/151‑1 (MQ=255)
                                                        tCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATg                                                                                   <  2:3188644/151‑1 (MQ=255)
                                                         cGGCGCTTGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCAGGc                                                                                  <  1:8344459/151‑1 (MQ=255)
                                                           gcgcTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCg                                                                                <  1:6088381/151‑1 (MQ=255)
                                                             gcTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGAt                                                                              <  2:3363797/151‑1 (MQ=255)
                                                             gcTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGAt                                                                              >  2:1486273/1‑151 (MQ=255)
                                                              cTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATa                                                                             <  1:1341755/151‑1 (MQ=255)
                                                                  ccGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCaa                                                                         >  1:8626402/1‑151 (MQ=255)
                                                                  ccGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCaa                                                                         >  1:2289275/1‑151 (MQ=255)
                                                                  ccGAGCAGGGCCAGGCGCTGCTGGTGGGGTCCGGGGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGATGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCaa                                                                         <  1:7238655/151‑1 (MQ=255)
                                                                   cGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAAc                                                                        <  2:1162559/151‑1 (MQ=255)
                                                                   cGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAAc                                                                        <  1:10379018/151‑1 (MQ=255)
                                                                   cGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAAc                                                                        <  1:5160723/151‑1 (MQ=255)
                                                                   cGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCGGGGTGCCGGAGTCGCATGCCGATAGCAAc                                                                        >  1:9371544/1‑151 (MQ=255)
                                                                    gAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACg                                                                       <  2:6146165/151‑1 (MQ=255)
                                                                      gctgggCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGATCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGAGAGCAACGcc                                                                     <  2:1367439/148‑1 (MQ=255)
                                                                      gCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGcc                                                                     >  1:2093351/1‑151 (MQ=255)
                                                                        aGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCgg                                                                   >  2:4806094/1‑151 (MQ=255)
                                                                         gggCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGc                                                                  >  2:4907024/1‑151 (MQ=255)
                                                                          ggCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGcc                                                                 <  2:9047488/151‑1 (MQ=255)
                                                                            ccAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTg                                                               <  1:6230483/151‑1 (MQ=255)
                                                                              aGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCACCGCCGGCCTGGt                                                             >  2:8422017/1‑151 (MQ=255)
                                                                               ggCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGtt                                                            >  2:9184618/1‑151 (MQ=255)
                                                                                 cgcTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTcg                                                          <  1:2809919/151‑1 (MQ=255)
                                                                                      ctgGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCgccgcc                                                     >  2:7802885/1‑151 (MQ=255)
                                                                                       tggtggTGTCCGGTGGGCGCAGTCCGGGGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCa                                                    <  2:8657956/151‑1 (MQ=255)
                                                                                        ggtggtGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCAc                                                   <  1:4662134/151‑1 (MQ=255)
                                                                                         gtggtgTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCAcc                                                  <  1:3374250/151‑1 (MQ=255)
                                                                                          tggggTCCGGTGGGCGCTGGCCGGTGGCCGTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGCGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCt                                                 <  2:5885500/151‑1 (MQ=255)
                                                                                              gtCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTc                                             >  1:4302627/1‑151 (MQ=255)
                                                                                               tCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTcc                                            <  2:3024995/151‑1 (MQ=255)
                                                                                                 cGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGt                                          <  1:5311187/151‑1 (MQ=255)
                                                                                                 cGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGt                                          >  2:8078502/1‑151 (MQ=255)
                                                                                                 cGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGt                                          <  1:7544269/151‑1 (MQ=255)
                                                                                                  ggTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTg                                         >  1:10138115/1‑151 (MQ=255)
                                                                                                  ggTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTg                                         <  1:9500866/151‑1 (MQ=255)
                                                                                                  ggTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTg                                         >  1:6649855/1‑151 (MQ=255)
                                                                                                   gTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTgg                                        <  1:6857740/151‑1 (MQ=255)
                                                                                                   gTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTgg                                        >  1:3500487/1‑151 (MQ=255)
                                                                                                     gggCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTGGCg                                      >  1:5004638/1‑151 (MQ=255)
                                                                                                            gTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTGGCGAggcggc                               <  2:6464708/151‑1 (MQ=255)
                                                                                                            gTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTGGCGAggcggc                               <  2:3760360/151‑1 (MQ=255)
                                                                                                              ccGGTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTGGCGAGGCGGCGa                             >  1:4481598/1‑151 (MQ=255)
                                                                                                                ggTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTGGCGAGGCGGCGAAg                           <  1:584878/151‑1 (MQ=255)
                                                                                                                ggTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTGGCGAGGCGGCGAAg                           <  2:9649176/151‑1 (MQ=255)
                                                                                                                 gTGGCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTGGCGAGGCGGCGAAgg                          >  2:7191786/1‑151 (MQ=255)
                                                                                                                    gCCTTCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTGGCGAGGCGGCGAAGgcgc                       <  2:3895019/151‑1 (MQ=255)
                                                                                                                        tCCTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTGGCGAGGCGGCGAAGGCGCGCtt                   >  2:4283737/1‑151 (MQ=255)
                                                                                                                         ccTCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTGGCGAGGCGGCGAAGGCGCGCTTc                  <  1:8412509/151‑1 (MQ=255)
                                                                                                                           tCGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTGGCGAGGCGGCGAAGGCGCGCTTCAt                >  1:2670190/1‑151 (MQ=255)
                                                                                                                           acGAAGCCTAGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTGGCGAGGCGGCGAAGGCGCGCTTCAt                <  1:1614017/150‑1 (MQ=255)
                                                                                                                            cGAAGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTGGCGAGGCGGCGAAGGCGCGCTTCATc               >  1:596713/1‑151 (MQ=255)
                                                                                                                              aaGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTGGCGAGGCGGCGAAGGCGCGCTTCATCgg             <  2:6292472/151‑1 (MQ=255)
                                                                                                                              aaGCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTGGCGAGGCGGCGAAGGCGCGCTTCATCgg             >  1:3115610/1‑151 (MQ=255)
                                                                                                                                gCCTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTGGCGAGGCGGCGAAGGCGCGCTTCATCGGcc           >  2:2356889/1‑151 (MQ=255)
                                                                                                                                  cTTGAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTGGCGAGGCGGCGAAGGCGCGCTTCATCGGCctc         >  2:454947/1‑151 (MQ=255)
                                                                                                                                   ttGAGCGTGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGAGAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTGGCGAGGCGGCGAAGGCGCGCTTCATCGGCctct        <  2:7390581/151‑1 (MQ=255)
                                                                                                                                     gAGCGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTGGCGAGGCGGCGAAGGCGCGCTTCATCGGCCTCTAc      <  2:10310921/151‑1 (MQ=255)
                                                                                                                                        cGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTGGCGAGGCGGCGAAGGCGCGCTTCATCGGCCTCTACCAg   >  1:747073/1‑151 (MQ=255)
                                                                                                                                        cGAGGAGCCGCTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTGGCGAGGCGGCGAAGGCGCGCTTCATCGGCCTCTACCAg   >  1:1490954/1‑151 (MQ=255)
                                                                                                                                                |                                                                                                                                               
GCAGCGCCGCCGAGCAGGCCCGCGGCCTGGCCGCCGAGGTCGCCGGCCGGTTGCGCTCGGCGCTGGCCGAGCAGGGCCAGGCGCTGCTGGTGGTGTCCGGTGGGCGCAGTCCGGTGGCCTTCCTCGAAGCCTTGAGCGAGGATCAGGTGGACTGGTCGCGGATCACCCTCAGCCTGGCCGACGAGCGCTGGGTGCCGGAGTCGCATGCCGATAGCAACGCCGGCCTGGTTCGCCGCCACCTGCTCCGTGGCGAGGCGGCGAAGGCGCGCTTCATCGGCCTCTACCAGC  >  NC_008463/2001667‑2001954

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