Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 3,051,413 T→C K369E (AAG→GAG)  PA14_RS13990 ← SfnB family sulfur acquisition oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084633,051,4130TC99.8% 2036.2 / NA 509K369E (AAG→GAG) PA14_RS13990SfnB family sulfur acquisition oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (271/237);  minor base A (1/0);  total (272/237)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 509 total aligned reads displayed.

GCTCATCGGGCGCTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTTGAGGTTGAACTTGGCGAGCGTGGCGCGGCTGCCGGCCAGCTTGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCACCGGCCGCGCGGCCACTTCGT  >  NC_008463/3051266‑3051559
                                                                                                                                                   |                                                                                                                                                  
gCTCATCGGGCGCTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCgg                                                                                                                                                  <  1:10306776/150‑1 (MQ=255)
gCTCATCGGGCGCTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCgg                                                                                                                                                  <  2:6736007/150‑1 (MQ=255)
 cTCATCGGGCGCTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGc                                                                                                                                                 >  1:7246015/1‑150 (MQ=255)
    aTCGGGCGCTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAg                                                                                                                                              >  2:611037/1‑150 (MQ=255)
     tCGGGCGCTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGc                                                                                                                                             >  2:10305373/1‑150 (MQ=255)
      cGGGCGCTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCg                                                                                                                                            <  2:1914352/150‑1 (MQ=255)
       gggCGCTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGt                                                                                                                                           <  1:5871220/150‑1 (MQ=255)
           gcTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGgcg                                                                                                                                       >  1:1754941/1‑150 (MQ=255)
           gcTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGgcg                                                                                                                                       <  1:5502136/150‑1 (MQ=255)
            ctcctcAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGgcgc                                                                                                                                      <  2:2108883/150‑1 (MQ=255)
            ctcctcAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGgcgc                                                                                                                                      <  2:783042/150‑1 (MQ=255)
            ctcctcAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGgcgc                                                                                                                                      <  1:4347936/150‑1 (MQ=255)
              cctcAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCgg                                                                                                                                    >  1:6110567/1‑150 (MQ=255)
              cctcAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCgg                                                                                                                                    >  2:1486247/1‑150 (MQ=255)
              cctcAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCgg                                                                                                                                    <  1:3910366/150‑1 (MQ=255)
               ctcAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGc                                                                                                                                   >  2:2948346/1‑150 (MQ=255)
                  aGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGc                                                                                                                                >  1:3756029/1‑150 (MQ=255)
                  aGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGc                                                                                                                                <  2:4570882/150‑1 (MQ=255)
                    aTCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCg                                                                                                                              >  1:1332136/1‑150 (MQ=255)
                     tCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCgg                                                                                                                             <  2:6365438/150‑1 (MQ=255)
                       cATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGcc                                                                                                                           <  1:9551905/150‑1 (MQ=255)
                       cATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGcc                                                                                                                           <  2:10601380/150‑1 (MQ=255)
                        aTGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCa                                                                                                                          <  1:7761254/150‑1 (MQ=255)
                         tGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAg                                                                                                                         >  1:1109312/1‑150 (MQ=255)
                          gAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGc                                                                                                                        <  2:4127707/150‑1 (MQ=255)
                           aGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCt                                                                                                                       <  2:3885629/150‑1 (MQ=255)
                            gTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCtc                                                                                                                      >  1:1043454/1‑149 (MQ=255)
                               aCGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGaa                                                                                                                   >  1:10031775/1‑150 (MQ=255)
                                cgcgCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCCGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAag                                                                                                                  <  1:4770512/150‑1 (MQ=255)
                                cgcgCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAag                                                                                                                  <  2:10436625/150‑1 (MQ=255)
                                 gcgcCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAaga                                                                                                                 >  1:9849835/1‑150 (MQ=255)
                                 gcgcCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAaga                                                                                                                 <  1:7180138/150‑1 (MQ=255)
                                 gcgcCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAaga                                                                                                                 <  1:167468/150‑1 (MQ=255)
                                  cgcCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAagag                                                                                                                <  1:2802782/150‑1 (MQ=255)
                                    ccGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCt                                                                                                              >  2:8014959/1‑150 (MQ=255)
                                    ccGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCt                                                                                                              <  1:2498854/150‑1 (MQ=255)
                                      ggACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTtc                                                                                                            >  1:7597971/1‑150 (MQ=255)
                                       gACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTtct                                                                                                           <  1:200877/150‑1 (MQ=255)
                                         cGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCg                                                                                                         >  1:3664811/1‑150 (MQ=255)
                                            gcgcCGGTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTg                                                                                                      <  2:8613965/150‑1 (MQ=255)
                                             cgcCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTgg                                                                                                     <  1:9939527/150‑1 (MQ=255)
                                             cgcCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTgg                                                                                                     >  1:7983912/1‑150 (MQ=255)
                                               ccGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTggcg                                                                                                   >  1:466059/1‑150 (MQ=255)
                                                  ttCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGGACCGGGTCGTGCAGTGTGTGGGTGCGGGCGTTGCGCCAGAGCCGGTCGAGGTTGAACTCGGCGAGCGTTGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGcc                                                                                                <  2:1097275/150‑1 (MQ=25)
                                                       gtagtaGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCt                                                                                           <  2:5616831/150‑1 (MQ=255)
                                                         agtagtCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGa                                                                                         <  1:5952839/150‑1 (MQ=255)
                                                          gtagtCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGAt                                                                                        >  1:6701484/1‑150 (MQ=255)
                                                           tagtCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGAtc                                                                                       >  2:3481759/1‑150 (MQ=255)
                                                            agtCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGAtct                                                                                      >  2:9782028/1‑150 (MQ=255)
                                                            agtCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCGCGAAGAGCTTCTCGCTGGCGGCCAGGCTGAtct                                                                                      >  1:7303233/1‑150 (MQ=255)
                                                              tCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCg                                                                                    >  2:4713542/1‑150 (MQ=255)
                                                               cGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCgg                                                                                   >  2:9610449/1‑150 (MQ=255)
                                                                gCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGt                                                                                  >  1:9349765/1‑150 (MQ=255)
                                                                gCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGAGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCCGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGt                                                                                  >  1:6805129/1‑150 (MQ=38)
                                                                   gACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCt                                                                               >  2:10164820/1‑150 (MQ=255)
                                                                       gcgTGGACCTTCCAGCGTACCGGGTCGTGCCGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGc                                                                           >  1:2209334/1‑150 (MQ=255)
                                                                       gcgTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGc                                                                           <  1:5679114/150‑1 (MQ=255)
                                                                        cgTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCa                                                                          <  1:8146352/150‑1 (MQ=255)
                                                                         gTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCAc                                                                         >  2:9296025/1‑150 (MQ=255)
                                                                          tGGACCGTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCAcc                                                                        <  2:283000/150‑1 (MQ=255)
                                                                            gACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCtt                                                                      >  2:4314738/1‑150 (MQ=255)
                                                                               cTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGc                                                                   >  1:3252601/1‑150 (MQ=255)
                                                                               cTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGc                                                                   >  1:8972246/1‑150 (MQ=255)
                                                                                ttCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGcc                                                                  <  2:1244526/150‑1 (MQ=255)
                                                                                 tCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCt                                                                 <  1:10468091/150‑1 (MQ=255)
                                                                                  ccAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTc                                                                <  2:8306807/150‑1 (MQ=255)
                                                                                   cAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCg                                                               >  1:6864395/1‑150 (MQ=255)
                                                                                     gCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGc                                                             >  1:531998/1‑150 (MQ=255)
                                                                                       gTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCTCCTTGGCCTCGGCCa                                                           >  2:7491797/1‑150 (MQ=255)
                                                                                       gTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCa                                                           <  1:4186336/150‑1 (MQ=255)
                                                                                         aCCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACg                                                         >  2:5119312/1‑150 (MQ=255)
                                                                                         aCCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACg                                                         <  2:7360862/150‑1 (MQ=255)
                                                                                          ccGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACgg                                                        <  1:5115419/150‑1 (MQ=255)
                                                                                             ggTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGa                                                     <  1:6161911/150‑1 (MQ=255)
                                                                                             ggTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGa                                                     <  1:8094545/150‑1 (MQ=255)
                                                                                              gTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGAt                                                    >  1:7402153/1‑150 (MQ=255)
                                                                                                 gTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGa                                                 >  1:8643807/1‑150 (MQ=255)
                                                                                                   gCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAgg                                               <  1:828892/150‑1 (MQ=255)
                                                                                                      gggTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGgcgc                                            >  2:9121965/1‑150 (MQ=255)
                                                                                                             ggTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTgccgc                                     <  1:6444220/150‑1 (MQ=255)
                                                                                                             ggTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTgccgc                                     >  2:3122161/1‑150 (MQ=255)
                                                                                                              gTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTgccgcc                                    <  2:4901310/150‑1 (MQ=255)
                                                                                                               tGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGccgccg                                   >  1:4879000/1‑150 (MQ=255)
                                                                                                                gCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGa                                  <  1:8613548/150‑1 (MQ=255)
                                                                                                                 cGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGAc                                 <  1:4444707/150‑1 (MQ=255)
                                                                                                                   ggCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTc                               <  2:6503734/150‑1 (MQ=255)
                                                                                                                   ggCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTc                               <  1:344010/150‑1 (MQ=255)
                                                                                                                      gTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGc                            <  1:9743121/150‑1 (MQ=255)
                                                                                                                         gcgcCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTc                         >  2:5412361/1‑150 (MQ=255)
                                                                                                                         gcgcCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTc                         >  2:1685258/1‑150 (MQ=255)
                                                                                                                             cAGTGCCGGTCGAGGTTGAACTCGGCGAGCGAGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCAcc                     <  2:6781944/150‑1 (MQ=255)
                                                                                                                              aGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGTTCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCACCg                    >  1:105526/1‑150 (MQ=255)
                                                                                                                              aGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCACCg                    >  1:7081392/1‑150 (MQ=255)
                                                                                                                                tGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCACCGGc                  <  1:4100762/150‑1 (MQ=255)
                                                                                                                                     gtcgaggttgaACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCACCGGCcgcgc             <  1:9101076/147‑1 (MQ=255)
                                                                                                                                      tcgaggttgaACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCACCGGCcgcgcg            >  1:3502778/3‑150 (MQ=255)
                                                                                                                                      tcgaggttgaACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCACCGGCcgcgcg            <  1:6547809/148‑1 (MQ=255)
                                                                                                                                          ggttgaACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCACCGGCCGCGCGGCCa        >  2:2642782/1‑150 (MQ=255)
                                                                                                                                            ttgaACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCACCGGCCGCGCGGCCACt      >  2:3731042/1‑150 (MQ=255)
                                                                                                                                               aaCTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCACCGGCCGCGCGGCCACTtcg   <  2:8841211/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCTCATCGGGCGCTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTTGAGGTTGAACTTGGCGAGCGTGGCGCGGCTGCCGGCCAGCTTGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCACCGGCCGCGCGGCCACTTCGT  >  NC_008463/3051266‑3051559

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: