Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,693,792 C→G 6.5% R87R (CGC→CGG)  PP_1489 → CheW domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,693,7920CG 6.5% 198.7 / 3.7 61R87R (CGC→CGG) PP_1489CheW domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (30/27);  new base G (2/2);  total (32/29)
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.91e-01

CGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTG‑GCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCT‑‑GGAGC‑AGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGACACAGGCGATGCGCCTTATCTGGGCCCGGTGCGC  >  NC_002947/1693653‑1693933
                                                                                                                                            |                                                                                                                                                
cgcgTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTG‑GCCAGCCGGCGGTGCAACGCACCAGCACCCGCctggtgctg                                                                                                                                         >  5:7575/1‑149 (MQ=255)
cgcgTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTG‑GCCAGCCGGCGGTGCAACGCACCAGCACCCGCctggtgctg                                                                                                                                         <  2:375493/149‑1 (MQ=255)
 gcgTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTG‑GCCAGCCGGCGGTGCAACGCACCAGCACCCGCctggtgctgg                                                                                                                                        <  2:98379/149‑1 (MQ=255)
  cgTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTG‑GCCAGCCGGCGGTGCAACGCACCAGCACCCGCctggtgctggt                                                                                                                                       <  8:135855/149‑1 (MQ=255)
   gTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTG‑GCCAGCCGGCGGTGCAACGCACCAGCACCCGCctggtgctggtg                                                                                                                                      >  2:78859/1‑149 (MQ=255)
    tCAGCCGCTCAAGCCGAGCGCCCAGGCACCTGCCTGGGGCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGCGCCAGCCGGCGGTGCAACGCACCAGCACCCGCctggtgctggtg                                                                                                                                      <  1:78859/149‑1 (MQ=255)
       gCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTG‑GCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCAc                                                                                                                                   >  6:119223/1‑148 (MQ=255)
       gCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTG‑GCCAGCCGGCGGTGCAACGCACCAGCACCAGGCgggtgctggtgcgtt                                                                                                                                  <  5:111616/149‑4 (MQ=255)
        ccGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTG‑GCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTa                                                                                                                                 <  2:21897/149‑1 (MQ=255)
         cGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTG‑GCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTAc                                                                                                                                <  5:320820/149‑1 (MQ=255)
           cTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTG‑GCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCg                                                                                                                              >  3:241773/1‑149 (MQ=255)
           cTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTG‑GCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCg                                                                                                                              <  6:7575/149‑1 (MQ=255)
              aaGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTG‑GCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGa                                                                                                                             <  8:272778/147‑1 (MQ=255)
              aaGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTG‑GCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAgg                                                                                                                           >  6:21782/1‑149 (MQ=255)
              aaGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTG‑GCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAgg                                                                                                                           >  6:4450/1‑149 (MQ=255)
              aaGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTG‑GCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAg                                                                                                                            <  3:105197/148‑1 (MQ=255)
                  cGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTG‑GCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAggcggc                                                                                                                       >  8:61506/1‑149 (MQ=255)
                  cGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTG‑GCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAggcggc                                                                                                                       <  7:61506/149‑1 (MQ=255)
                   gATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTG‑GCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGcc                                                                                                                      >  3:222837/1‑149 (MQ=255)
                    aTCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTG‑GCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCt                                                                                                                     <  3:73156/149‑1 (MQ=255)
                       gCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTG‑GCCAGCCGGCGGTGCAACGCACCAGCCCCCGCCTGGTGCTGGTGCACTACCGAggcg                                                                                                                         >  1:377321/1‑142 (MQ=255)
                       acccAGGCAGGTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTG‑GCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTg                                                                                                                    >  7:282093/2‑147 (MQ=255)
                        cccAGGCAGGTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTG‑GCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAg                                                                                                                 <  8:282093/149‑1 (MQ=255)
                        cccAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTG‑GCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACt                                                                                                                                  >  8:106366/1‑132 (MQ=255)
                        cccAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTG‑GCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACt                                                                                                                                  <  7:106366/132‑1 (MQ=255)
                         ccAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTG‑GCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGc                                                                                                                <  4:241773/149‑1 (MQ=255)
                               aCCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTG‑GCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGcc                                                                                                          >  6:279415/1‑149 (MQ=255)
                                cctgcctgGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTG‑GCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCt                                                                                                         <  8:112708/149‑1 (MQ=255)
                                    cctgGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGGCATCGACCTTGCCGCCTTGAGCTTTG‑GCCAGCCGGCGGTGCAACGCACCAGCACCAGGCgggtgctggtgcgttgc                                                                                                                                <  7:19397/122‑6 (MQ=38)
                                     ctgGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTG‑GCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATcc                                                                                                    <  4:222837/149‑1 (MQ=255)
                                      tgGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTG‑GCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCATGGCCTGATCCt                                                                                                   >  1:5416/1‑149 (MQ=255)
                                        ggTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTG‑GCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCATGGCCTGATCCT‑‑g                                                                                                <  2:5416/148‑1 (MQ=255)
                                          tCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTG‑GCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAggcggc                                                                                                                       <  4:116956/125‑1 (MQ=255)
                                          tCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTG‑GCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAggcggc                                                                                                                       >  3:116956/1‑125 (MQ=255)
                                              gggCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTG‑GCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCT‑‑GGAGC‑Agg                                                                                        >  4:366778/1‑149 (MQ=255)
                                                        gCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTG‑GCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCT‑‑GGAGC‑AGGCCACCGAAAc                                                                              >  4:55204/1‑149 (MQ=255)
                                                                  gCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTG‑GCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCT‑‑GGAGC‑AGGCCACCGAAACCCTGCGTTGc                                                                    <  3:55204/149‑1 (MQ=255)
                                                                       cTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTG‑GCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCT‑‑GGAGC‑AGGCCACCGAAACCCTGCGTTGCCAAcc                                                               >  7:194009/1‑149 (MQ=255)
                                                                         ggTTCCGGTCATCGGCCTTGCCGCCTTGAGCTTTG‑GCCAGCCGGCGGTGCAACGCACCTGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCT‑‑GGAGC‑AGGCCACCGAAACCCTGCGTTGCCAACCgg                                                             >  7:326145/1‑149 (MQ=255)
                                                                         ggTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTG‑GCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCT‑‑GGAGC‑AGGCCACCGAAACCCTGCGTTGCCAACCg                                                              >  5:149747/1‑148 (MQ=255)
                                                                               ggTCATCGACCTTGCCGCCTTGAGCTTTG‑GCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCt                                                                                                         <  3:311143/102‑1 (MQ=255)
                                                                               ggTCATCGACCTTGCCGCCTTGAGCTTTG‑GCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCt                                                                                                         >  4:311143/1‑102 (MQ=255)
                                                                                        ccTTGCCGCCTTGAGCTTTG‑GCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCT‑‑GGAGC‑AGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCTAGCCCt                                              <  5:21782/149‑1 (MQ=255)
                                                                                        ccTTGCCGCCTTGAGCTTTG‑GCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCT‑‑GGAGC‑AGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCt                                              <  6:149747/149‑1 (MQ=255)
                                                                                         cTTGCCGCCTTGAGCTTTG‑GCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCTTGGGGGCGAGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGcc                                                >  1:126624/1‑149 (MQ=255)
                                                                                          ttGCCGCCTTGAGCTTTG‑GCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCT‑‑GGAGC‑AGGCCACCGAAACCCTGCGTTg                                                                     >  7:165012/1‑124 (MQ=255)
                                                                                          ttGCCGCCTTGAGCTTTG‑GCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCT‑‑GGAGC‑AGGCCACCGAAACCCTGCGTTg                                                                     <  8:165012/124‑1 (MQ=255)
                                                                                              cgccTTGAGCTTTG‑GCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCT‑‑GGAGC‑AGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATgggg                                        >  5:326773/1‑149 (MQ=255)
                                                                                                       cTTTG‑GCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCT‑‑GGAGC‑AGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGACACAg                               >  3:32321/1‑149 (MQ=255)
                                                                                                          tG‑GCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCT‑‑GGAGC‑AGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGACACAgg                              <  6:326773/147‑1 (MQ=255)
                                                                                                             gCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCT‑‑GGAGC‑AGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGACACAGGCGAt                          >  8:192673/1‑149 (MQ=255)
                                                                                                              ccAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCT‑‑GGAGC‑AGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGACACAGGCGATg                         <  3:366778/149‑1 (MQ=255)
                                                                                                                 gCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCT‑‑GGAGC‑AGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGACACAGGCGATgcgc                      <  5:279415/149‑1 (MQ=255)
                                                                                                                           cAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCT‑‑GGAGC‑AGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGACACAGGCGATGCGCCTTATCTgg             <  5:119223/148‑1 (MQ=255)
                                                                                                                            aaCGCACCAGCACCCGCCTGGTGCTGGTGCACTCCCGACGCGGCCTGCAGCTTCGCCTGATCCT‑‑CGAGC‑AGGCCCCCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGCTGGACACAGGCGATGCGACTTATCTGGGc           >  1:219864/1‑149 (MQ=255)
                                                                                                                            aaCGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCT‑‑GGAGC‑AGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGACACAGGCGATGCGCCTTAt                 >  3:204653/1‑143 (MQ=255)
                                                                                                                             aCGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCT‑‑GGAGC‑AGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGACACAGGCGATGCGCCTTATCTGGGcc          >  3:153149/1‑149 (MQ=255)
                                                                                                                                caccagcaccagGCGGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCT‑‑GGAGC‑AGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGACACAGGCGATGCGCCTTATCTGGGCCCgg       >  7:313433/1‑149 (MQ=255)
                                                                                                                                  ccAGCACCAGGCGGGTGCTGGTGCGCTACCGAGGCGGCCTGCAGCTTGGCCTGATCCT‑‑GGAGC‑AGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGACACAGGCGATGCGCCTTATCTGGGCCCGGTg     >  6:103648/1‑149 (MQ=255)
                                                                                                                                   cAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCT‑‑GGAGC‑AGGCCACCGAAACCCTGCGTTGCCACCCGGACGAGTTCCAGCCCTATGGGGTGGACACAGGCGATGCGCCTTATCTGGGCCCGGTgc    >  5:187408/1‑149 (MQ=255)
                                                                                                                                     gCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCT‑‑GGAGC‑AGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGACACAGGCGATGCGCCTTATCTGGGCCCGGTgcgc  >  5:134409/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                                
CGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTG‑GCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTGATCCT‑‑GGAGC‑AGGCCACCGAAACCCTGCGTTGCCAACCGGACGAGTTCCAGCCCTATGGGGTGGACACAGGCGATGCGCCTTATCTGGGCCCGGTGCGC  >  NC_002947/1693653‑1693933

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: ‑