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% 1939.6 / NA 499K369E (AAG→GAG) PA14_RS13990SfnB family sulfur acquisition oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (245/253);  minor base G (0/1);  total (245/254)
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 500 total aligned reads displayed.

CTCATCGGGCGCTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTTGAGGTTGAACTTGGCGAGCGTGGCGCGGCTGCCGGCCAGCTTGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCACCGGCCGCGCGGCCACTTCGT  >  NC_008463/3051267‑3051559
                                                                                                                                                  |                                                                                                                                                  
cTCATCGGGCGCTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGc                                                                                                                                                 >  1:7455498/1‑150 (MQ=255)
cTCATCGGGCGCTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGc                                                                                                                                                 >  1:10004462/1‑150 (MQ=255)
     cGGGCGCTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCg                                                                                                                                            <  2:10699413/150‑1 (MQ=255)
     cGGGCGCTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCg                                                                                                                                            >  2:5936410/1‑150 (MQ=255)
      gggCGCTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGt                                                                                                                                           >  1:8915639/1‑150 (MQ=255)
        gcgcTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTgg                                                                                                                                         <  1:1335370/150‑1 (MQ=255)
        gcgcTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGGAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGTACTCGGCGAGCGTgg                                                                                                                                         >  2:2272269/1‑150 (MQ=255)
         cgcTCCTCAGATCCTTGAGTGACGCGCCGGTGGGGCGCCGTTCTGGTAGTAGGCGCCGACCGCGTGGAGGTTCCAGCGTTCCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGGCAGTGCCGGTCGAGGTTGTACTCGGCGAGGGTGgc                                                                                                                                        <  2:8215867/150‑1 (MQ=25)
         cgcTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGgc                                                                                                                                        >  2:7185981/1‑150 (MQ=255)
         cgcTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGgc                                                                                                                                        >  1:3255207/1‑150 (MQ=255)
           ctcctcAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGgcgc                                                                                                                                      <  1:10325464/150‑1 (MQ=255)
            tcctcAGATCCATGAGTGACGCGCCGGCCGGGCGCCGTTCAGGTAGTAGTCGCCGAACGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGcgcg                                                                                                                                     >  1:5858152/1‑150 (MQ=255)
              ctcAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGc                                                                                                                                   <  1:9323628/150‑1 (MQ=255)
               tcAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCt                                                                                                                                  >  2:4770658/1‑150 (MQ=255)
               tcAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCt                                                                                                                                  >  1:9905506/1‑150 (MQ=255)
                cAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTg                                                                                                                                 >  1:7944547/1‑150 (MQ=255)
                 aGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGc                                                                                                                                <  2:5630887/150‑1 (MQ=255)
                 aGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGc                                                                                                                                >  2:9086676/1‑150 (MQ=255)
                  gATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGcc                                                                                                                               >  1:5643386/1‑150 (MQ=255)
                   aTCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCg                                                                                                                              >  2:5320912/1‑150 (MQ=255)
                    tCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCgg                                                                                                                             <  2:6101785/150‑1 (MQ=255)
                     ccATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGc                                                                                                                            >  1:72201/1‑150 (MQ=255)
                      cATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGcc                                                                                                                           >  1:3416652/1‑150 (MQ=255)
                      cATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGcc                                                                                                                           <  2:10541375/150‑1 (MQ=255)
                       aTGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCa                                                                                                                          >  1:2082297/1‑150 (MQ=255)
                        tGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAg                                                                                                                         >  2:7619695/1‑150 (MQ=255)
                         gAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGc                                                                                                                        <  2:6071047/150‑1 (MQ=255)
                           gTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCtc                                                                                                                      >  1:6595911/1‑149 (MQ=255)
                            tGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCtcg                                                                                                                     <  2:8436972/150‑3 (MQ=255)
                            tGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCtcg                                                                                                                     <  1:1665116/150‑3 (MQ=255)
                                 cgcCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAagag                                                                                                                <  2:6939766/150‑1 (MQ=255)
                                 cgcCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAagag                                                                                                                >  1:6072059/1‑150 (MQ=255)
                                  gcCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGc                                                                                                               <  2:4687921/150‑1 (MQ=255)
                                  gcCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGc                                                                                                               >  2:7156311/1‑150 (MQ=255)
                                      gACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTtct                                                                                                           <  2:299104/150‑1 (MQ=255)
                                         gggCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCATCCAGCGAACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGATGCCGGCCAGCTCGAAGAGCTTCTCGc                                                                                                        >  2:1400030/1‑150 (MQ=25)
                                          ggCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCt                                                                                                       <  2:4271744/150‑1 (MQ=255)
                                            cgcCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTgg                                                                                                     <  1:8842411/150‑1 (MQ=255)
                                             gcCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTggc                                                                                                    >  1:7801476/1‑150 (MQ=255)
                                             gcCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGGACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTggc                                                                                                    <  2:1803059/150‑1 (MQ=255)
                                               cGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTggcgg                                                                                                  <  2:2152243/150‑1 (MQ=255)
                                                 ttCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGcc                                                                                                >  2:285052/1‑150 (MQ=255)
                                                 ttCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGcc                                                                                                >  2:10110908/1‑150 (MQ=255)
                                                   cAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAg                                                                                              >  1:2705220/1‑150 (MQ=255)
                                                    aGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAgg                                                                                             >  1:9143530/1‑150 (MQ=255)
                                                    aGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAgg                                                                                             >  1:3988597/1‑150 (MQ=255)
                                                     ggTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGc                                                                                            >  2:1675720/1‑150 (MQ=255)
                                                      gtagtaGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCt                                                                                           >  1:368073/1‑150 (MQ=255)
                                                        agtagtCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGa                                                                                         <  1:9429857/150‑1 (MQ=255)
                                                         gtagtCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGGGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGAt                                                                                        <  2:9315685/150‑1 (MQ=255)
                                                          tagtCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGTtc                                                                                       >  2:94683/1‑150 (MQ=255)
                                                          tagtCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGAtc                                                                                       >  1:5501131/1‑150 (MQ=255)
                                                           agtCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGAtct                                                                                      <  2:4974462/150‑1 (MQ=255)
                                                              cGCCGACCGCGTGGACCTTCCAGCGTAGCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCgg                                                                                   <  2:676054/150‑1 (MQ=255)
                                                                 cGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGc                                                                                >  1:7539682/1‑150 (MQ=255)
                                                                 cGACCGCGTGGACCTTCCAGCGGACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCTGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCGGGCGGCCAGGCTGATCTCGGTGc                                                                                <  2:2495466/150‑1 (MQ=37)
                                                                   aCCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTc                                                                              >  2:4036073/1‑150 (MQ=255)
                                                                    ccGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCa                                                                             >  1:7006200/1‑150 (MQ=255)
                                                                     cgcgTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAg                                                                            <  1:7658438/150‑1 (MQ=255)
                                                                      gcgTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGc                                                                           <  2:8706283/150‑1 (MQ=255)
                                                                       cgTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCGTCTCGCTGGCGGCCAGGCTGATCTCGGTGCGCCGCa                                                                          >  2:3514738/1‑150 (MQ=37)
                                                                          ggACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCt                                                                       <  1:9254115/150‑1 (MQ=255)
                                                                          ggACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCt                                                                       <  1:8546217/150‑1 (MQ=255)
                                                                           gACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCtt                                                                      <  2:6696531/150‑1 (MQ=255)
                                                                           gACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCtt                                                                      <  2:7758687/150‑1 (MQ=255)
                                                                            aCCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTg                                                                     <  1:9509727/150‑1 (MQ=255)
                                                                             ccTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTgg                                                                    >  1:10651134/1‑150 (MQ=255)
                                                                              cTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGc                                                                   <  1:4362099/150‑1 (MQ=255)
                                                                               gtCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGcc                                                                  <  1:515142/149‑1 (MQ=255)
                                                                                 ccAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTc                                                                >  1:6800750/1‑150 (MQ=255)
                                                                                       tACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCAc                                                          <  2:1282040/150‑1 (MQ=255)
                                                                                        aCCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACg                                                         <  2:3150446/150‑1 (MQ=255)
                                                                                         ccGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACgg                                                        <  1:10391745/150‑1 (MQ=255)
                                                                                          cGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGc                                                       >  1:4509840/1‑150 (MQ=255)
                                                                                             gTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGAt                                                    <  2:11176261/150‑1 (MQ=255)
                                                                                             gTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGAt                                                    <  1:8136149/150‑1 (MQ=255)
                                                                                             gTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGAt                                                    <  2:10940953/150‑1 (MQ=255)
                                                                                               cGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCg                                                  <  1:2132663/150‑1 (MQ=255)
                                                                                                gTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGa                                                 >  1:461864/1‑150 (MQ=255)
                                                                                                  gCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAgg                                               <  1:10931418/150‑1 (MQ=255)
                                                                                                   cAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGgc                                              <  1:5775088/150‑1 (MQ=255)
                                                                                                    aGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCGAGGCTGATCTCGGTGGTCAGCACCGTGGCGTCGGCCACGGCGATCGAGgcg                                             >  2:7893161/1‑150 (MQ=37)
                                                                                                    aGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGgcg                                             <  1:1896518/150‑1 (MQ=255)
                                                                                                     gggTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGgcgc                                            >  2:1557255/1‑150 (MQ=255)
                                                                                                      ggTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGcgcg                                           >  2:10303395/1‑150 (MQ=255)
                                                                                                                cGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGAc                                 >  2:6426827/1‑150 (MQ=255)
                                                                                                                    cGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGACGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCgg                             <  2:2602166/150‑1 (MQ=255)
                                                                                                                     gTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGc                            <  2:9615848/150‑1 (MQ=255)
                                                                                                                          gcCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCa                       >  1:4770666/1‑150 (MQ=255)
                                                                                                                           ccAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCAc                      <  2:3713849/150‑1 (MQ=255)
                                                                                                                           ccAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCAc                      >  2:8374995/1‑150 (MQ=255)
                                                                                                                           ccAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCAc                      <  1:2954930/150‑1 (MQ=255)
                                                                                                                             aGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCACCg                    <  1:2511074/150‑1 (MQ=255)
                                                                                                                               tGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCACCGGc                  <  1:9966799/150‑1 (MQ=255)
                                                                                                                               tGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCACCGGc                  <  1:9961119/150‑1 (MQ=255)
                                                                                                                                    gtcgaggttgaACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCACCGGCcgcgc             <  1:304591/147‑1 (MQ=255)
                                                                                                                                        aggttgaACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACTGCGAGCGAGGCGCGTGCCGCCGACTCGGCGTCCACCGGCCGCGCGGcc         >  1:9033706/1‑150 (MQ=255)
                                                                                                                                        aggttgaACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCACCGGCCGCGCGGcc         >  1:4256971/1‑150 (MQ=255)
                                                                                                                                             gaACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCACCGGCCGCGCGGCCACTtc    <  1:1165588/150‑1 (MQ=255)
                                                                                                                                               aCTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCACCGGCCGCGCGGCCACTtcgt  <  2:11085273/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CTCATCGGGCGCTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTTGAGGTTGAACTTGGCGAGCGTGGCGCGGCTGCCGGCCAGCTTGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCACCGGCCGCGCGGCCACTTCGT  >  NC_008463/3051267‑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: