Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,108,695 C→G 5.9% A315P (GCC→CCC) ‡ PP_5359 ← CobW/P47K family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,108,6950CG 5.9% 330.4 / 8.9 103A315P (GCC→CCC) ‡PP_5359CobW/P47K family protein
Reads supporting (aligned to +/- strand):  ref base C (45/50);  new base G (3/3);  total (49/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

Only 100 of 103 total aligned reads displayed.

CACAGCGTGCAGATCGAACTAGCGCCCTACCCGTTCGAAATCAAGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGACCACCAGCCCTCCATGCTGTGGATAACTCCCTTGG  >  NC_002947/6108549‑6108841
                                                                                                                                                  |                                                                                                                                                   
cacaGCGTGCAGATCGAACTAGCGCCCTACCCGTTCGAAATCAAGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCtt                                                                                                                                                   <  1:79508/149‑1 (MQ=255)
        gCAGATCGAACTAGCGCCCTACCCGTTCGAAATCAAGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAAcgcg                                                                                                                                          >  5:408969/1‑149 (MQ=255)
           gATCGAACTAGCGCCCTACCCGTTCGAAATCAAGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGggt                                                                                                                                       >  6:31544/1‑149 (MQ=255)
             tCGAACTAGCGCCCTACCCGTTCGAAATCAAGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCggg                                                                                                                                        >  3:349933/1‑146 (MQ=255)
             tCGAACTAGCGCCCTACCCGTTCGAAATCAAGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCggg                                                                                                                                        <  4:349933/146‑1 (MQ=255)
             tCGAACTAGCGCCCTACCCGTTCGAAATCAAGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGggttg                                                                                                                                     <  5:250129/149‑1 (MQ=255)
             tCGAACTAGCGCCCTACCCGTTCGAAATCAAGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTAGACACGCGggttg                                                                                                                                     >  6:406274/1‑149 (MQ=255)
               gAACTAGCGCCCTACCCGTTCGAAATCAAGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGggttgcg                                                                                                                                   <  3:72903/149‑1 (MQ=255)
               gAACTAGCGCCCTACCCGTTCGAAATCAAGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGggttgc                                                                                                                                    <  3:199590/148‑1 (MQ=255)
               gAACTAGCGCCCTACCCGTTCGAAATCAAGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGCGTTGCAAGCCGggttgaa                                                                                                                                   <  2:80751/149‑3 (MQ=255)
                   tAGCGCCCTACCCGTTCGAAATCAAGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGggttgc                                                                                                                                    <  6:280561/144‑1 (MQ=255)
                   tAGCGCCCTACCCGTTCGAAATCAAGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGggttgc                                                                                                                                    >  5:280561/1‑144 (MQ=255)
                     gcgcCCTACCCGTTCGAAATCAAGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGgttgcggttgcg                                                                                                                             >  6:193429/1‑149 (MQ=255)
                     gcgcCCTACCCGTTCGAAATCAAGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGgttgcggttgcg                                                                                                                             <  5:21469/149‑1 (MQ=255)
                     gcgcCCTACCCGTTCGAAATCAAGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACAGACGCGACATTCATTCAACCCGGCTTGCAACGCGggttgcggttg                                                                                                                               <  1:292776/147‑1 (MQ=255)
                      ctcccTTCCCGTTCGAAATCAAGACCCCACCACTCGCCGACAACCTCATAGAGATGCAGCAATGGCGCCAGGAGCTTAACCGCTGGCGCAACTGATCAACTGCCGCGACATTCATTCAACCCGGCTTGCAACGCGGGttgcggttgcgg                                                                                                                            <  7:224827/147‑1 (MQ=255)
                        cccTACCCGTTCGAAATCAAGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCt                                                                                                                          <  1:411030/149‑1 (MQ=255)
                          cTACCCGTTCGAAATCAAGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGa                                                                                                                        <  5:406274/149‑1 (MQ=255)
                          cTACCCGTTCGAAATCAAGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGa                                                                                                                        <  1:226545/149‑1 (MQ=255)
                             cccGTTCGAAATCAAGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTCGCAACGCGGGTTGCGGTTGCGGCTGATg                                                                                                                      <  1:232865/149‑1 (MQ=255)
                                 ttCGAAATCAAGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGaaaa                                                                                                                 >  3:152667/1‑149 (MQ=255)
                                    gAAATCAAGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAAt                                                                                                                >  1:212030/1‑147 (MQ=255)
                                    gAAATCAAGACCCCACCAAGCGCCGACAACGTCATAGAGAAGCAGCAATGGCGCCAGCAGAAAAACACCTGCCGCAACTAATCAACAGACCCGACATTCATTCACCCCGGCTTGCAACGCCGGTTGCGGTTGCGGCTGAAGGAAAAt                                                                                                                >  8:320651/1‑147 (MQ=255)
                                     aaaTCAAGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAg                                                                                                             >  5:52352/1‑149 (MQ=255)
                                     aaaTCAAGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAg                                                                                                             <  5:31544/149‑1 (MQ=255)
                                          aaGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGAt                                                                                                                       <  4:207793/134‑1 (MQ=255)
                                          aaGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGAt                                                                                                                       >  3:207793/1‑134 (MQ=255)
                                           aGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGc                                                                                                                           >  6:92885/1‑129 (MQ=255)
                                           aGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGc                                                                                                                           <  5:92885/129‑1 (MQ=255)
                                            gACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAAGTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATa                                                                                                      >  1:130460/1‑149 (MQ=255)
                                            gACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATa                                                                                                      <  5:193429/149‑1 (MQ=255)
                                               cccACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGc                                                                                                                                                <  2:349755/104‑1 (MQ=255)
                                               cccACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACgg                                                                                                   <  7:171590/149‑1 (MQ=255)
                                                 caccacGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCt                                                                                                 <  1:260390/149‑1 (MQ=255)
                                                  accacGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTg                                                                                                >  7:136297/1‑149 (MQ=255)
                                                  accacGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTg                                                                                                >  1:65108/1‑149 (MQ=255)
                                                    cacGCGCAGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTc                                                                                              <  8:378408/149‑1 (MQ=255)
                                                            gACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTgcgc                                                                                       >  4:104764/1‑148 (MQ=255)
                                                               aaCGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGt                                                                                   <  4:152667/149‑1 (MQ=255)
                                                                    cATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACgg                                                                                                   <  7:198255/128‑1 (MQ=255)
                                                                    cATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACgg                                                                                                   >  8:198255/1‑128 (MQ=255)
                                                                       agagATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTc                                                                                  >  5:304938/1‑142 (MQ=255)
                                                                        gagaTGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCt                                                                          >  8:121773/1‑149 (MQ=255)
                                                                        gagaTGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCAGCGACTGGCGCAgggac                                                                          <  6:52352/149‑3 (MQ=255)
                                                                             gcagcaATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGcc                                                                       >  3:24740/1‑147 (MQ=255)
                                                                              cagcaATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCaa                                                                     <  8:109210/148‑1 (MQ=255)
                                                                              cagcaATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAg                                                                    <  7:121773/149‑1 (MQ=255)
                                                                               agcaATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGc                                                                   <  2:130460/149‑1 (MQ=255)
                                                                               agcaATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGCGTTGCAAGCCGggttgaa                                                                                                                                   <  5:152723/85‑3 (MQ=38)
                                                                               agcaATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGCGTTGCAAGCCGggttgaa                                                                                                                                   >  6:152723/1‑83 (MQ=38)
                                                                                gcaATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCg                                                                  <  2:65108/149‑1 (MQ=255)
                                                                                 caATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCgg                                                                 >  1:396369/1‑149 (MQ=255)
                                                                                  aaTGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCggg                                                                <  8:136297/149‑1 (MQ=255)
                                                                                   aTGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGc                                                               <  3:109224/149‑1 (MQ=255)
                                                                                     ggCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACgg                                                                                                   <  4:73846/111‑1 (MQ=255)
                                                                                     ggCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACgg                                                                                                   >  3:73846/1‑111 (MQ=255)
                                                                                      gcgcCAGGAGCATAACCGCTGGCGCAACTGAGCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCat                                                             <  7:320651/148‑1 (MQ=255)
                                                                                        gcCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATc                                                          >  5:373855/1‑149 (MQ=255)
                                                                                             gAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGcc                                                                       >  7:295108/1‑131 (MQ=255)
                                                                                             gAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGcc                                                                       <  8:295108/131‑1 (MQ=255)
                                                                                             gAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCAt                                                     >  5:245623/1‑149 (MQ=255)
                                                                                               gCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGc                                                               <  4:329726/137‑1 (MQ=255)
                                                                                               gCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGc                                                               >  3:329726/1‑137 (MQ=255)
                                                                                                 aTAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAg                                                 >  7:113828/1‑149 (MQ=255)
                                                                                                   aaCCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGTTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCa                                                  >  7:68457/1‑146 (MQ=255)
                                                                                                   aaCCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGTTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCa                                                  <  8:68457/146‑1 (MQ=255)
                                                                                                      cGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCAt                                            <  2:212030/149‑1 (MQ=255)
                                                                                                        cTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCAtt                                           >  7:88321/1‑148 (MQ=255)
                                                                                                                aCTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGacca                                  >  8:293683/1‑149 (MQ=255)
                                                                                                                       aaCTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCTGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCAATTGATTGACCACCCGccc                           >  2:100371/1‑149 (MQ=255)
                                                                                                                       aaCTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGACCACCAGccc                           >  3:148396/1‑149 (MQ=255)
                                                                                                                       aaCTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCGGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGACCACCAGcc                            >  8:150485/1‑148 (MQ=255)
                                                                                                                        aCTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGACCACCAGCCCt                          <  6:348349/149‑1 (MQ=255)
                                                                                                                        aCTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGACCACCAGCCCt                          >  5:348349/1‑149 (MQ=255)
                                                                                                                         cTGACGCGACATTCATTCAACCCGGCTTGCAACGAGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCACCATTCAGCCCATTGAATGACCACCAGCCCTc                         >  7:160608/1‑149 (MQ=255)
                                                                                                                          tGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGc                                                               <  5:81663/110‑1 (MQ=255)
                                                                                                                          tGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGc                                                               >  6:81663/1‑110 (MQ=255)
                                                                                                                              gcgACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGc                                                                        <  4:194887/97‑1 (MQ=255)
                                                                                                                              gcgACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGc                                                                        >  3:194887/1‑97 (MQ=255)
                                                                                                                                gACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGACCACCAGCCCTCCATGCtg                  <  7:150485/149‑1 (MQ=255)
                                                                                                                                 aCATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGACCACCAGCCCTCCATGCtgt                 <  7:293683/149‑1 (MQ=255)
                                                                                                                                   attcattcAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGACCACCAGCCCTCCATGCtgtg                >  1:172021/1‑148 (MQ=255)
                                                                                                                                     tcattcaACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGAt                                                                                                       <  4:306792/59‑1 (MQ=255)
                                                                                                                                     tcattcaACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGAt                                                                                                       >  3:306792/1‑59 (MQ=255)
                                                                                                                                       attcaACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGACCACCAGTCCTCCATGCTGTGGATaa           <  6:245623/149‑1 (MQ=255)
                                                                                                                                       attcaACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGACCACCAGCCCTCCATGCTGTGGATaa           <  4:148396/149‑1 (MQ=255)
                                                                                                                                       attcaACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGACCACCAGCCCTCCATGCTGTGGATaa           <  6:373855/149‑1 (MQ=255)
                                                                                                                                       attaaACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGACCACCAGCCCTCCATGCTGTGGATa            >  4:380421/1‑148 (MQ=255)
                                                                                                                                        ttcaACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGACCACCAGCCCTCCATGCTGTGtata            >  4:89888/1‑147 (MQ=255)
                                                                                                                                        ttcaACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGACCACCAGCCCTCCATGCTGTGGATaa           >  7:195723/1‑148 (MQ=255)
                                                                                                                                        cgcaACCCGCGTTGCAAGCCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGACCACCAGCCCTCCATGCTGTGGATAAc          >  6:71829/3‑149 (MQ=255)
                                                                                                                                         tcaACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGaa                                        >  2:302156/1‑118 (MQ=255)
                                                                                                                                         tcaACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGaa                                        <  1:302156/118‑1 (MQ=255)
                                                                                                                                           aaCCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGacc                                   <  5:277336/121‑1 (MQ=255)
                                                                                                                                           aaCCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGacc                                   >  6:277336/1‑121 (MQ=255)
                                                                                                                                           aaCCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGACCACCAGCCCTCCATGCTGTGGATAACTc        >  3:7285/1‑148 (MQ=255)
                                                                                                                                           aaCCCGCGTTGCAAGCCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGACCACCAGCCCTCCATGCTGTGGATAACTcc       >  7:8195/1‑149 (MQ=255)
                                                                                                                                             cccGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGACTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGACCACCAGCCCt                          <  2:189937/128‑1 (MQ=255)
                                                                                                                                             cccGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGACTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGACCACCAGCCCt                          >  1:189937/1‑128 (MQ=255)
                                                                                                                                                ggCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGACCACCAGCCCTCCATGCTGTGGATAACTCCCTTgg  <  4:24740/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CACAGCGTGCAGATCGAACTAGCGCCCTACCCGTTCGAAATCAAGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGACCACCAGCCCTCCATGCTGTGGATAACTCCCTTGG  >  NC_002947/6108549‑6108841

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: