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% 1708.9 / NA 439intergenic (+31/+37)PA14_RS23800/mdpAABC transporter substrate‑binding protein/metallopeptidase MdpA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (198/240);  minor base A (0/1);  total (198/241)
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 441 total aligned reads displayed.

GATGCCCTGGTTGCCGCTGGCGCATCCGTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCTTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGATCCTCCAGGCGAATGCCGAACTCGCC  >  NC_008463/5198257‑5198552
                                                                                                                                                    |                                                                                                                                                   
gATGCCCTGGTTGCCGCTGGCGCATCCGTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCctt                                                                                                                                                   <  2:4824339/151‑1 (MQ=255)
 aTGCCCTGGTTGCCGCTGGCGCATCCGTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTg                                                                                                                                                  <  1:2130176/151‑1 (MQ=255)
   gCCCTGGTTGCCGCTGGCGCATCCGTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGcc                                                                                                                                                <  1:1916377/151‑1 (MQ=255)
    cccTGGTTGCCGCTGGCGCATCCGTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCt                                                                                                                                               >  1:386952/1‑151 (MQ=255)
    cccTGGTTGCCGCTGGCGCATCCGTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGAGGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGccg                                                                                                                                               >  2:8891275/1‑150 (MQ=255)
      cTGGTTGCCGCTGGCGCATCCGTTGATGTTCGACGCGCAACGCAAGCACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGCCTGTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCTCGCGGGAAGCCCTGGCCCTTGCCTTc                                                                                                                                             >  2:4497554/1‑151 (MQ=255)
      cTGGTTGCCGCTGGCGCATCCGTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTc                                                                                                                                             <  1:9748506/151‑1 (MQ=255)
        ggTTGCCGCTGGCGCATCCGGTGATGTTCGACGCGCAACGCAAGACCGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTGTAGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGACGCCCTGGCCCTTGCCTTccc                                                                                                                                           >  2:6089636/1‑151 (MQ=255)
         gTTGCCGCTGGCGCATCCGTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCcg                                                                                                                                          <  1:1254337/151‑1 (MQ=255)
         gTTGCCGCTGGCGAATCCGTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTAGGCGCGGGACTATTCGCGGGTGAAGCTGTACTGATCGCTGCGCCGCACGAGGGAAGCCCTGGCCCTGGCCTTCCcg                                                                                                                                          >  2:8652687/1‑151 (MQ=255)
          ttGCCGCTGGCGCATCCGTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCcgc                                                                                                                                         <  2:3598196/151‑1 (MQ=255)
           tGCCGCTGGCGCATCCGTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCcgcg                                                                                                                                        <  1:530313/151‑1 (MQ=255)
              cGCTGGCGCATCCGTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAATCCCTGGCCCTTGCCTTCCCGCGTgc                                                                                                                                     >  1:2295259/1‑151 (MQ=255)
               gCTGGCGCATCCGTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTgcg                                                                                                                                    >  1:1786712/1‑151 (MQ=255)
                cTGGCGCATCCGTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTgcgg                                                                                                                                   >  1:7210319/1‑151 (MQ=255)
                cTGGCGCATCCGTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTgcgg                                                                                                                                   <  1:8647868/151‑1 (MQ=255)
                  ggCGCATCCGTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTgcggcg                                                                                                                                 <  1:2688793/151‑1 (MQ=255)
                      cATCCGTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTAt                                                                                                                             >  2:6829028/1‑151 (MQ=255)
                        tCCGTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCg                                                                                                                           >  2:2247409/1‑151 (MQ=255)
                         ccGTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGt                                                                                                                          >  2:6554959/1‑151 (MQ=255)
                          cGTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGtt                                                                                                                         <  1:3413140/151‑1 (MQ=255)
                          cGTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCACGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGtt                                                                                                                         <  2:2933528/151‑1 (MQ=255)
                           gTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTc                                                                                                                        >  1:9069037/1‑151 (MQ=255)
                           gTTGATGTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATAGTTc                                                                                                                        >  1:4080384/1‑151 (MQ=255)
                                 gTTCGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCATGCCTTCCCGCGTGCGGCGGTATCGTTcggccgg                                                                                                                  >  1:9623690/1‑151 (MQ=255)
                                    cGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTc                                                                                                               <  2:7644886/151‑1 (MQ=255)
                                    cGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTc                                                                                                               <  2:128503/151‑1 (MQ=255)
                                    cGACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGGTCGGCCGGCTc                                                                                                               >  1:298577/1‑151 (MQ=255)
                                     gACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCa                                                                                                              <  1:9586147/151‑1 (MQ=255)
                                     gACGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCa                                                                                                              <  2:8101007/151‑1 (MQ=255)
                                      aCGCGCAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGGGCGGCGTTATCGTTCGGCCGGCTccg                                                                                                             >  1:688894/1‑149 (MQ=255)
                                        gcgcAACGCAAGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAgcc                                                                                                           <  2:7516834/151‑1 (MQ=255)
                                           caacgcaaGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAgccgcc                                                                                                        <  1:10341156/151‑1 (MQ=255)
                                            aacgcaaGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGGTCGGCCGGCTCAGccgccg                                                                                                       >  1:6181603/1‑151 (MQ=255)
                                               gcaaGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGCgg                                                                                                    >  1:1683138/1‑151 (MQ=255)
                                               gcaaGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAgg                                                                                                    <  1:2428494/151‑1 (MQ=255)
                                                caaGAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCAGGGCCGTGCCGTCCCGCGTGCGGGGGTTTAGTTAGGGCGGCGCAGCCGCCGAGGc                                                                                                   >  1:8423013/1‑151 (MQ=255)
                                                   gAACGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGa                                                                                                <  1:5490788/151‑1 (MQ=255)
                                                     aCGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAg                                                                                              >  2:7346504/1‑151 (MQ=255)
                                                      cGTCAGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAgg                                                                                             <  1:7472816/151‑1 (MQ=255)
                                                          aGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGAtc                                                                                         <  2:10758724/151‑1 (MQ=255)
                                                          aGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGAtc                                                                                         >  1:9988134/1‑151 (MQ=255)
                                                          aGCGGCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGAtc                                                                                         <  2:9368501/151‑1 (MQ=255)
                                                              gCTACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGAtcgtcg                                                                                     <  1:2844016/151‑1 (MQ=255)
                                                                tACCGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGAt                                                                                   <  1:10587815/151‑1 (MQ=255)
                                                                  ccGCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATgg                                                                                 <  1:4270902/151‑1 (MQ=255)
                                                                    gCATGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATggac                                                                               <  2:10390551/151‑1 (MQ=255)
                                                                      aTGAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATggacgg                                                                             <  2:4124830/151‑1 (MQ=255)
                                                                        gAGCCCGATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCCCGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATggacggac                                                                           <  2:8398530/151‑1 (MQ=255)
                                                                              gATGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCgggt                                                                     >  1:10208028/1‑149 (MQ=255)
                                                                               aTGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTTCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCTGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTg                                                                    >  1:11113861/1‑151 (MQ=255)
                                                                                tGTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTgg                                                                   <  2:11222052/151‑1 (MQ=255)
                                                                                 gTCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTggg                                                                  <  1:7114843/151‑1 (MQ=255)
                                                                                  tCGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGt                                                                 <  1:5789557/151‑1 (MQ=255)
                                                                                   cGGCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTg                                                                >  2:11108062/1‑151 (MQ=255)
                                                                                    ggCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCATTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGa                                                               <  2:1867081/151‑1 (MQ=255)
                                                                                    ggCGCGGGACTTTTCGCGGGTGAAGCTGGACTGATCCCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCGCAGCCGCCGAGGCCGAAGGGATCGGCGATGGACGGACTCGGGTGGGTGa                                                               >  2:4257418/1‑151 (MQ=255)
                                                                                      cgcgGGACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAAc                                                             <  2:9187996/151‑1 (MQ=255)
                                                                                          ggACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGc                                                         <  2:8193645/151‑1 (MQ=255)
                                                                                          ggACTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGc                                                         <  1:7593939/151‑1 (MQ=255)
                                                                                            aCTTTTCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGt                                                       <  1:6503422/151‑1 (MQ=255)
                                                                                              tgttCGCGGGTGAAGCTGGACTGATCGCTGCGCCTCACGCGGGAAGCCCTTGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTgg                                                     <  2:3274723/149‑1 (MQ=255)
                                                                                               tttCGCGGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGGTGGGTGAACCAGCGTGGt                                                    >  2:7078255/1‑151 (MQ=255)
                                                                                                   gcgGGTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGtccgt                                                >  2:5781626/1‑149 (MQ=255)
                                                                                                     gggTGAAGCTGGACTGATCGCTGCGCCGCCCGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGGTATCGGTCGGGCGGCTCAGCCGCCGAGGCCGCAGGGATCGTCGATGGACGGACTCGGGTGGGTGAACCAGCGGGGtccttcc                                              >  2:6643460/1‑151 (MQ=255)
                                                                                                     gggTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGtccttcc                                              <  2:6154407/151‑1 (MQ=255)
                                                                                                     gggTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGtccttcc                                              <  1:9380834/151‑1 (MQ=255)
                                                                                                      ggTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGtccttcct                                             <  2:5516246/151‑1 (MQ=255)
                                                                                                      ggTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGGTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGATGGGATCGGCGATGGACGGACGCGGTTGGGTGAACCAGCGTGGtccttcct                                             >  2:4613789/1‑151 (MQ=255)
                                                                                                       gTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTtcctc                                            <  1:3755083/151‑1 (MQ=255)
                                                                                                       gTGAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTtcctc                                            <  2:9562221/151‑1 (MQ=255)
                                                                                                         gAAGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCg                                          <  2:5590114/151‑1 (MQ=255)
                                                                                                           aGCTGGACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTc                                        <  1:6762967/151‑1 (MQ=255)
                                                                                                               ggACTGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGt                                    <  1:7756942/151‑1 (MQ=255)
                                                                                                                   tGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGaa                                <  1:3086316/151‑1 (MQ=255)
                                                                                                                   tGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGAGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGaa                                <  2:9830520/151‑1 (MQ=255)
                                                                                                                   tGATCGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCGTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGaa                                <  1:4923599/151‑1 (MQ=255)
                                                                                                                       cGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGATGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGa                            <  1:5859169/151‑1 (MQ=255)
                                                                                                                       cGCTGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGa                            <  2:8541247/151‑1 (MQ=255)
                                                                                                                       cGCTGCGCCGCACGCGGGAAGCCCAGGCCCATGCCATACCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCCAAGGGCACGTCGATGGACCGACACGGGTGGGTGAACCAGCGTGGTCCGTCCGCCGTCATGTAGAAGTGa                            >  1:9139354/1‑151 (MQ=255)
                                                                                                                          tGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGAtcc                         >  1:10019169/1‑151 (MQ=255)
                                                                                                                          tGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGAtcc                         >  2:10966762/1‑151 (MQ=255)
                                                                                                                          tGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGAtcc                         <  1:2592499/151‑1 (MQ=255)
                                                                                                                          tGCGCCGCACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACGCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGAtcc                         >  2:2467792/1‑151 (MQ=255)
                                                                                                                            cgccgcACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGAtcctc                       <  1:8723249/151‑1 (MQ=255)
                                                                                                                            cgccgcACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGAtcctc                       <  1:8071169/151‑1 (MQ=255)
                                                                                                                            cgccgcACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGAtcctc                       <  1:4551157/151‑1 (MQ=255)
                                                                                                                             gccgcACGCGGGAAGCCCTGGCCCTTGCCGTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGAtcctcc                      <  1:4847037/151‑1 (MQ=255)
                                                                                                                                gcACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGATCCTCCAgg                   <  2:335823/151‑1 (MQ=255)
                                                                                                                                 cACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGATCCTCCAGGc                  <  1:6617971/151‑1 (MQ=255)
                                                                                                                                 cACGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGATCCTCCAGGc                  >  1:3340654/1‑151 (MQ=255)
                                                                                                                                  aCGCGGGAAGCCCTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGATCCTCCAGGCg                 <  1:6126368/151‑1 (MQ=255)
                                                                                                                                            cccTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGATCCTCCAGGCGAATGCCGAAc       <  2:1741644/151‑1 (MQ=255)
                                                                                                                                            cccTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGATCCTCCAGGCGAATGCCGAAc       <  2:4995125/151‑1 (MQ=255)
                                                                                                                                             ccTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGATCCTCCAGGCGAATGCCGAACt      >  2:1450641/1‑151 (MQ=255)
                                                                                                                                             ccTGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGATCCTCCAGGCGAATGCCGAACt      >  1:516699/1‑151 (MQ=255)
                                                                                                                                               tGGCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGATCCTCCAGGCGAATGCCGAACTCg    >  2:7910910/1‑151 (MQ=255)
                                                                                                                                                 gCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGATCCTCCAGGCGAATGCCGAACTCGcc  >  1:10741575/1‑151 (MQ=255)
                                                                                                                                                 gCCCTTGCCTTCCCGCGTGCGGCGTTATCGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTGGTCCTTCCTCCGTCATGTAGAAGTGATCCTCCAGGCGAATGCCGAACTCGcc  >  2:9021686/1‑151 (MQ=255)
                                                                                                                                                 gCCCTTGCCTTCCCGCGTGCGGCGTTAACGTTCGGCCGGCTCAGCCGCCGAGGCCGAAGGGATCGTCGATGGACGGACTCGGTTGGGTGAACCAGCGTTGTCCTTCCTCCGTCATCTAAAAGTGATCCTCCAGGCGAATGCCGAACTCGcc  >  2:2764784/1‑151 (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: