Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 5,198,405 T→C intergenic (+31/+37) PA14_RS23800 → / ← mdpA ABC transporter substrate‑binding protein/metallopeptidase MdpA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084635,198,4050TC99.8% 1498.8 / NA 393intergenic (+31/+37)PA14_RS23800/mdpAABC transporter substrate‑binding protein/metallopeptidase MdpA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (168/224);  minor base G (0/1);  total (168/225)
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 394 total aligned reads displayed.

GATGCCCTGGTTGCCGCTGGCGCATCCGTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCTTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGATCCTCCAGGCGAATGCCGAACTCGCC  >  NC_008463/5198257‑5198552
                                                                                                                                                    |                                                                                                                                                   
 cagCCCTGGTTGCCGCTGGCGCATCAGTTGATGTTCGACGCGCACCGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGAACGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTg                                                                                                                                                  <  2:4225189/149‑1 (MQ=255)
 aTGCCCTGGTTGCCGCTGGCGCATCCGTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTg                                                                                                                                                  >  2:4908813/1‑151 (MQ=255)
  tGCCCTGGTTGCCGCTGGCGCATCCGTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGc                                                                                                                                                 >  2:3597616/1‑151 (MQ=255)
     ccTGGTTGCCGCTGGCGCATCCGTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCtt                                                                                                                                              >  2:3372712/1‑151 (MQ=255)
       tGGTTGCCGCTGGCGCATCCGTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTcc                                                                                                                                            >  1:9706832/1‑151 (MQ=255)
            gCCGCTGGCGCATCCGTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGt                                                                                                                                       <  1:784871/151‑1 (MQ=255)
              cGCTGGCGCATCCGTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTgc                                                                                                                                     <  2:2323173/151‑1 (MQ=255)
               gCTGGCGCATCCGTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTgcg                                                                                                                                    >  1:3441029/1‑151 (MQ=255)
                cTGGCGCATCCGTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCgcgggcgg                                                                                                                                   >  1:2437869/1‑151 (MQ=255)
                cTGGCGCATCCGTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTgcgg                                                                                                                                   <  2:10359094/151‑1 (MQ=255)
                   gcgcATCCGTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGt                                                                                                                                <  1:5389490/151‑1 (MQ=255)
                   gcgcATCCGTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGt                                                                                                                                >  2:3504877/1‑151 (MQ=255)
                     gcATCCGTTGCTGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTa                                                                                                                              >  1:8095646/1‑151 (MQ=255)
                     gcATCCGTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCCCGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTa                                                                                                                              <  2:5634693/151‑1 (MQ=255)
                       ccccGCTGATGCTCGACTCGCCACGCACGCACGTCAGCGGCTCCCGCAAGAGCCCGAGGTCGGCGCGGGACATTACGCGGGAGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTGGCCTTCCCGCGTGCGGCGTTATc                                                                                                                            <  2:6656222/149‑1 (MQ=255)
                        tCCGTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCg                                                                                                                           <  1:1632901/151‑1 (MQ=255)
                         gcGTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGGCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGTAGCCCGGGCCCTTGCCTTCCCGCGTGCGGCGtgttcgg                                                                                                                          >  1:1099214/2‑145 (MQ=255)
                         ccGTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGt                                                                                                                          >  2:4015602/1‑151 (MQ=255)
                         ccGTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGt                                                                                                                          >  2:8574483/1‑151 (MQ=255)
                           gTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTc                                                                                                                        >  2:365493/1‑151 (MQ=255)
                           gTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTc                                                                                                                        >  2:2799116/1‑151 (MQ=255)
                           gTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCATGCCTTCCCGCGTGCGGCGTTATCGGTc                                                                                                                        >  1:3201340/1‑151 (MQ=255)
                            ttGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTcg                                                                                                                       <  1:6915724/151‑1 (MQ=255)
                             tGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTcgg                                                                                                                      <  2:10580866/151‑1 (MQ=255)
                               aTGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATAGTTcgggc                                                                                                                    >  1:2785395/1‑149 (MQ=255)
                                tGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGTAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGGTcggccg                                                                                                                   >  1:2634423/1‑151 (MQ=255)
                                tGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGGCCTGGCCCTTGCCTTCCCGCGGGCGGCGTTATCGTTcggccg                                                                                                                   >  1:9478555/1‑151 (MQ=255)
                                 gTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTcggccgg                                                                                                                  >  2:7238412/1‑151 (MQ=255)
                                  ttCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTcggccggc                                                                                                                 <  1:6172102/151‑1 (MQ=255)
                                    cGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCGTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTc                                                                                                               >  1:7660110/1‑151 (MQ=255)
                                     gACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCa                                                                                                              <  1:6695928/151‑1 (MQ=255)
                                     gACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCa                                                                                                              <  2:8139514/151‑1 (MQ=255)
                                        gcgcAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCGTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAgcc                                                                                                           <  2:10687202/151‑1 (MQ=255)
                                        gcgcAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCGTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAgcc                                                                                                           <  2:1583425/151‑1 (MQ=255)
                                           caacgcaaGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAgccgcc                                                                                                        <  2:7997593/151‑1 (MQ=255)
                                           caacgcaaGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCGTGCCGTCCCGCGGGCGGCGGTATCGTTCGGCCGGCTCAgccgcc                                                                                                        >  1:2850950/1‑151 (MQ=255)
                                            aacgcaaGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGccgccg                                                                                                       >  2:8441532/1‑151 (MQ=255)
                                              cgcaaGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGACGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAg                                                                                                     <  2:10445396/151‑1 (MQ=255)
                                              cgcaaGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAg                                                                                                     <  1:5173794/151‑1 (MQ=255)
                                              cgcaaGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCGTGCCGTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAg                                                                                                     >  1:3322785/1‑151 (MQ=255)
                                               gcaaGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCGTCCCGCGTGCGGCGTTATCGTTAGGCCGGCTCAGCCGCCGAgg                                                                                                    >  1:1416545/1‑151 (MQ=255)
                                                caaGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGc                                                                                                   <  2:10134555/151‑1 (MQ=255)
                                                 aaGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGcc                                                                                                  <  2:9887441/151‑1 (MQ=255)
                                                 aaGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGcc                                                                                                  <  1:9064387/151‑1 (MQ=255)
                                                   gAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGGTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGa                                                                                                >  2:7465318/1‑151 (MQ=255)
                                                   gAACGTCAGCGGCTACCGCATGAGCCCGATGACGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCAGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGa                                                                                                <  1:8440929/151‑1 (MQ=255)
                                                     aCGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAg                                                                                              >  2:6815474/1‑151 (MQ=255)
                                                     aCGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAg                                                                                              <  1:10658094/151‑1 (MQ=255)
                                                         cAGCGGCTACCGCATGATCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCGCGCGTGCGGCGGTATCGTTAGGCCGGCTCTGCCGCCGAGGCCGTCGGGAt                                                                                          >  1:5936161/1‑151 (MQ=255)
                                                          aGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGAtc                                                                                         <  2:3033406/151‑1 (MQ=255)
                                                          aGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGAtc                                                                                         <  1:9864631/151‑1 (MQ=255)
                                                            cGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGAtcgt                                                                                       >  1:7289290/1‑151 (MQ=255)
                                                            cGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCGTGCCTTCCCGCGTGCGGGGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGAtcgt                                                                                       >  1:1742771/1‑151 (MQ=255)
                                                             ggCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGAtcgtc                                                                                      <  2:5560869/151‑1 (MQ=255)
                                                             ggCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGAtcgtc                                                                                      >  2:9324424/1‑151 (MQ=255)
                                                              gCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGAtcgtcg                                                                                     <  2:10222389/151‑1 (MQ=255)
                                                                tACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGAt                                                                                   >  1:4457976/1‑151 (MQ=255)
                                                                tACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCGTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGAt                                                                                   >  1:8881580/1‑151 (MQ=255)
                                                                 aCCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATg                                                                                  <  1:4032581/151‑1 (MQ=255)
                                                                   cGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATgga                                                                                <  1:7182117/151‑1 (MQ=255)
                                                                    gCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATggac                                                                               <  2:7913402/151‑1 (MQ=255)
                                                                          gCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTc                                                                         <  1:1306577/151‑1 (MQ=255)
                                                                           cccGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCg                                                                        <  1:4961354/151‑1 (MQ=255)
                                                                           cccGATGTCGGCGCGGGACTTTTCGCGGGCGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCg                                                                        >  2:3844451/1‑151 (MQ=255)
                                                                                tGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTgg                                                                   <  1:10540952/151‑1 (MQ=255)
                                                                                  tCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTTCCTTCCCGCGGGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGTCGGACGCGGTTgggg                                                                 >  1:10093836/1‑150 (MQ=255)
                                                                                       gcgGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGGGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGCTCGTCGATGGACGGACTCGGGTGGGTGAAcc                                                            >  2:1406776/1‑151 (MQ=255)
                                                                                        cgtgACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGATGCCCTGGCCGTTGCCTTCCCGCGTGCGGCGTTATCGGTCGGCCGGATCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGTGTGAACCa                                                           <  2:8790572/148‑1 (MQ=255)
                                                                                        cgGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGGTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCa                                                           >  2:837410/1‑151 (MQ=255)
                                                                                         gggACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAg                                                          <  1:3098583/151‑1 (MQ=255)
                                                                                          ggACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGc                                                         <  1:7013026/151‑1 (MQ=255)
                                                                                               tttCGCGGGTGAAGCTGGACTGATCGCTGCGACGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGt                                                    <  2:262162/151‑1 (MQ=255)
                                                                                                   gcgGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGtcctt                                                <  2:4752875/151‑1 (MQ=255)
                                                                                                     gggTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGtccttac                                              >  2:212081/1‑149 (MQ=255)
                                                                                                        tGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTtcctcc                                           <  2:11032187/151‑1 (MQ=255)
                                                                                                        tGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTtcctcc                                           <  1:2280639/151‑1 (MQ=255)
                                                                                                           aGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTc                                        <  1:10285756/151‑1 (MQ=255)
                                                                                                                gACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTa                                   <  1:4214221/151‑1 (MQ=255)
                                                                                                                  cTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCATCCCGCGGGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGTTCGTCGATGGACGGTCTCGGTTGGGTGACCCCGCGTGGTGCGTCCTCCGTCATGTAgt                                 >  2:2905665/1‑150 (MQ=255)
                                                                                                                   tGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGGCCTTCCTCCGTCATGTAGaa                                >  2:3141405/1‑151 (MQ=255)
                                                                                                                     aTCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGt                              <  1:5827720/151‑1 (MQ=255)
                                                                                                                     aTCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGt                              <  1:5554950/151‑1 (MQ=255)
                                                                                                                     aTCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGCGGTCCTTCCTCCGTCATGTAGAAGt                              <  2:6331458/151‑1 (MQ=255)
                                                                                                                              ccgcACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGATCCTCCa                     <  2:7012762/151‑1 (MQ=255)
                                                                                                                              ccgcACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGATCCTCCa                     >  1:8174555/1‑151 (MQ=255)
                                                                                                                              ccgcACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGATCCTCCa                     <  1:4304006/151‑1 (MQ=255)
                                                                                                                                gcACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGATCCTCCAgg                   <  2:7770391/151‑1 (MQ=255)
                                                                                                                                 cACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGATCCTCCAGGc                  <  2:9575268/151‑1 (MQ=255)
                                                                                                                                 cACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGATCCTCCAGGc                  >  1:6788135/1‑151 (MQ=255)
                                                                                                                                 cACGCGGGAAGCCCTGGCCCTAGCCTTCCCGCGTGCTGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGGTGGGTGAACCAGCGTGGACCTTCCTCCGTCATGTAGAAGTTATCCTCCAGGc                  >  2:4319283/1‑151 (MQ=255)
                                                                                                                                     cgGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCTGCGTGGTCCTTCCTCCGTCATGTAGAAGTGATCCTCCAGGCGAAt              <  2:8919824/151‑1 (MQ=255)
                                                                                                                                     cgGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGATCCTCCAGGCGAAt              <  2:5291925/151‑1 (MQ=255)
                                                                                                                                     cgGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGATCCTCCAGGCGAAt              >  2:8492209/1‑151 (MQ=255)
                                                                                                                                      gggAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGATCCTCCAGGCGAATg             <  1:8574483/151‑1 (MQ=255)
                                                                                                                                         aaGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGATCCTCCAGGCGAATGCCg          >  2:4155555/1‑151 (MQ=255)
                                                                                                                                              cTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGATCCTCCAGGCGAATGCCGAACTc     <  1:3977013/151‑1 (MQ=255)
                                                                                                                                              cTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGATCCTCCAGGCGAATGCCGAACTc     >  2:1834107/1‑151 (MQ=255)
                                                                                                                                              cTGGCCCTTGCCGTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCTGCCGCCGATGACGAAGGGATCGTCGTTGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGATCCTCCAGGCGAATGCCGAACTc     <  2:903387/151‑1 (MQ=255)
                                                                                                                                                 gCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGATCCTCCAGGCGAATGCCGAACTCGcc  >  1:2212746/1‑151 (MQ=255)
                                                                                                                                                 gCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGATCCTCCAGGCGAATGCCGAACTCGcc  <  2:618051/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GATGCCCTGGTTGCCGCTGGCGCATCCGTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCTTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGATCCTCCAGGCGAATGCCGAACTCGCC  >  NC_008463/5198257‑5198552

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