Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 3,051,401 T→C N373D (AAC→GAC)  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,4010TC99.8% 1801.7 / NA 464N373D (AAC→GAC) PA14_RS13990SfnB family sulfur acquisition oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (266/197);  minor base G (0/1);  total (266/198)
Fisher's exact test for biased strand distribution p-value = 4.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 491 total aligned reads displayed.

AGCGTGTGCGAGGCTCATCGGGCGCTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTTGAGGTTGAACTTGGCGAGCGTGGCGCGGCTGCCGGCCAGCTTGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCACCGGCCGCG  >  NC_008463/3051254‑3051547
                                                                                                                                                   |                                                                                                                                                  
 gCGTGTGCGAGGCTCATCGGGCGCTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCggtcgag                                                                                                                                                 >  1:9372050/1‑150 (MQ=255)
 gCGTGTGCGAGGCTCATCGGGCGCTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCggtcgag                                                                                                                                                 <  2:3111386/150‑1 (MQ=255)
 gCGTGTGCGAGGCTCATCGGGCGCTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCggtcgag                                                                                                                                                 >  1:5840887/1‑150 (MQ=255)
 gCGTGTGCGAGGCTCATCGGGCGCTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGGCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCggtcgag                                                                                                                                                 >  1:9566677/1‑150 (MQ=255)
    tgtgCGAGGCTCATCGGGCGCTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCggtcgaggtt                                                                                                                                              >  2:7371547/1‑150 (MQ=255)
    tgtgCGAGGCTCATCGGGCGCTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCggtcgaggtt                                                                                                                                              <  2:3165117/150‑1 (MQ=255)
       gCGAGGCTCATCGGGCGCTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGaa                                                                                                                                           >  2:3614575/1‑150 (MQ=255)
         gAGGCTCATCGGGCGCTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACt                                                                                                                                         <  2:5379348/150‑1 (MQ=255)
           ggCTCATCGGGCGCTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACtcg                                                                                                                                       <  2:4838960/150‑3 (MQ=255)
            gCTCATCGGGCGCTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCgg                                                                                                                                      <  1:2083954/150‑1 (MQ=255)
             cTCATCGGGCGCTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGc                                                                                                                                     <  1:7819847/150‑1 (MQ=255)
             cTCATCGGGCGCTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGc                                                                                                                                     <  1:8975352/150‑1 (MQ=255)
                   gggCGCTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGt                                                                                                                               >  2:2304599/1‑150 (MQ=255)
                    ggCGCTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTg                                                                                                                              >  2:2887168/1‑150 (MQ=255)
                    ggCGCTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTg                                                                                                                              <  2:1309671/150‑1 (MQ=255)
                     gcgcTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTgg                                                                                                                             >  1:1528440/1‑150 (MQ=255)
                      cgcTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGgc                                                                                                                            <  1:6757905/150‑1 (MQ=255)
                       gcTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGgcg                                                                                                                           >  2:5192489/1‑150 (MQ=255)
                         tcctcAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGcgcg                                                                                                                         >  1:2299210/1‑150 (MQ=255)
                         tcctcAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGcgcg                                                                                                                         >  2:4757262/1‑150 (MQ=255)
                           ctcAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGc                                                                                                                       >  1:5438670/1‑150 (MQ=255)
                              aGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGc                                                                                                                    >  2:7772579/1‑150 (MQ=255)
                               gATCCTTGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGcc                                                                                                                   >  2:1540602/1‑150 (MQ=255)
                                aTCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCg                                                                                                                  <  2:11058518/150‑1 (MQ=255)
                                 tCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCgg                                                                                                                 <  1:1701699/150‑1 (MQ=255)
                                   cATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGcc                                                                                                               <  1:7173150/150‑1 (MQ=255)
                                      gAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGc                                                                                                            >  2:1943246/1‑150 (MQ=255)
                                       aGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCt                                                                                                           >  2:6316336/1‑150 (MQ=255)
                                       aGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCt                                                                                                           <  2:8573158/150‑1 (MQ=255)
                                        gTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCtc                                                                                                          >  2:1162790/1‑149 (MQ=255)
                                          gACGCGCCGGACGGGCGCCGTTCAGGTATTAGTCGCAGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGa                                                                                                        >  1:124555/1‑150 (MQ=25)
                                          gACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGa                                                                                                        <  1:8682595/150‑1 (MQ=255)
                                             gcgcCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAaga                                                                                                     >  2:1631920/1‑150 (MQ=255)
                                              cgcCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAagag                                                                                                    >  1:10933942/1‑150 (MQ=255)
                                               gcCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGc                                                                                                   <  1:10569156/150‑1 (MQ=255)
                                                ccGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGTGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCt                                                                                                  >  1:4892794/1‑150 (MQ=255)
                                                ccGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCt                                                                                                  >  2:9851667/1‑150 (MQ=255)
                                                  ggACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTtc                                                                                                >  2:3419897/1‑150 (MQ=255)
                                                  ggACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTtc                                                                                                >  2:7917535/1‑150 (MQ=255)
                                                   gACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTtct                                                                                               >  1:7549037/1‑150 (MQ=255)
                                                    aCGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTtctc                                                                                              <  2:4488675/150‑1 (MQ=255)
                                                      gggCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGc                                                                                            <  1:2164424/150‑1 (MQ=255)
                                                       ggCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCt                                                                                           <  1:653450/150‑1 (MQ=255)
                                                       ggCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCt                                                                                           >  1:10008247/1‑150 (MQ=255)
                                                           ccGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTggcg                                                                                       >  2:2584716/1‑150 (MQ=255)
                                                           ccGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTggcg                                                                                       >  2:8754931/1‑150 (MQ=255)
                                                            cGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTggcgg                                                                                      >  1:6105719/1‑150 (MQ=255)
                                                             gTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTggcggc                                                                                     >  1:10826733/1‑150 (MQ=255)
                                                               tCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCa                                                                                   >  1:3721964/1‑150 (MQ=255)
                                                               tCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCa                                                                                   <  2:5595509/150‑1 (MQ=255)
                                                                     agtagtCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGa                                                                             <  2:768415/150‑1 (MQ=255)
                                                                      gtagtCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGag                                                                            >  2:10658036/1‑149 (MQ=255)
                                                                      gtagtCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGAt                                                                            >  1:6409851/1‑150 (MQ=255)
                                                                        agtCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGAtct                                                                          <  1:249177/150‑1 (MQ=255)
                                                                           cGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCgg                                                                       <  2:10486234/150‑1 (MQ=255)
                                                                             ccGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTg                                                                     >  2:6447136/1‑150 (MQ=255)
                                                                                aCCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGAGCTCGGTGcgc                                                                  >  1:4645729/1‑148 (MQ=255)
                                                                                  cgcgTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAg                                                                <  2:969179/150‑1 (MQ=255)
                                                                                    cgTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCa                                                              >  1:9921186/1‑150 (MQ=255)
                                                                                      tGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCAcc                                                            >  2:1467736/1‑150 (MQ=255)
                                                                                          ccTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTgg                                                        >  2:3752197/1‑150 (MQ=255)
                                                                                             tCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCt                                                     <  1:103695/150‑1 (MQ=255)
                                                                                             tCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCt                                                     >  1:10491220/1‑150 (MQ=255)
                                                                                              ccAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTc                                                    >  2:5966304/1‑150 (MQ=255)
                                                                                               cAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCg                                                   >  1:4115707/1‑150 (MQ=255)
                                                                                                 gCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGc                                                 <  1:1125508/150‑1 (MQ=255)
                                                                                                   gTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCa                                               >  1:5240207/1‑150 (MQ=255)
                                                                                                    tACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCAc                                              >  2:4273493/1‑150 (MQ=255)
                                                                                                    tACCGGGTCGTGCAGGGTGGGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCAc                                              <  2:7477048/150‑1 (MQ=255)
                                                                                                       cGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGc                                           >  1:5751902/1‑150 (MQ=255)
                                                                                                       cGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGc                                           <  2:9317399/150‑1 (MQ=255)
                                                                                                        gggTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCg                                          >  1:5992925/1‑150 (MQ=255)
                                                                                                        gggTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCg                                          >  2:3246857/1‑150 (MQ=255)
                                                                                                               gCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAgg                                   >  2:6779099/1‑150 (MQ=255)
                                                                                                                  gggTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGGTGATCTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGTGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCCCCTTGGCCTCGGCCACGGCGATCGAGgcgc                                >  1:876252/1‑150 (MQ=37)
                                                                                                                    gtgtGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGt                              <  2:9518161/150‑1 (MQ=255)
                                                                                                                      gtgGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTgc                            >  2:2652472/1‑150 (MQ=255)
                                                                                                                         ggTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCTAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTgccgc                         >  1:4492470/1‑150 (MQ=255)
                                                                                                                         ggTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTgccgc                         <  2:9028301/150‑1 (MQ=255)
                                                                                                                         ggTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTgccgc                         >  1:5609583/1‑150 (MQ=255)
                                                                                                                         ggTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTgccgc                         >  1:3859322/1‑150 (MQ=255)
                                                                                                                          gTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTgccgcc                        >  1:8173253/1‑150 (MQ=255)
                                                                                                                              gggCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACt                    >  1:4360272/1‑150 (MQ=255)
                                                                                                                                gCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCg                  <  1:8042913/150‑1 (MQ=255)
                                                                                                                                gCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCg                  <  1:9050578/150‑1 (MQ=255)
                                                                                                                                gCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCg                  <  2:6569954/150‑1 (MQ=255)
                                                                                                                                   ttGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCg               <  2:4593984/150‑1 (MQ=255)
                                                                                                                                      cgcCAGTGCCGGTCGAGGGTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGAACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGTCTCGGCGTcc            >  2:7212715/1‑150 (MQ=39)
                                                                                                                                       gcCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCa           <  2:5700784/150‑1 (MQ=255)
                                                                                                                                        ccAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCAc          >  1:2708535/1‑150 (MQ=255)
                                                                                                                                        ccAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCAc          >  1:85416/1‑150 (MQ=255)
                                                                                                                                        ccAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCAc          >  2:8832604/1‑150 (MQ=255)
                                                                                                                                        ccAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCAc          <  2:8994656/150‑1 (MQ=255)
                                                                                                                                         cAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCAcc         >  1:3079212/1‑150 (MQ=255)
                                                                                                                                           gTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCACCgg       <  1:3374643/150‑1 (MQ=255)
                                                                                                                                           gTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCACCgg       >  1:8919680/1‑150 (MQ=255)
                                                                                                                                             gCCGGTCGAGGGTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCACCGGcc     <  1:2890289/150‑1 (MQ=255)
                                                                                                                                              ccGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCACCGGCcg    >  1:398010/1‑150 (MQ=255)
                                                                                                                                              ccGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCACCGGCcg    >  2:4443435/1‑150 (MQ=255)
                                                                                                                                                ggtcgaggttgaACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCACCGGCcgcg  <  2:8193797/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AGCGTGTGCGAGGCTCATCGGGCGCTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTTGAGGTTGAACTTGGCGAGCGTGGCGCGGCTGCCGGCCAGCTTGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCACCGGCCGCG  >  NC_008463/3051254‑3051547

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: