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.5% 1462.4 / NA 388intergenic (+31/+37)PA14_RS23800/mdpAABC transporter substrate‑binding protein/metallopeptidase MdpA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (181/204);  minor base G (0/2);  total (182/206)
Fisher's exact test for biased strand distribution p-value = 5.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 390 total aligned reads displayed.

GATGCCCTGGTTGCCGCTGGCGCATCCGTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCTTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGATCCTCCAGGCGAATGCCGAACTCGCC  >  NC_008463/5198257‑5198552
                                                                                                                                                     |                                                                                                                                                   
gATGCCCTGGTTGCCGCTGGCGCATCCGTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCctt                                                                                                                                                   >  2:2437599/1‑151 (MQ=255)
gATGCCCTGGTTGCCGCTGGCGCATCCGTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCctt                                                                                                                                                   <  2:771314/151‑1 (MQ=255)
 aTGCCCTGGTTGCCGCTGGCGCATCCGTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTg                                                                                                                                                  <  1:4298816/151‑1 (MQ=255)
  tGCCCTGGTTGCCGCTGGCGCATCCGTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGc                                                                                                                                                 <  1:8929540/151‑1 (MQ=255)
  tGCCCTGGTTGCCGCTGGCGCATCCGTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGc                                                                                                                                                 <  2:4349508/151‑1 (MQ=255)
   gCCCTGGTTGCCGCTGGCGCATCCGTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGcc                                                                                                                                                <  1:3597550/151‑1 (MQ=255)
     ccTGGTTGCCGCTGGCGCATCCGTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCGTGCCtt                                                                                                                                              >  1:1366022/1‑151 (MQ=255)
      cTGGTTGCCGCTGGCGCATCCGTTTATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCAAGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGACGCCCTGGCCCTTGCCTTc                                                                                                                                             >  2:5781703/1‑151 (MQ=255)
       tGGTTGCCGCTGGCGCATCCGTTGATGTTCGACGCGCAACGCAAGTACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTcc                                                                                                                                            >  2:9606927/1‑151 (MQ=255)
       tGGTTGCCGCTGGCGCATCCGTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCGTGCCTTcc                                                                                                                                            >  1:2153010/1‑151 (MQ=255)
         gTTGCCGCTGGCGCATCCGTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCcg                                                                                                                                          <  1:99993/151‑1 (MQ=255)
             ccGCTGGCGCATCCGTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTg                                                                                                                                      >  1:8729894/1‑151 (MQ=255)
               gCTGGCGCATCCGTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTgcg                                                                                                                                    <  1:9551169/151‑1 (MQ=255)
               gCTGGCGCATCCGTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTgcg                                                                                                                                    >  2:5049239/1‑151 (MQ=255)
               gCTGGCGCATCCGTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCGTCCCGCGGgcg                                                                                                                                    >  1:1902416/1‑151 (MQ=255)
                cTGGCGCATCCGTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCgcgggcgg                                                                                                                                   >  1:702392/1‑151 (MQ=255)
                  ggCGCATCCGTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTgcggcg                                                                                                                                 <  2:8756159/151‑1 (MQ=255)
                       aTCCGTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCGACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGACGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATc                                                                                                                            <  2:4536598/151‑1 (MQ=255)
                         ccGTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCATGCCTTCCCGCGTGCGGCGTTATCGt                                                                                                                          >  1:2324579/1‑151 (MQ=255)
                            ttGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTcg                                                                                                                       >  2:524044/1‑151 (MQ=255)
                            ttGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGAGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTcg                                                                                                                       <  2:4032425/151‑1 (MQ=255)
                             tGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTcgg                                                                                                                      >  2:3760375/1‑151 (MQ=255)
                              gATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCGTCCCGCGTGCGGCGTTATCGTTcggc                                                                                                                     >  1:4890363/1‑151 (MQ=255)
                                  ttCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTcggccggc                                                                                                                 >  2:10609261/1‑151 (MQ=255)
                                    cGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTc                                                                                                               <  2:6755339/151‑1 (MQ=255)
                                     gACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCa                                                                                                              <  2:2275957/151‑1 (MQ=255)
                                       cgcgCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAgc                                                                                                            <  1:2664667/151‑1 (MQ=255)
                                        gcgcAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAgcc                                                                                                           <  2:10660910/151‑1 (MQ=255)
                                         cgcaacgcaaGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAgccg                                                                                                          <  1:10450404/151‑1 (MQ=255)
                                            aacgcaaGTACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCGGccgccg                                                                                                       >  1:10560719/1‑151 (MQ=255)
                                              cgcaaGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGTAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCAGGCCCTTGCCGGCCCGCGTGCGGCGGTATCGTTCGGCCGGCTCTGCCGCCGAg                                                                                                     >  2:1264975/1‑151 (MQ=255)
                                              cgcaaGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCGTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCCAg                                                                                                     >  1:10255376/1‑151 (MQ=255)
                                                   gAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGa                                                                                                >  2:8565848/1‑151 (MQ=255)
                                                     aCGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAg                                                                                              <  2:4242663/151‑1 (MQ=255)
                                                      cGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGGGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAgg                                                                                             >  1:8429125/1‑151 (MQ=255)
                                                        tCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCCTGCCATCCCGCGTGCGGCGTTATCGTTCGGCCGGCGCAGCCGCCGCGGCCGGAGGGa                                                                                           >  1:3826310/1‑151 (MQ=255)
                                                          aGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGGAGCCCTGGCCCTTGCCGTCCCGCGTGCGGCGTTATAGTTCGGGCGGCTGAGCCGCCGAGGCCGAAGGGAtc                                                                                         >  1:7291058/1‑151 (MQ=255)
                                                          aGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGAtc                                                                                         >  2:6904703/1‑151 (MQ=255)
                                                          aGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACCCGGGAAGCCCAGGCCCATGCCATCCCCCGTGCGGCGGTATCGTTCGGGCGGCTCAGCCCCCGAGGCCGAAGGGAtc                                                                                         >  1:160843/1‑151 (MQ=255)
                                                            cGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGAtcgt                                                                                       <  1:1193401/151‑1 (MQ=255)
                                                            cGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCATGGCTTCCCGCGTGCGGCGTTTTAGTTAGGCCGGCGCAGCCGCCGAGGCCGAAGGGAtcgg                                                                                       >  1:7416662/1‑150 (MQ=255)
                                                             ggCTACCGCATGAGCCCGATGTCGGCGCGGGAGTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGGGGGAAGCCCTGGCCCGTGCCTTCCCGCGTGCGGGGGTATCGTTCGGCCGGCTCAGCCGCCGTGGCCGAAGGGAtcgtc                                                                                      >  1:5706580/1‑151 (MQ=255)
                                                                tACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGGACGCGGGTAGCCCAGGCCCTTGCCTTCCCGGGGGCGGGGGTATAGGTAGGGCGGCTCAGCCCCCGGTGCCGAAGGGAtagtcgct                                                                                   >  1:107259/1‑149 (MQ=255)
                                                                tACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCGTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGAt                                                                                   <  2:3096734/151‑1 (MQ=255)
                                                                tACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCATGCCATCCCGCGGGCGGCGTTATCGTTAGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGAt                                                                                   >  1:5077533/1‑151 (MQ=255)
                                                                   cGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATgga                                                                                <  1:5241243/151‑1 (MQ=255)
                                                                   cGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCATCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATgga                                                                                >  2:2844735/1‑151 (MQ=255)
                                                                    gCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATggac                                                                               <  2:4874368/151‑1 (MQ=255)
                                                                     cATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATggacg                                                                              >  2:1090083/1‑151 (MQ=255)
                                                                     cATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATggacg                                                                              >  2:199769/1‑151 (MQ=255)
                                                                         aGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACt                                                                          <  1:4064442/151‑1 (MQ=255)
                                                                                   cGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACCCGGGAAGCCCAGGGCCATTGCTTCCCCCGTGCGGCGGTTACGTTCGGCCGGCGCACCCCCCGAGGCCGAACGGATCGTCGATGGTCGGACACGggtgggtg                                                                >  1:5517383/1‑151 (MQ=255)
                                                                                    ggCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGa                                                               <  2:8490633/151‑1 (MQ=255)
                                                                                     gcgcGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGaa                                                              <  1:3291468/151‑1 (MQ=255)
                                                                                      cgcgGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGGTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAAc                                                             >  1:810281/1‑151 (MQ=255)
                                                                                       gcgGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAAcc                                                            >  2:642523/1‑151 (MQ=255)
                                                                                       gcgGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAAcc                                                            >  2:8931295/1‑151 (MQ=255)
                                                                                             cTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTg                                                      >  2:2609750/1‑151 (MQ=255)
                                                                                              ttttCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTgg                                                     <  2:9521266/151‑1 (MQ=255)
                                                                                              ttttCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTgg                                                     <  1:7126814/151‑1 (MQ=255)
                                                                                                 tCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGtcc                                                  <  1:9645532/151‑1 (MQ=255)
                                                                                                  cgcgGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCGTCCCCCGTGCGGCGTTATAGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGtcct                                                 >  1:5943195/1‑151 (MQ=255)
                                                                                                    cgGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGtccttc                                               <  2:1911309/151‑1 (MQ=255)
                                                                                                    cgGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCAGGCCCGTGCCGTACCCCGGGCGGCGTTATAGTTAGGCCGGCGAAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGTCACGGGTTGGTGAAACACCGTGGtccgtc                                               >  1:3145281/1‑151 (MQ=255)
                                                                                                     gggTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGtccttcc                                              <  2:5870096/151‑1 (MQ=255)
                                                                                                        tGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTtcctcc                                           <  1:6155085/151‑1 (MQ=255)
                                                                                                         gAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGGCCTtcctcgg                                          >  2:8340796/1‑149 (MQ=255)
                                                                                                          aaGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGGCCTTGCCTTCCCGCCTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGGTCGGCGATGGACGGACTCGGGTGGGTGACCCCGCGTGGTCCTTCCTCCgg                                         >  1:650117/1‑150 (MQ=255)
                                                                                                          aaGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTtcctcctt                                         <  2:8266215/151‑3 (MQ=255)
                                                                                                            gCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCGGCGTGCGGCGTTTTCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCa                                       <  2:9420778/151‑1 (MQ=255)
                                                                                                            gCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCa                                       <  1:3887608/151‑1 (MQ=255)
                                                                                                            gCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCa                                       <  1:9770092/151‑1 (MQ=255)
                                                                                                             cTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCAt                                      <  1:1390985/151‑1 (MQ=255)
                                                                                                              tGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATg                                     <  1:9323113/151‑1 (MQ=255)
                                                                                                              tGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATg                                     <  2:8029486/151‑1 (MQ=255)
                                                                                                                 aCTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAg                                  <  2:6107953/151‑1 (MQ=255)
                                                                                                                 aCTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAg                                  <  2:6273441/151‑1 (MQ=255)
                                                                                                                  cTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGa                                 <  1:759028/151‑1 (MQ=255)
                                                                                                                  cTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGa                                 <  2:7852754/151‑1 (MQ=255)
                                                                                                                     aTCGCTGCGCCGCACGCGGGAAGCCCAGGCCCGTGCCATCCCGCGGGCGGCGGTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACACGGGTTGGTGAACCAGCGTGGTCCGTCCAACGTCATGTTGAAGt                              >  2:1726731/1‑151 (MQ=255)
                                                                                                                         cTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGAtc                          <  2:9790078/151‑1 (MQ=255)
                                                                                                                           gcgcCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGAtcct                        <  2:7044296/151‑1 (MQ=255)
                                                                                                                            cgccgcACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGAtcctc                       >  1:7934301/1‑151 (MQ=255)
                                                                                                                             gccgcACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGAtcctcc                      <  1:8828644/151‑1 (MQ=255)
                                                                                                                             gccgcACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGAtcctcc                      <  2:1418909/151‑1 (MQ=255)
                                                                                                                               cgcACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGATCCTCCAg                    <  1:3455540/151‑1 (MQ=255)
                                                                                                                               cgcACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGATCCTCCAg                    <  1:8282879/151‑1 (MQ=255)
                                                                                                                                gcACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGATCCTCCAgg                   <  1:6637233/151‑1 (MQ=255)
                                                                                                                                     cgGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGATCCTCCAGGCGAAt              >  2:3640207/1‑151 (MQ=255)
                                                                                                                                     cgGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGATCCTCCAGGCGAAt              <  1:4426029/151‑1 (MQ=255)
                                                                                                                                        gAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGATCCTCCAGGCGAATGcc           >  2:3215952/1‑151 (MQ=255)
                                                                                                                                        gAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGGGGTCCTTCCTCCGTCATGTAGAAGTGATCCTCCAGGCGAATGcc           >  2:7625650/1‑151 (MQ=255)
                                                                                                                                          aGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGATCCTCCAGGCGAATGCCGa         >  1:10128743/1‑151 (MQ=255)
                                                                                                                                           gCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGATCCTCCAGGCGAATGCCGaa        <  1:8189802/151‑1 (MQ=255)
                                                                                                                                           gCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGATCCTCCAGGCGAATGCCGaa        >  2:7253048/1‑151 (MQ=255)
                                                                                                                                            cccTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGATCCTCCAGGCGAATGCCGAAc       >  1:6602976/1‑151 (MQ=255)
                                                                                                                                            cccTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGATCCTCCAGGCGAATGCCGAAc       <  2:5479787/151‑1 (MQ=255)
                                                                                                                                            cccTGGCCCTTGCCGTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGATCCTCCAGGCGAATGCCGAAc       <  2:3342693/151‑1 (MQ=255)
                                                                                                                                                tGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGATCCTCCAGGCGAATGCCGAACTCg    <  2:8650356/151‑1 (MQ=255)
                                                                                                                                                 ggCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGATCCTCCAGGCGAATGCCGAACTCGc   <  2:9208042/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: