Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,140,776 T→C 6.5% V171A (GTC→GCC) ‡ PP_1897 → DNA internalization‑related competence protein ComEC/Rec2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,140,7760TC 6.5% 176.1 / 2.5 61V171A (GTC→GCC) ‡PP_1897DNA internalization‑related competence protein ComEC/Rec2
Reads supporting (aligned to +/- strand):  ref base T (32/25);  new base C (2/2);  total (34/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01

TGGCGGCTGGCCGTGACCTTGCAGCGCCCGAACGGGCTGCTGAACCCGCATGGGCCCGACCGTGAAGCGCAGTTGCTGGCGCGACGGGTGGGGGCTACGGGTACGGTCAAGGCGGGGCAGTTGCTGGCGCCAGTCGCCGGTGGTTGGCGCGATGCGCTGCGTCAGCGCTTGCTGATGGTCGAAGCCAATGGCCGGCAGGCAGCGCTGGTGGCACTGGTGCTTGGCGACGGTGCGGGCCTGGCCCGTGAGGACTGGCAGACGTTGCAGGCCACCGGCACGGTG  >  NC_002947/2140643‑2140924
                                                                                                                                     |                                                                                                                                                    
tGGCGGCTGGCCGTGACCTTGCAGCGCCCGAACGGGCTGCTGAACCCGCATGGGCCCGACCGTGAAGCGCAGTTGCTGGCGCGACGGGTGGGGGCTACGGGTACGGACAAGGCGGGGCAGTTGCTGGCGCCAGTCGCCGGTGGTTGgcg                                                                                                                                       >  1:256144/1‑149 (MQ=255)
tGGCGGCTGGCCGTGACCTTGCAGCGCCCGAACGGGCTGCTGAACCCGCATGGGCCCGACCGTGAAGCGCAGGTGCTGGCGCGACGGGTGGGGGCTACGGGTACGGTCAAGGCGGGGCAGTTGCTGGCGCCAGTCGCCGGTGGTTGgcg                                                                                                                                       <  5:11553/149‑1 (MQ=255)
       tGGCCGTGACCTTGCAGCGCCCGAACGGGCTGCTGAACCCGCATGGGCCCGACCGTGAAGCGCAGTTGCTGGCGCGACGGGTGGGGGCTACGGGTACGGTCAAGGCGGGGCAGTTGCTGGCGCCAGTCGCCGGTGGTTGGCGCGATgcg                                                                                                                                >  2:126824/1‑149 (MQ=255)
         gCCGTGACCTTGCAGCGCCCGAACGGGCTGCTGAACCCGCATGGGCCCGACCGTGAAGCGCAGTTGCTGGCGCGACGGGTGGGGGCTACGGGTACGGTCAAGGCGGGGCAGTTGCTGGCGCCAGTCGCCGGTGGTTGGCGCGATGCGCt                                                                                                                              <  5:43262/149‑1 (MQ=255)
             tGACCTTGCAGCGCCCGAACGGGCTGCTGAACCCGCATGGGCCCGACCGTGAAGCGCAGTTGCTGGCGCGACGGGTGGGGGCTACGGGTACGGTCAAGGCGGGGCAGTTGCTGGCGCCAGTCGCCGGTGGTTGGCGCGATGCGCTGCGt                                                                                                                          <  7:38862/149‑1 (MQ=255)
               agcTTGCAGCGCCGGATCGGGCTGCTGAACCCGCATGGGCGCTACGGTGAAGCGGTGTTGCTTGCGCGACGGGTGTGGGCTACGGGTACGGTCAAGGGGGGGCAGTTGCTGGCGCCTGTCGGCGGTGGTTGGCGCGATGCGCTGCGTCa                                                                                                                        <  3:193138/147‑1 (MQ=255)
               aCCTTGCAGCGCCCGAACGGGCTGCTGAACCCGCATGGGCCCGACCGTGAAGCGCAGTTGCTGGCGCGACGGGTGGGGGCTACGGGTACGGTCAAGGCGGGGCAGTTGCTGGCGCCAACCACCGGCGACTGgcgcca                                                                                                                                    >  5:103553/1‑135 (MQ=255)
               aCCTTGCAGCGCCCGAACGGGCTGCTGAACCCGCATGGGCCCGACCGTGAAGCGCAGTTGCTGGCGCGACGGGTGGGGGCTACGGGTACGGTCAAGGCGGGGCAGTTGCTGGCGCCAACCACCGGCGACTGgcgcca                                                                                                                                    >  5:99157/1‑135 (MQ=255)
               aCCTTGCAGCGCCCGAACGGGCTGCTGAACCCGCATGGGCCCGACCGTGAAGCGCAGTTGCTGGCGCGACGGGTGGGGGCTACGGGTACGGTCAAGGCGGGGCAGTTGCTGGCGCCAACCACCGGCGACTGgcgcca                                                                                                                                    <  6:103553/137‑3 (MQ=255)
               aCCTTGCAGCGCCCGAACGGGCTGCTGAACCCGCATGGGCCCGACCGTGAAGCGCAGTTGCTGGCGCGACGGGTGGGGGCTACGGGTACGGTCAAGGCGGGGCAGTTGCTGGCGCCAACCACCGGCGACTGgcgcca                                                                                                                                    <  6:99157/137‑3 (MQ=255)
                      aGCGCCCGAACGGGCTGCTGAACCCGCATGGGCCCGACCGTGAAGCGCAGTTGCTGGCGCGACGGGTGGGGGCTACGGGTACGGTCAAGGCGGGGCAGTTGCTGGCGCCAGTCGCCGGTGGTTGGCGCGATGCGCTGCGTCAGCGCTTg                                                                                                                 >  4:112951/1‑149 (MQ=255)
                      aGCGCCCGAACGGGCTGCTGAACCCGCATGGGCCCGACCGTGAAGCGCAGTTGCTGGCGCGACGGGTGGGGGCTACGGGTACGGTCAAGGCGGGGCAGTTGCTGGCGCCAGTCGCCGGTGGTTGGCGCGATGCGCTGACGCAgcgcatc                                                                                                                 <  6:122370/149‑4 (MQ=255)
                          cccGAACGGGCTGCTGAACCCGCATGGGCCCGACCGTGAAGCGCAGTTGCTGGCGCGACGGGTGGGGGCTACGGGTACGGTCAAGGCGGGGCAGTTGCTGGCGCCAGTCGCCGGTGGTTGGCGCGATGCGCTGCGTCAGCGCTTGCTGa                                                                                                             >  7:153277/1‑149 (MQ=255)
                                   gctgctGAACCCGCATGGGCCCGACCGTGAAGCGCAGTTGCTGGCGCGACGGGTGGGGGCTACGGGTACGGTCAAGGCGGGGCAGTTGCTGGCGCCAGTCGCCGGTGGTTGGCGCGATGCGCTGCGTCAGCGCTTGCTGATGGTCGAAg                                                                                                    >  3:246650/1‑149 (MQ=255)
                                     tgctgAACCCGCATGGGCCCGACCGTGAAGCGCAGTTGCTGGCGCGACGGGTGGGGGCTACGGGTACGGTCAAGGCGGGGCAGTTGCTGGCGCCAGTCGCCGGTGGTTGGCGCGATGCGCTGCGTCAGCGCTTGCTGATGGTCGAAGc                                                                                                   >  3:131203/1‑148 (MQ=255)
                                      gctgAACCCGCATGGGCCCGACCGTGAAGCGCAGTTGCTGGCGCGACGGGTGGGGGCTACGGGTACGGTCAAGGCGGGGCAGTTGCTGGCGCCAGTCGCCGGTGGTTGGCGCGATGCGCTGCGTCAGCGCTTGCTGATGGTCGAAGcc                                                                                                  >  5:163770/1‑148 (MQ=255)
                                       ctgAACCCGCATGGGCCCGACCGTGAAGCGCAGTTGCTGGCGCGACGGGTGGGGGCTACGGGTACGGTCAAGGCGGGGCAGTTGCTGGCGCCAGTCGCCGGTGGTTGGCGCGATGCGCTGCGTCAGCGCTTGCTGATGGTCGAAGCCaa                                                                                                >  4:83005/1‑149 (MQ=255)
                                        tgAACCCGCATGGGCCCGACCGTGAAGCGCAGTTGCTGGCGCGACGGGTGGGGGCTACGGGTACGGTCAAGGCGGGGCAGTTGCTGGCGCCAGTCGCCGGTGGTTGGCGCGATGCGCTGCGTCAGCGCTTGCTGATGGTCGAAGCCAAt                                                                                               <  5:86904/149‑1 (MQ=255)
                                         gAACCCGCATGGGCCCGACCGTGAAGCGCAGTTGCTGGCGCGACGGGTGGGGGCTACGGGTACGGTCAAGGCGGGGCAGTTGCTGGCGCCAGTc                                                                                                                                                     <  2:256154/94‑1 (MQ=255)
                                         gAACCCGCATGGGCCCGACCGTGAAGCGCAGTTGCTGGCGCGACGGGTGGGGGCTACGGGTACGGTCAAGGCGGGGCAGTTGCTGGCGCCAGTc                                                                                                                                                     >  1:256154/1‑94 (MQ=255)
                                           aCCCGCATGGGCCCGACCGTGAAGCGCAGTTGCTGGCGCGACGGGTGGGGGCTACGGGTACGGTCAAGGCGGGGCAGTTGCTGGCGCCAGTCGCCGGTGGTTGGCGCGATGCGCTGCGTCAGCGCTTGCTGATGGTCGAAGCCAATGGc                                                                                            <  3:83005/149‑1 (MQ=255)
                                           aCCCGCATGGGCCCGACCGTGAAGCGCAGTTGCTGGCGCGACGGGTGGGGGCTACGGGTACGGTCAAGGCGGGGCAGTTGCTGGCGCCAGTCGCCGGTGGTTGGCGCGATGCGCTGCGTCAGCGCTTGCTGATGGTCGAAGCCAATGGc                                                                                            <  8:153277/149‑1 (MQ=255)
                                             ccGCATGGGCCCGACCGTGAAGCGCAGTTGCTGGCGCGACGGGTGGGGGCTACGGGTACGGTCAAGGCGGGGCAGTTGCTGGCGCCAGTCGCCGGTGGTTGGCGCGATGCGCTGCGTCAGCGCTTGCTGATGGTCGAAGCCAATGGCCg                                                                                          >  5:64922/1‑149 (MQ=255)
                                              cGCATGGGCCCGACCGTGAAGCGCAGTTGCTGGCGCGACGGGTGGGGGCTACGGGTACGGTCAAGGCGGGGCAGTTGCTGGCGCCAGTCGCCGGTGGTTGGCGCGATGCGCTGCGTCAGCGCTTGCTGATGGTCGAAGCCAATGGCCgg                                                                                         <  3:112951/149‑1 (MQ=255)
                                               gCATGGGCCCGCCCGTGAAGCGCAGTTGCTGGCGCGACGGGTGGGGGCTACGGGTACGGTCAAGGCGGGGCAGTTGCTGGCGCCAGTCGCCGGTGGTTGGCGCGATGCGCTGCGTCAGCGCTTGCTGATGGTCGAAGCCAATGGCCgg                                                                                         >  6:24423/1‑148 (MQ=255)
                                               gCATGGGCCCGACCGTGAAGCGCAGTTGCTGGCGCGACGGGTGGGGGCTACGGGTACGGTCAAGGCGGGGCAGTTGCTGGCGCCAGTCGCCGGTGGTTGGCGCGATGCGCTGCGTCAGCGCTTGCTGATGGTCGAAGCCAATGGCCgg                                                                                         >  4:125852/1‑148 (MQ=255)
                                               gCATGGGCCCGACCGTGAAGCGCAGTTGCTGGCGCGACGGGTGGGGGCTACGGGTACGGTCAAGGCGGGGCAGTTGCTGGCGCCAGTCGCCGGTGGTTGGCGCGATGCGCTGCGTCAGCGCTTGCTGATGGTCGAAGCCAATGGCCgg                                                                                         >  5:49101/1‑148 (MQ=255)
                                                cATGGGCCCGACCGTGAAGCGCAGTTGCTGGCGCGACGGGTGGGGGCTACGGGTACGGTCAAGGCGGGGCAGTTGCTGGCGCCAGTCGCCGGTGGTTGGCGCGATGCGCTGCGTCAGCGCTTGCTGATGGTCGAAGCCAATGGCCgg                                                                                         >  4:52222/1‑147 (MQ=255)
                                                      cccGACCGTGAAGCGCAGTTGCTGGCGCGACGGGTGGGGGCTACGGGTACGGTCAAGGCGGGGCAGTTGCTGGCGCCAGTCGCCGGTGGTTGGCGCGATGCGCTGCGTCAGCGCTTGCTGATGGTCGAAGCCAATGGCCGGCAGGCAgc                                                                                 >  5:110311/1‑149 (MQ=255)
                                                      cccGACCGTGAAGCGCAGTTGCTGGCGCGACGGGTGGGGGCTACGGGTACGGTCAAGGCGGGGCAGTTGCTGGCGCCAGTCGCCGGTGGTTGGCGCGATGCGCTGCGTCAGCGCTTGCTGATGGTCGAAGCCAATGGCCGGCAGGCAgc                                                                                 >  1:10630/1‑149 (MQ=255)
                                                       ccGACCGTGAAGCGCAGTTGCTGGCGCGACGGGTGGGGGCTACGGGTACGGTCAAGGCGGGGCAGTTGCTGGCGCCAGTCGCCGGTGGTTGGCGCGATGCGCTGCGTCCGCGCTTGCTGATGGTCGACGCAAATGGCCGGCAGGCAgcg                                                                                >  6:173781/1‑149 (MQ=255)
                                                       ccGACCGTGAAGCGCAGTTGCTGGCGCGACGGGTGGGGGCTACGGGTACGGTCAAGGCGGGGCAGTTGCTGGCGCCAGTCGCCGGTGGTTGGCGCGATGCGCTGCGTCAGCGCTTGCTGATGGTCGAAGCCAATGGCCGGCAGGCAgcg                                                                                <  4:131203/149‑1 (MQ=255)
                                                            cGTGAAGCGCAGTTGCTGGCGCGACGGGTGGGGGCTACGGGTACGGTCAAGGCGGGGCAGTTGCTGGCGCCAGTCGCCGGTGGTTGGCGCGATGCGCTGCGTCAGCGCTTGCTGATGGTCGAAGCCAATGGCCGGCAGGCAGCGCtggt                                                                           <  2:256144/149‑1 (MQ=255)
                                                                aaGCGCAGTTGCTGGCGCGACGGGTGGGGGCTACGGGTACGGTCAAGGCGGGGCAGTTGCTGGCGCCAGTCGCCGGTGGTTGGCGCGATGCGCTGCGTCAGCGCTTGCTGATGGTCGAAGCCAATGGCCGGCAGGCAGCGCTGGTGGCa                                                                       <  6:49101/149‑1 (MQ=255)
                                                                aaGCGCAGTTGCTGGCGCGACGGGTGGGGGCTACGGGTACGGTCAAGGCGGGGCAGTTGCTGGCGCCAGTCGCCGGTGGTTGGCGCGATGCGCTGCGTCAGCGCTTGCTGATGGTCGAAGCCAATGGCCGGCAGGCAGCGCTGGTGGCa                                                                       >  1:135586/1‑149 (MQ=255)
                                                                       gTTGCTGGCGCGACGGGTGGGGGCTACGGGTACTGTCAAGGCGGGGCAGTTGCTGGCGCCAGTCGCCGGTGGTTGGCGCGATGCGCTGCGTCAGCGCTTGCTGATGGTCGAAGCCAATGGCCGGCAGGCAGCGCTGGTGGCACTGGTGc                                                                >  6:16248/1‑149 (MQ=255)
                                                                            tGGCGCGACGGGTGGGGGCTACGGGTACGGTCAAGGCGGGGCAGTTGCTGGCGCCAGTCGCCGGTGGTTGGCGCGATGCGCTGCGTCAGCGCTTGCTGATGGTCGAAGCCAATGGCCGGCAGGCAGCGCTGGTGGCACTGGTGCTTGGc                                                           >  2:152528/1‑149 (MQ=255)
                                                                              gcgcGACGGGTGGGGGCTACGGGTACGGTCAAGGCGGGGCAGTTGCTGGCGCCAGTCGCCGGTGGTTGGCGCGATGCGCTGCGTCAGCGCTTGCTGATGGTCGAAGCCAATGGCCGGCAGGCAGCGCTGGTGGCACTGGTGCTTGGc                                                           <  3:125852/147‑1 (MQ=255)
                                                                                     gggTGGGGGCTACGGGTACGGTCAAGGCGGGGCAGTTGCTGGCGCCAGTCGCCGGTGGTTGGCGCGATGCGCTGCGTCAGCGCTTGCTGATGGTCGAAGCCAATGGCCGGCAGGCAGCGCTGGTGGCACTGGTGCTTGGCGACGGTGCg                                                  <  1:126824/149‑1 (MQ=255)
                                                                                          ggggCTACGGGTACGGTCAAGGCGGGGCAGTTGCTGGCTCCAGTCGCCGGTGGTTGGCGCGATGCGCTGCGTCAGCGCTTGCTGATGGTCGAAGCCAATGGCCGGCAGGCAGCGCTGGTGGCACTGGTGCTTGGCGACGGTGCGGGcc                                              <  5:16248/148‑1 (MQ=255)
                                                                                             gCTACGGGTACGGTCAAGGCGGGGCAGTTGCTGGCGCCAGTCGCCGGTGGTTGGCGCGATGCGCTGCGTCAGCGCTTGCTGATGGTCGAAGCCAATGGCCGGCAGGCAGCGCTGGTGGCACTGGTGCTTGGCGACGGTGCGGGCCTGGc                                          <  6:163770/149‑1 (MQ=255)
                                                                                                   ggTACGGTCAAGGCGGGGCAGTTGCTGGCGCCAGTCGCCGGTGGTTGGCGCGATGCGCTGCGTCAGCGCTTGCTGATGGTCGAAGCCAATGGCCGGCAGGCAGCGCTGGTGGCACTGGTGCTTGGCGACGGTGCGGGCCTGGCCCGTGa                                    <  2:10630/149‑1 (MQ=255)
                                                                                                   ggTACGGTCAAGGCGGGGCAGTTGCTGGCGCCAGTCGCCGGTGGTTGGCGCGATGCGCTGCGTCAGCGCTTGCTGATGGTCGAAGCCAATGGCCGGCAGGCAGCGCTGGTGGCACTGGTGCTTGGCGACGGTGCGGGCCTGGCCCGTGa                                    <  3:52222/149‑1 (MQ=255)
                                                                                                         gTCAAGGCGGGGCAGTTGCTGGCGCCAGTCGCAGGTGGTTGGCGCGATGCGCTGCGTCAGCGCTTGCTGATGGTCGAAGCCAATGGCCGTCAGGCAGCGCTGGTGGCACTGGTGCTTGGCGACGGTGCGGGCCTGGCCCGTGATTACTg                              <  5:173781/149‑1 (MQ=255)
                                                                                                          tCAAGGCGGGGCAGTTGCTGGCGCCAGTCGCCGGTGGTTGGCGCGATGCGCTGCGTCAGCGCTTGCTGATGGTCGAAGCCAATGGCCGGCAGGCAGCGCTGGTGGCACTGGTGCTTGGCGACGGTGCGGGCCTGGCCCGTGAGGACTgg                             >  4:261000/1‑149 (MQ=255)
                                                                                                           cAAGGCGGGGCAGTTGCTGGCGCCAGTCGCCGGTGGTTGGCGCGATGCGCTGCGTCAGCGCTTGCTGATGGTCGAAGCCAATGGCCGGCAGGCAGCGCTGGTGGCACTGGTGCTTGGCGACGGTGCGGGCCTGGCCCGTGAGGACTGGc                            <  6:110311/149‑1 (MQ=255)
                                                                                                            aaGGCGGGGCAGTTGCTGGCGCCAGTCGCCGGTGGTTGGCGCGATGCGCTGCGTCAGGGCTTGCTGATGGTCGAAGCCAATGGCCGGCAGGCAGCGCTGGTGGCACTGGTGCTTGGCGACGGTGCGGGCCTGGCCCGTGAGGACTGGc                            <  2:135586/148‑1 (MQ=255)
                                                                                                            aaGGCGGGGCAGTTGCTGGCGCCAGTCGCCGGTGGTTGGCGCGATGCGCTGCGTCAGCGCTTGCTGATGGTCGAAGCCAATGGCCGGCAGGCAGCGCTGGTGGCACTGGTGCTTGGCGACGGTGCGGGCCTGGCCCGTGAGGACTgg                             >  3:69269/1‑147 (MQ=255)
                                                                                                              ggCGGGGCAGTTGCTGGCGCCAGTCGCCGGTGGTTGGCGCGATGCGCTGCGTCAGCGCTTGCTGATGGTCGAAGCCAATGGCAGGCAGGCAGCGCTGGTGGCACTGGTGCTTGGCGACGGTGCGGGCCTGGCCCGTGAGGACTGGCAGa                         >  4:258932/1‑149 (MQ=255)
                                                                                                               gCGGGGCAGTTGCTGGCGCCAGTCGCCGGTGGTTGGCGCGATGCGCTGCGTCAGCGCTTGCTGATGGTCGAAGCCAATGGCCGGCAGGCAGCGCTGGTGGCACTGGTGCTTGGCGACGGt                                                     >  8:69293/1‑120 (MQ=255)
                                                                                                               gCGGGGCAGTTGCTGGCGCCAGTCGCCGGTGGTTGGCGCGATGCGCTGCGTCAGCGCTTGCTGATGGTCGAAGCCAATGGCCGGCAGGCAGCGCTGGTGGCACTGGTGCTTGGCGACGGt                                                     <  7:69293/120‑1 (MQ=255)
                                                                                                               gCGGGGCAGTTGCTGGCGCCAGTCGCCGGTGGTTGGCGCGATGCGCTGCGTCAGCGCTTGCTGATGGTCGAAGCCAATGGCCGGCAGGCAGCGCTGGTGGCACTGGTGCTTGGCGACGGTGCGGGCCTGGCCCGTGAGGACTGGCAGAc                        >  8:92187/1‑149 (MQ=255)
                                                                                                                    gCAGTTGCTGGCGCCAGTCGCCGGTGGTTGGCGCGATGCGCTGCGTCAGCGCTTGCTGATGGTCGAAGCCAATGGCCGGCAGGCAGCGCTGGTGGCACTGGTGCTTGGCGACGGTGCGGGCCTGGCCCGTGAGGACTGGCAGACGTTg                    >  7:2071/1‑148 (MQ=255)
                                                                                                                       gTTGCTGGCGCCAGTCGCCGGTGGTTGGCGCGATGCGCTGCGTCAGCGCTTGCTGATGGTCGAAGCCAATGGCCGGCAGGCAGCGCTGGTGGCACTGGTGCTTGGCGACGGTGCGGGCCTGGCCCGTGAGGACTGGCAGACGTTGCAgg                >  5:251368/1‑149 (MQ=255)
                                                                                                                        ttGCTGGCGCCAGTCGCCGGTGGTTGGCGCGATGCGCTGCGTCAGCGCTTGCTGATGGTCGAAGCCAATGGCCGGCAGGCAGCGCTGGTGGCACTGGTGCTTGGCGACGGTGCGGGCCTGGCCCGTGAGGACTGGCAGACGTTGCAGGc               <  6:64922/149‑1 (MQ=255)
                                                                                                                        ttGCTGGCGCCAGTCGCCGGTGGTTGGAGCGAAGCGCTGCGTCAGCGCTTGCTGATGGTCGAAGCCAATGGCCGGCACGCAGCGCTGGTGGCACTGGTGCTTGGCGACGGTGCGGGCCTGGCCCGTGAGGACTGGCAGACGTTGCAGGc               >  6:16803/1‑149 (MQ=255)
                                                                                                                          gCTGGCGCCATTCGCCGGTGGTTGGCGCGATGCGCTGCGTCAGCGCTTGCTGATGGTCGAAGCCAATGGCCGGCAGGCAGCGCtggtgg                                                                         <  6:187990/89‑1 (MQ=255)
                                                                                                                          gCTGGCGCCATTCGCCGGTGGTTGGCGCGATGCGCTGCGTCAGCGCTTGCTGATGGTCGAAGCCAATGGCCGGCAGGCAGCGCtggtgg                                                                         >  5:187990/1‑89 (MQ=255)
                                                                                                                           cTGGCGCCAGTCGCCGGTGGTTGGCGCGATGCGCTGCGTCAGCGCTTGCTGATGGTCGAAGCCAATGGCCGGCAGGCAGCGCTGGTGGCACTGGTGCTTGGCGACGGTGCGGGCCTGGCCCGTGAGGACTGGCAGACGTTGCAGGCCAc            >  8:33536/1‑149 (MQ=255)
                                                                                                                                gcCAGTCGCCGGTGGTTGGCGCGATGCGCTGCGTCAGCGCTTGCTGATGGTCGAAGCCAATGGCCGGCAGGCAGCGCTGGTGGCACTGGTGCTTGGCGACGGTGCGGGCCTGGCCCGTGAGGACTGGCAGACGTTGCAGGCCACCGGCa       >  7:183532/1‑149 (MQ=255)
                                                                                                                                    gTCGCCGGTGGTTGGCGCGATGCGCTGCGTCAGCGCTTGCTGATGGTCGAAGCCAATGGCCGGCAGGCAGCGCTGGTGGCACTGGTGCTTGGCGACGGTGCGGGCCTGGCCCGTGAGGACTGGCAGACGTTGCAGGCCACCGGCACGGt   >  5:222877/1‑149 (MQ=255)
                                                                                                                                     tCGCCGGTGGTTGGCGCGATGCGCTGCGTCAGCGCTTGCTGATGGTCGAAGCCAATGGCCGGCAGGCAGCGCTGGTGGCACTGGTGCTTGGCGACGGTGCGGGCCTGGCCCGTGAGGACTGGCAGACGTTGCAGGCCACCGGCACGGTg  <  7:33536/149‑1 (MQ=255)
                                                                                                                                     tCGCCGGTGGTTGGCGCGATGCGCTGCGTCAGCGCTTGCTGATGGTCGAAGCCAATGGCCGGCAGGCAGCGCTGGTGGCACTGGTGCTTGGCGACGGTGCGGGCCTGGCCCGTGAGGACTGGCAGACGTTGCAGGCCACCGGCACGGTg  <  4:246650/149‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                                    
TGGCGGCTGGCCGTGACCTTGCAGCGCCCGAACGGGCTGCTGAACCCGCATGGGCCCGACCGTGAAGCGCAGTTGCTGGCGCGACGGGTGGGGGCTACGGGTACGGTCAAGGCGGGGCAGTTGCTGGCGCCAGTCGCCGGTGGTTGGCGCGATGCGCTGCGTCAGCGCTTGCTGATGGTCGAAGCCAATGGCCGGCAGGCAGCGCTGGTGGCACTGGTGCTTGGCGACGGTGCGGGCCTGGCCCGTGAGGACTGGCAGACGTTGCAGGCCACCGGCACGGTG  >  NC_002947/2140643‑2140924

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