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% 2122.2 / NA 544K369E (AAG→GAG) PA14_RS13990SfnB family sulfur acquisition oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (281/262);  minor base A (0/1);  total (281/263)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 544 total aligned reads displayed.

GCTCATCGGGCGCTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTTGAGGTTGAACTTGGCGAGCGTGGCGCGGCTGCCGGCCAGCTTGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCACCGGCCGCGCGGCCACTTCGT  >  NC_008463/3051266‑3051559
                                                                                                                                                   |                                                                                                                                                  
   cATCGGGCGCTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGa                                                                                                                                               >  1:4315985/1‑150 (MQ=255)
    aTCGGGCGCTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAg                                                                                                                                              >  1:6180245/1‑150 (MQ=255)
    aTCGGGCGCTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAg                                                                                                                                              >  1:8631692/1‑150 (MQ=255)
    aTCGGGCGCTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAg                                                                                                                                              >  1:7147911/1‑150 (MQ=255)
     tCGGGCGCTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGc                                                                                                                                             <  2:698075/150‑1 (MQ=255)
       gggCGCTCCTCAGCTCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGt                                                                                                                                           <  2:9792607/150‑1 (MQ=255)
        ggCGCTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTg                                                                                                                                          >  1:1730937/1‑150 (MQ=255)
        ggCGCTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTg                                                                                                                                          >  1:3710935/1‑150 (MQ=255)
        ggCGCTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTg                                                                                                                                          <  2:309345/150‑1 (MQ=255)
          cgcTCCTCAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGgc                                                                                                                                        <  2:8839569/150‑1 (MQ=255)
            ctcctcAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGgcgc                                                                                                                                      <  2:10261159/150‑1 (MQ=255)
              cctcAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCgg                                                                                                                                    <  2:6513615/150‑1 (MQ=255)
                tcAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCt                                                                                                                                  >  2:1909034/1‑150 (MQ=255)
                tcAGATCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCt                                                                                                                                  <  2:10805108/150‑1 (MQ=255)
                    aTCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCg                                                                                                                              >  1:7617783/1‑150 (MQ=255)
                     tCCATGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCgg                                                                                                                             <  1:6874151/150‑1 (MQ=255)
                         tGAGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAg                                                                                                                         <  1:9807415/150‑1 (MQ=255)
                           aGTGACGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCt                                                                                                                       <  2:4012519/150‑1 (MQ=255)
                               aCGCGCCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGaa                                                                                                                   >  2:4452396/1‑150 (MQ=255)
                               aCGCGCCGGACGGGCGCCGTTCAGGTAGAAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGaa                                                                                                                   >  1:6069640/1‑150 (MQ=255)
                                   gcCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGc                                                                                                               >  1:4460911/1‑150 (MQ=255)
                                   gcCGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGc                                                                                                               >  2:5398673/1‑150 (MQ=255)
                                    ccGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCt                                                                                                              <  2:1189740/150‑1 (MQ=255)
                                    ccGGACGGGCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCt                                                                                                              <  1:7026276/150‑1 (MQ=255)
                                      tgACGGGCGCCGTTCTGGTAGTAGTGGCAGCCCGCGTGGATCTTCCAGCGTACGGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGCTGAACTCGGCGTTCGTGGCGGGGCTGCCGTCCAGCTCTAAGAGCTtc                                                                                                            <  2:8465292/149‑1 (MQ=25)
                                        aCGGGCGCCGTACAGGTAGTAGTCGCCGACCGCGTGCACCTTCCAGCGTACCGGGTCGAGCAGGGTGAGGGTGCGGGCGATGCGCCAGAGCCGGACGAGGTAGAACTCGGCGAGCGTGGCGCGGCTCCCGGCCAGCTCGAAGAGCTactc                                                                                                          >  2:83118/1‑150 (MQ=21)
                                          gggCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGc                                                                                                        <  2:6913861/150‑1 (MQ=255)
                                           ggCGCCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCt                                                                                                       >  2:914198/1‑150 (MQ=255)
                                             cgcCGTTCAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTgg                                                                                                     <  2:7495020/150‑1 (MQ=255)
                                                  ttCAGGGAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGcc                                                                                                <  1:10770502/150‑1 (MQ=255)
                                                   tCAGGTAGTAGTCGCCGACCGCGTTGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCa                                                                                               <  2:7993655/150‑1 (MQ=255)
                                                    cAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAg                                                                                              <  1:2964000/150‑1 (MQ=255)
                                                    cAGGTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAg                                                                                              >  1:2393393/1‑150 (MQ=255)
                                                      ggTAGTAGTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGc                                                                                            >  1:3153264/1‑150 (MQ=255)
                                                        tagtagTCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTg                                                                                          <  2:9430349/150‑1 (MQ=255)
                                                             gtCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCCGGGTGTGGGTGCGGGCGTTGCGCCAGTGCAGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGAtctc                                                                                     >  1:2646155/1‑150 (MQ=38)
                                                             gtCGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGAtctc                                                                                     >  2:3616115/1‑150 (MQ=255)
                                                               cGCCGACCGCGTGGACCTTCCAGCGTTCCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCgg                                                                                   <  2:5280936/150‑1 (MQ=255)
                                                               cGCCGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCgg                                                                                   <  1:4589747/150‑1 (MQ=255)
                                                                  cGACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGc                                                                                >  2:2812501/1‑150 (MQ=255)
                                                                   gACCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCt                                                                               >  1:7360883/1‑150 (MQ=255)
                                                                    aCCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTc                                                                              >  2:2553399/1‑150 (MQ=255)
                                                                    aCCGCGTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCATCTCGCTGGCGGCCAGGCTGATCTCGGTGCTc                                                                              >  1:10006214/1‑150 (MQ=255)
                                                                      cgcgTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAg                                                                            >  1:6546586/1‑150 (MQ=255)
                                                                      cgcgTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAg                                                                            >  2:5002904/1‑150 (MQ=255)
                                                                       gcgTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGc                                                                           <  2:5588265/150‑1 (MQ=255)
                                                                       gcgTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGc                                                                           <  2:6386700/150‑1 (MQ=255)
                                                                       gcgTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGc                                                                           >  1:999608/1‑150 (MQ=255)
                                                                        cgTGGACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCa                                                                          >  1:8829845/1‑150 (MQ=255)
                                                                           ggACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCt                                                                       >  1:5420934/1‑150 (MQ=255)
                                                                           ggACCTTCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCt                                                                       <  1:8359771/150‑1 (MQ=255)
                                                                                 tCCAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGGTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCt                                                                 <  1:8869234/150‑1 (MQ=255)
                                                                                   cAGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCg                                                               <  2:1465843/150‑1 (MQ=255)
                                                                                    aGCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGGGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCGGCACCTTGGCCTCgg                                                              >  2:3798601/1‑150 (MQ=255)
                                                                                     gCGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGc                                                             >  2:10399472/1‑150 (MQ=255)
                                                                                      cGTACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGcc                                                            <  2:6789154/150‑1 (MQ=255)
                                                                                        tACCGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCAc                                                          >  1:9405472/1‑150 (MQ=255)
                                                                                          ccGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACgg                                                        <  1:8011750/150‑1 (MQ=255)
                                                                                           cGGGTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGc                                                       <  2:10101753/150‑1 (MQ=255)
                                                                                            gggTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCg                                                      >  2:3214665/1‑150 (MQ=255)
                                                                                              gTCGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGAt                                                    <  1:2064473/150‑1 (MQ=255)
                                                                                                cGTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCg                                                  >  2:6445644/1‑150 (MQ=255)
                                                                                                 gTGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGa                                                 <  2:2686444/150‑1 (MQ=255)
                                                                                                  tGCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAg                                                >  2:10520868/1‑150 (MQ=255)
                                                                                                   gCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAgg                                               >  2:7571887/1‑150 (MQ=255)
                                                                                                   gCAGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAgg                                               >  1:10221456/1‑150 (MQ=255)
                                                                                                     aGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGgcg                                             >  2:2981169/1‑150 (MQ=255)
                                                                                                     aGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGgcg                                             >  1:9562815/1‑150 (MQ=255)
                                                                                                     aGGGTGTGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGgcg                                             <  2:411431/150‑1 (MQ=255)
                                                                                                        gtgtGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGt                                          >  1:10600243/1‑150 (MQ=255)
                                                                                                        gtgtGGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGt                                          <  1:6362188/150‑1 (MQ=255)
                                                                                                         tgtgGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTg                                         >  1:10473220/1‑150 (MQ=255)
                                                                                                           tgGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTgcc                                       >  1:4963182/1‑150 (MQ=255)
                                                                                                           tgGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTgcc                                       <  1:1180544/150‑1 (MQ=255)
                                                                                                           tgGGTGCGGGCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTgcc                                       <  1:7849246/150‑1 (MQ=255)
                                                                                                                  gggCGTTGCGCCAGTGCCGGTCGCGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGCTCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACt                                <  2:2120326/150‑1 (MQ=38)
                                                                                                                   ggCGTTGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTc                               >  1:10098219/1‑150 (MQ=255)
                                                                                                                       ttGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCg                           <  1:5978484/150‑1 (MQ=255)
                                                                                                                       ttGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCg                           >  2:11010342/1‑150 (MQ=255)
                                                                                                                       ttGCGCCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCg                           <  1:5649602/150‑1 (MQ=255)
                                                                                                                           gcCAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCa                       >  1:10954072/1‑150 (MQ=255)
                                                                                                                            ccAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCAc                      >  1:8389816/1‑150 (MQ=255)
                                                                                                                             cAGTGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCAcc                     <  2:5064001/150‑1 (MQ=255)
                                                                                                                                tGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCACCGGc                  <  1:9180582/150‑1 (MQ=255)
                                                                                                                                tGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCACCGGc                  >  1:9672787/1‑150 (MQ=255)
                                                                                                                                tGCCGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGAGCGAGGCGCGTGCCGCCGACTCGGCGTCCACCGGc                  >  2:6701408/1‑150 (MQ=255)
                                                                                                                                  ccGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCACCGGCcg                >  2:2409023/1‑150 (MQ=255)
                                                                                                                                   cGGTCGAGGTTGAACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCACCGGCcgc               >  1:5798153/1‑150 (MQ=255)
                                                                                                                                     gtcgaggttgaACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGTTCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCACCGGCcgcgc             >  2:8743002/4‑150 (MQ=255)
                                                                                                                                      tcgaggttgaACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCACCGGCcgcgcg            >  1:2291759/3‑150 (MQ=255)
                                                                                                                                       cgaggttgaACTCGGCGAGCGTGGCGCGGCTGCCGGACATCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCACCGGCCGCGCgg           <  2:5967098/149‑1 (MQ=255)
                                                                                                                                        gaggttgaACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCACCGGCCGCGCGGc          <  1:3241985/150‑1 (MQ=255)
                                                                                                                                         aggttgaACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCACCGGCCGCGCGGcc         >  2:4658752/1‑150 (MQ=255)
                                                                                                                                          ggttgaACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCACCGGCCGCGCGGCCa        >  2:7790846/1‑150 (MQ=255)
                                                                                                                                          ggttgaACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCACCGGCCGCGCGGCCa        <  2:1652616/150‑1 (MQ=255)
                                                                                                                                           gttgaACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCACCGGCCGCGCGGCCAc       <  1:1844651/150‑1 (MQ=255)
                                                                                                                                           gttgaACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCACCGGCCGCGCGGCCAc       >  1:9897150/1‑150 (MQ=255)
                                                                                                                                            ttgaACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCACCGGCCGCGCGGCCACt      <  2:6212161/150‑1 (MQ=255)
                                                                                                                                             tgaACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCACCGGCCGCGCGGCCACtt     <  1:291560/150‑1 (MQ=255)
                                                                                                                                             tgaACTCGGCGAGCGTGGCGCGGCTGCCGGCCAGCTCGAAGAGCTTCTCGCTGGCGGCCAGGCTGATCTCGGTGCTCAGCACCTTGGCCTCGGCCACGGCGATCGAGGCGCGTGCCGCCGACTCGGCGTCCACCGGCCGCGCGGCCACtt     <  1:5887694/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: