Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 859,283 A→T 5.3% L144Q (CTG→CAG)  PP_0740 ← MerR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947859,2830AT 5.3% 222.3 / 2.3 76L144Q (CTG→CAG) PP_0740MerR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (36/35);  new base T (2/2);  total (39/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.15e-01

GCCAGGCGCACAACCATAAGTCCGGGCTGAAGGTGCGTTCGCTGTCTTCGGCCAGCAACACGGCAGGCCCTTGCAGCAGCTGGTTGTCGTGGTACACGCGGGCACCGAGCTTGCTGCGCAGCCAGCTGTGAAAAAACACCTGCTCCAGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGTCGAGCATGTCGAGGAGCGGCAGCAGCAATTGCTCGCACAGGGTGACCGCCGGGTACAGCGCCATGGCCTGGTTGAGTTGCTGGTCGAGGGCGCGCTGTGAAAGGCCGGCGATG  >  NC_002947/859137‑859429
                                                                                                                                                  |                                                                                                                                                   
gCCAGGCGCACAACCATAAGTCCGGGCTGAAGGTGCGTTCGCTGTCTTCGGCCAGCAACACGGCAGGCCCTTGCAGCAGCTGGTTGTCGTGGTACACGCGGGCACCGAGCTTGCTGCGCAGCCAGCTGTGAAAAAACACCTGCTCCAGc                                                                                                                                                   <  3:284125/149‑1 (MQ=255)
  cAGGCGCACAACCATAAGTCCGGGCTGAAGGTGCGTTCGCTGTCTTCGGCCAGCAACACGGCAGGCCCTTGCAGCAGCTGGTTGTCGTGGTACACGCGGGCACCGAGCTTGCTGCGCAGCCAGCTGTGAAAAAACACCTGCTCCAGcc                                                                                                                                                  >  4:339896/1‑148 (MQ=255)
     gcgcACAACCATAAGTCCGGGCTGAAGGTGCGTTCGCTGTCTTCGGCCAGCAACACGGCAGGCCCTTGCAGCAGCTGGTTGTCGTGGTACACGCGGGCACCGAGCTTGCTGCGCAGCCAGCTGTGAAAAAACACCTGCTCCAGCCGAGc                                                                                                                                              >  6:135156/1‑149 (MQ=255)
        cacaACCATAAGTCCGGGCTGAAGGTGCGTTCGCTGTCTTCGGCCAGCAACACGGCAGGCCCTTGCAGCAGCTGGTTGTCGTGGTACACGCGGGCACCGAGCTTGCTGCGCAGCCAGCTGTGAAAAAACACCTGCTCCAGCCGAGCGtt                                                                                                                                           <  1:98637/149‑1 (MQ=255)
        cacaACCATAAGTCCGGGCTGAAGGTGCGTTCGCTGTCTTCGGCCAGCAACACGGCAGGCCCTTGCAGCAGCTGGTTGTCGTGGTACACGCGGGCACCGAGCTTGCTGCACAGCCAGCTGTGAAAAAACACCTGCTCCAGCCGAGCGtt                                                                                                                                           >  3:286118/1‑149 (MQ=255)
           aaCCATAAGTCCGGGCTGAAGGTGCGTTCGCTGTCTTCGGCCAGCAACACGGCAGGCCCTTGCAGCAGCTGGTTGTCGTGGTACACGCGGGCACCGAGCTTGCTGCGCAGCCAGCTGTGAAAAAACACCTGCTCCAGCCGAGCGTTGaa                                                                                                                                        <  1:245291/149‑1 (MQ=255)
           aaCCATAAGTCCGGGCTGAAGGTGCGTTCGCTGTCTTCGGCCAGCAACACGGCAGGCCCTTGCAGCAGCTGGTTGTCGTGGTACACGCGGGCACCGAGCTTGCTGCGCAGCCAGCTGTGAAAAAACACCTGCTCCAGCCGAGCGTTGaa                                                                                                                                        <  5:181390/149‑1 (MQ=255)
             ccATAAGTCCGGGCTGAAGGTGCGTTCGCTGTCTTCGGCCAGCAACACGGCAGGCCCTTGCAGCAGCTGGTTGTCGTGGTACACGCGGGCACCGAGCTTGCTGCGCAGCCAGCTGTGAAAAAACACCTGCTCCAGCCGAGCGTTGAAg                                                                                                                                       >  3:303781/1‑148 (MQ=255)
                         gCTGAACGTCCGATCGCTGTCTTCTGCCAGCAACACGGCAGGCCCTTGCACCAGCTGGTTGTCGTCGAACACGCTCGCACCGAGCTTTCTCCGCATCCACCTGTGAAAAAACACCTGCTCCAGCCGAGCGGTGAACTAATTGCGCCAgc                                                                                                                          >  6:25390/1‑149 (MQ=255)
                           tGAAGGTGCGTTCGCTGTCTTCGGCCAGCAACACGGCAGGCCCTTGCAGCAGCTGGTTGTCGTGGTACACGCGGGCACCGAGCTTGCTGCGCAGCCAGCTGTGAAAAAACACCTGCTCCAGCCGAGCGTTGAAGTAATTGCGCCAgcg                                                                                                                         >  4:86528/1‑148 (MQ=255)
                               ggTGCGTTCGCTGTCTTCGGCCAGCAACACGGCAGGCCCTTGCAGCAGCTGGTTGTCGTGGTACACGCGGGCACCGAGCTTGCTGCGCAGCCAGCTGTGAAAAAACACCTGCTCCAGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGt                                                                                                                    <  5:135156/149‑1 (MQ=255)
                                gTGCGTTCGCTGTCTTCGGCCAGCAACACGGCAGGCCCTTGCAGCAGCTGGTTGTCGTGGTACACGCGGGCACCGAGCTTGCTGCGCAGCCAGCTGTGAAAAAACACCTGCTCCAGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGTc                                                                                                                   >  1:130011/1‑149 (MQ=255)
                                   cGTTCGCTGTCTTCGGCCAGCAACACGGCAGGCCCTTGCAGCAGCTGGTTGTCGTGGTACACGCGGGCACCGAGCTTGCTGCGCAGCCAGCTGTGAAAAAACACCTGCTCCAGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGTCGAg                                                                                                                >  6:58413/1‑149 (MQ=255)
                                    gTTCGCTGTCTTCTGCCAGCAACACGGCATGCCCTTGCAGCAGCTGGTTGTCGTGGTACACGCGGGCACCGAGCTTGCTGCGCAGCCAGCTGTGAAAAAACACCTGCTCCATCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGTCGAGc                                                                                                               <  5:269415/149‑1 (MQ=255)
                                        gCTGTCTTCGGCCAGCAACACGGCAGGCCCTTGCAGCAGCTGGTTGTCGTGGTACACGCGGGCACCGAGCTTGCTGCGCAGCCAGCTGTGAAAAAACACCTGCTCCAGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGTCGAGCATGt                                                                                                           <  2:134377/149‑1 (MQ=255)
                                          tGTCTTCGGCCAGCAACACGGCAGGCCCTTGCAGCAGCTGGTTGTCGTGGTACACGCGGGCACCGAGCTTGCTGCGCAGCCAGCTGTGAAAAAACACCTGCTCCAGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGTCGAg                                                                                                                <  8:225952/142‑1 (MQ=255)
                                          tGTCTTCGGCCAGCAACACGGCAGGCCCTTGCAGCAGCTGGTTGTCGTGGTACACGCGGGCACCGAGCTTGCTGCGCAGCCAGCTGTGAAAAAACACCTGCTCCAGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGTCGAg                                                                                                                >  7:225952/1‑142 (MQ=255)
                                             cTTCGGCCAGCAACACGGCAGGCCCTTGCAGCAGCTGGTTGTCGTGGTACACGCGGGCACCGAGCTTGCTGCGCAGCCAGCTGTGAAAAAACACCTGCTCCAGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGTCGAGCATGTCgagg                                                                                                      <  3:36189/149‑1 (MQ=255)
                                             cTTCGGCCAGCAACACGGCAGGCCCTTGCAGCAGCTGGTTGTCGTGGTACACGCGGGCACCGAGCTTGCTGCGCAGCCAGCTGTGAAAAAACACCTGCTCCAGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGTCGAGCATGTCgag                                                                                                       >  4:371812/1‑148 (MQ=255)
                                              ttCGGCCAGCAACACGGCAGGCCCTTGCAGCAGCTGGTTGTCGTGGTACACGCGGGCACCGAGCTTGCTGCGCAGCCAGCTGTGAAAAAACACCTGCTCCAGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGTCGAGCATGTCgag                                                                                                       >  3:140372/1‑147 (MQ=255)
                                                  gCCAGCAACACGGCAGGCCCTTGCAGCAGCTGGTTGTCGTGGTACACGCGGGCACCGAGCTTGCTGCGCAGCCAGCTGTGAAAAAACACCTGCTCCAGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGTCGAGCATGTCGAGGAGCgg                                                                                                 <  2:360628/149‑1 (MQ=255)
                                                   ccAGCAACACGGCAGGCCCTTGCAGCAGCTGGTTGTCGTGGTACACGCGGGCACCGAGCTTGCTGCGCAGCCAGCTGTGAAAAAACACCTGCTCCAGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGTCGAGCATGTCGAGGAGCGgc                                                                                                <  8:270753/149‑1 (MQ=255)
                                                     aGCAACACGGCAGGCCCTTGCAGCAGCTGGTTGTCGTGGTACACGCGGGCACCGAGCTTGCTGCGCAGCCAGCTGTGAAACAACACCTGCTCCAGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGTCGAGCATGTCGAGGAGCGgcag                                                                                              >  8:219833/1‑149 (MQ=255)
                                                       cAACACGGCAGGCCCTTGCAGCAGCTGGTTGTCGTGGTACACGCGGGCACCGAGCTTGCTGCGCAGCCAGCTGTGAAAAAACACCTGCTCCAGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGTCGAGCATGTCGAGGAGCGgcagca                                                                                            <  5:104216/149‑1 (MQ=255)
                                                         acacGGCAGGCCCTTGCAGCAGCTGGTTGTCGTGGTACACGCGGGCACCGAGCTTGCTGCGCAGCCAGCTGTGAAAAAACACCTGCTCCAGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGTCGAGCATGTCGAGGAGCGgcagcag                                                                                           >  1:59918/1‑148 (MQ=255)
                                                          cacGGCAGGCCCTTGCAGCAGCTGGTTGTCGTGGTACACGCGGGCACCGAGCTTGCTGCGCAGCCAGCTGTGAAAAAACACCTGCTCCAGCCGAGCGTTGAAGTAATTGAGCCAGCGCAGGTCGAGCATGTCGAGGAGCGgcagcagc                                                                                          >  6:287786/1‑148 (MQ=255)
                                                              gCAGGCCCTTGCAGCAGCTGGTTGTCGTGGTACACGCGGGCACCGAGCTTGCTGCGCAGCCAGCTGTGAAAAAACACCTGCTCCAGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGTCGAGCATGTCGAGGAGCGGCAGCAGCAATTg                                                                                    >  2:101958/1‑149 (MQ=255)
                                                                    ccTTGCAGCAGCTGGTTGTCGTGGTACACGCGGGCACCGAGCTTGCTGCGCAGCCAGCTGTGAAAAAACACCTGCTCCTGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGTCGAGCATGTCgag                                                                                                       <  7:261992/125‑1 (MQ=255)
                                                                    ccTTGCAGCAGCTGGTTGTCGTGGTACACGCGGGCACCGAGCTTGCTGCGCAGCCAGCTGTGAAAAAACACCTGCTCCTGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGTCGAGCATGTCgag                                                                                                       >  8:261992/1‑125 (MQ=255)
                                                                         cagcagCTGGTTGTCGTGGTACACGCGGGCACCGAGCTTGCTGCGCAGCCAGCTGTGAAAAAACACCTGCTCCAGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGTCGAGCATGTCGAGGAGCGGCAGCAGCAATTGCTCGCACAgg                                                                          <  8:530/148‑1 (MQ=255)
                                                                         cagcagCTGGTTGTCGTGGTACACGCGGGCACCGAGCTTGCTGCGCAGCCAGCTGTGAAAAAACACCTGCTCCAGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGTCGAGCATGTCGAGGAGCGGCAGCAGCAATTGCTCGCACAgg                                                                          <  5:25390/148‑1 (MQ=255)
                                                                          agcagcTGGTTGTCGTGGTACACGCGGGCACCGAGCTTGCTGCGCAGCCAGCTGTGAAAAAACACCTGCTCCAGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGTCGAGCATGTCGAGGAGCGGCAGCAGCAATTGCTCGCACAGGGt                                                                        <  5:58413/149‑1 (MQ=255)
                                                                             agcTGGTTGTCGTGGTACACGCGGGCACCGAGCTTGCTGCTCAGCCAGCTGTGAAAAAACACCTGCTCCAGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGTCGAGCATGTCGAGGAGCGGCAGCAGCAATTGCTCGCACAGGGTGAc                                                                     <  4:286118/149‑1 (MQ=255)
                                                                              gcTGGTTGTCGTGGTACACGCGGGCACCGAGCTTGCTGCGCAGCCAGCTGTGAAAAAACACCTGCTCCAGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGTCGAGCATGTCGAGGAGCGGCAGCAGCAATTGCTCGCACAGGGTGAcc                                                                    >  7:187238/1‑149 (MQ=255)
                                                                                  gTTGTCGTGGTACACGCGGGCACCGAGCTTGCTGCGCAGCCAGCTGTGAAAAAACACCTGCTCCAGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGTCGAGCATGTCGAGGAGCGGCAGCAGCAATTGCTCGCACAGGGTGAccgccg                                                                >  3:47904/1‑149 (MQ=255)
                                                                                   ttGTCGTGGTACACGCGGGCACCGAGCTTGCTGCGCAGCCAGCTGTGAAAAAACACCTGCTCCAGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGTCGAGCATGTCGAGGAGCGGCAGCAGCAATTGCTCGCACAGGGTGACCGCCgg                                                               >  8:78092/1‑149 (MQ=255)
                                                                                     gTCGTGGTACACGCGGGCACCGAGCTTGCTGCGCAGCCAGCTGTGAAAAAACACCTGCTCCAGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGTCGAGCATGTCGAGGAGCGGCAGCAGCAATTGCTCGCACAGGGTGACCGCCggg                                                              <  4:338996/148‑1 (MQ=255)
                                                                                     gTCGTGGTACACGCGGGCACCGAGCTTGCTGCGCAGCCAGCTGTGAAAAAACACCTGCTCCAGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGTCGAGCATGTCGAGGAGCGGCAGCAGCAATTGCTCGCACAGGGTGACCGCCggg                                                              <  3:371812/148‑1 (MQ=255)
                                                                                      tCGTGGTACACGCGGGCACCGAGCTTGCTGCGCAGCCAGCTGTGAAAAAACACCTGCTCCAGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGTCGAGCATGTCGAGGGGCGGCAGCAGCAATTGCTCGCACAGGGTGACCGCCGGGTa                                                            >  4:138806/1‑149 (MQ=255)
                                                                                        gTGGTACACGCGGGCACCGAGCTTGCTGCGCAGCCAGCTGTGAAAAAACACCTGCTCCAGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGTCGAGCATGTCGAGGAGCGGCAGCAGCAATTGCTCGCACAGGGTGACCGCCGGGTAc                                                           <  2:52033/148‑1 (MQ=255)
                                                                                          ggTACACGCGGGCACCGAGCTTGCTGCGCAGCCAGCTGTGAAAAAACACCTGCTCCAGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGTCGAGCATGTCGAGGAGCGGCAGCAGCAATTGCTCGCACAGGGTGACCGCCGGGTACAgc                                                        >  6:282338/1‑149 (MQ=255)
                                                                                           gTACACGCGGGCACCGAGCTTGCTGCGCAGCCAGCTGTGAAAAAACACCTGCTCCAGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGCTCGAGCATGTCGAGGAGCGGCAGCAGCAATTGCTCGCACAGGGTGACCGCCGGGTACAgcg                                                       >  3:182082/1‑149 (MQ=255)
                                                                                                cgcgGGCACCGAGCTTGCTGCGCAGCCAGCTGTGAAAAAACACCTGCTCCAGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGTCGAGCATGTCGAGGAGCGGCAGCAGCAATTGCTCGCACAGGGTGACCGCCGGGTACAGCGCCATg                                                  >  3:283776/1‑149 (MQ=255)
                                                                                                       aCCGAGCTTGCTGCGCAGCCAGCTGTGAAAACACACCTGCTCCCGCCGAGCGTGGAAGTAATTGCGCCAGCGCCGGTCGAGCAGGTCGAGGCGCGGCAGCCGCAATTGCTCGCCCAGGGTGCCCGCCGGGGACAGCGCCATGGCCAGGt                                           >  3:11250/1‑149 (MQ=255)
                                                                                                       aCCGAGCTTGCTGCGCAGCCAGCTGTGAAAAAACACCTGCTCCAGCCGAGCGTTGaa                                                                                                                                        >  4:113581/1‑57 (MQ=255)
                                                                                                        ccGAGCTTGCTGCGCAGCCAGCTGTGAAAAAACACCTGCTCCAGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGTCGAGCATGTCGAGGAGCGGCAGCAGCAATTGCTCGCACAGGGTGACCGCCGGGTACAGCGCCATGGCCTGGtt                                          <  3:85027/149‑1 (MQ=255)
                                                                                                         cGAGCTTGCTGCGCAGCCAGCTGTGAAAAAACACCTGCTCCTGCCGAGCGTTGAAGTAATTGCGTCAGCGCAGGTCGAGCATGTCGAGGAGCGGCAGCAGCAATTGCTCGCACAGGGTGACCGCCGGGTACAgc                                                        <  2:55373/134‑1 (MQ=255)
                                                                                                         cGAGCTTGCTGCGCAGCCAGCTGTGAAAAAACACCTGCTCCTGCCGAGCGTTGAAGTAATTGCGTCAGCGCAGGTCGAGCATGTCGAGGAGCGGCAGCAGCAATTGCTCGCACAGGGTGACCGCCGGGTACAgc                                                        >  1:55373/1‑134 (MQ=255)
                                                                                                          gAGCTTGCTGCGCAGCCAGCTGTGAAAAAACACCTGCTCCAGCCGAGCGTTGAAGTACTTGCCCCAGCGCAGGTCGAGCATGTCGAGGAGCGTCAGCAGCAATTGCTCTCACAGGGTGACCGCCTGGTACAGCGCCATGGCCTGGTTGa                                        <  8:1605/149‑1 (MQ=255)
                                                                                                                gctgcGCAGCCAGCTGTGAAAAAACACCTGCTCCAGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGTCGAGCATGTCGAGGAGCGGCAGCAGCAATTGCTCGCACAGGGTGACCGCCGGGTACAGCGCCATGGCCTGGTTGAGTTGCt                                  >  5:118459/1‑149 (MQ=255)
                                                                                                                 ctgcGCAGCCAGCTGTGAAAAAACACCTGCTCCAGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGTCGAGCATGTCGAGGAGCGGCAGCAGCAACCTGCTCGCACAGGGTGACCGCCGGGTACAGCGCCATGGCCTGGTTGAGTTGCt                                  <  1:101958/149‑1 (MQ=255)
                                                                                                                 ctgcGCAGCCAGCTGTGAAAAAACACCTGCTCCAGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGTCGAGCATGTCGAGGAGCGGCAGCAGCAATTGCTCGCACAGGGTGACCGCCGGGTACAGCGCCATGGCCTGGTTGAGTTGCt                                  <  5:32510/148‑1 (MQ=255)
                                                                                                                      cagccagcTGTGAAAAAACACCTGCTCCAGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGTCGAGCATGTCGAGGAGCGGCAGCAGCAATTGCTCGCACAGGGTGACCGCCGGGTACAGCGCCATGGCCTGGTTGAGTTGCTGGTc                              >  1:231861/1‑147 (MQ=255)
                                                                                                                        gccagcTGTGAAAAAACACCTGCTCCAGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGTCGAGCATGTCGAGGAGCGGCAGCAGCAATTGCTCGCACAGGGTGACCGCCGGGTACAGCGCCATGGCCTGGTTGAGTTGCTGGTCGAgg                          <  4:303781/149‑1 (MQ=255)
                                                                                                                          cagcTGTGAAAAACCCCCTGCGCCAGCCGAGCGTTGAAGTAATAGCGCCAGCGCAGGTCGAGCATCTCGAGGAGCGGCAGCAGCACTTGCGCGCACAGGGGGACCGCCGGGTACAGCGCCCTGGCCTGGTTGAGTTGCTGGTCGAGGgc                        <  6:118459/149‑1 (MQ=255)
                                                                                                                          cagcTGTGAAAAAACACCTGCTCCAGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGTCGAGCATGTCGAGGAGCGGCAGCAGCAATTGCTCGCACAGGGTGACCGCCGGGTACAGCGCCATGGCCTGGTTGAGTTGCTGGTCGAGGgc                        <  8:187238/149‑1 (MQ=255)
                                                                                                                           agcTGTGAAAAAACACCTGCTCCAGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGTCGAGCATGTCGAGGAGCGGCAGCAGCAATTGCTCGCACAGGGTGACCGCCGGGTACAGCGCCATGGCCTGGTTGAGTTGCTGGTCGAGGgcg                       <  4:140372/149‑1 (MQ=255)
                                                                                                                             cTGTGAAAAAACACCTGCTCCAGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGTCGAGCATGTCGAGGAGCGGCAGCAGCAATTGCTCGCACAGGGTGACCGCCGGGTACAGCGCCATGGCCTGGTTGAGTTGCTGGTCGAGGgcgcg                     >  3:127821/1‑149 (MQ=255)
                                                                                                                              tgtgAAAAAACACCTGCTCCAGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGTCGAGCATGTCGAGGAGCGGCAGCAGCAATTGCTCGCACAGGGTGACCGCCGGGTACAGCGCCATGGCCTGGTTGAGTTGCTGGTCGAGGgcgcgc                    >  1:356663/1‑149 (MQ=255)
                                                                                                                              tgtgAAAAAACACCTGCTCCAGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGTCGAGCATGTCGAGGAGCGGCAGCAGCAATTGCTCGCACAGGGTGACCGCCGGGTACAGCGCCATGGCCTGGTTGAGTTGCTGGTCGAGGgcgcgc                    >  8:317243/1‑149 (MQ=255)
                                                                                                                                  aaaaaaCACCTGCTCCAGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGTCGAGCATGTCGAGGAGCGGCAGCATCAATTGCTAGCACAGGGTGACCGCCTGGTACAGCGCCATGGCCTTGGTGAGTTGCTGGTCGAGGGCGCGCtgg                 >  4:252908/1‑147 (MQ=255)
                                                                                                                                  aaaaaaCACCTGCTCCAGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGTCGAGCATGTCGAGGAGCGGCAGCAGCAATTGCTCGCACAGGGTGACCGCCGGGTACAGCGCCATGGCCTGGTTGAGTTGCTGGTCGAGGGCGCGCtgt                 <  8:122395/148‑1 (MQ=255)
                                                                                                                                   aaaaaCACCTGCTCCAGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGTCGAGCATGTCGAGGAGCGGCAGCAGCAATTGCTCGCACAGGGTGACCGCCGGGTACAGCGCCATGGCCTGGTTGAGTTGCTGGTCGAGGGCGCGCtgtg                >  4:48074/1‑148 (MQ=255)
                                                                                                                                    aaaaCACCTGCTCCAGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGTCGAGCATGTCGAGGAGCGGCAGCAGCAATTGCTCGCACAGGGTGACCACCGGGTACAGCGCCATGGCCTGGTTGAGTTGCTGGTCGAGGGCGCGCTGTGaa              <  3:335186/149‑1 (MQ=255)
                                                                                                                                         acCTGCTCCAGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGTCGAGCATGTCGAGGAGCGGCAGCAGCAATTGCTCGCACAGGGTGACCGCCGGGTACAGCGCCATGGCCTGGTTGAGTTGCTGGTCGAGGGCGCGCTGTGAAAGGcc         >  7:276561/1‑149 (MQ=255)
                                                                                                                                          ccTGCTCCAGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGTCGAGCATGTCGAGGAGCGGCAGCAGCAATTGCTCGCACAGGGTGACCGCCGGGTACAGCGCCATGGCCTGGTTGAGTTGCTGGTCGAGGGCGCGCTGTGAAAGGCCg        <  3:86528/149‑1 (MQ=255)
                                                                                                                                           cTGCTCCAGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGTCGAGCATGTCGAGGAGCGGCAGCAGCAATTGCTCGCACAGGGTGACCGCCGGGTACAGCGCCATGGCCTGGTTGAGTTGCTGGTCGAGGGCGCGCTGTGAAAGGCCgg       >  2:47902/1‑149 (MQ=255)
                                                                                                                                           cTGCTCCAGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGTCGAGCATGTCGAGGAGCGGCAGCAGCAATTGCTCGCACAGGGTGACCGCCGGGTACAGCGCCATGGCCTGGTTGAGTTGCTGGTCGAGGGCGCGCTGTGAAAGGCCgg       >  4:389120/1‑149 (MQ=255)
                                                                                                                                           cTGCTCCAGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGTCGAGCATGTCGAGGAGCGGCAGCAGCAATTGCTCGCACAGGGTGACCGCCGGGTACAGCGCCATGGCCTGGTTGAGTTGCTGGTCGAGGGCGCGCTGTGAAAGGCCgg       <  5:287786/149‑1 (MQ=255)
                                                                                                                                               tCCAGGCGAGCGTTGAAGTAATTGCGCGAGCTCAGTTCGAGAATGTCGAGGATCGGCAGGAGTAATTGCTCGCACAGTGTGACCGACGGGTACAGCGGCATGGCCTGGTTGAGTTGCTGGTCGAGGTCGAGCTGTGAAAGGGCGGCGAt   <  4:182082/149‑1 (MQ=255)
                                                                                                                                               tCCAGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGTCGAGCATGTCGAGGAGCGGCAGCAGCAATTGCTCGCACAGGGTGACCGCCGGGTACAGCGCCATGGCCTGGTTGAGTTGCTGGTCGAGGGCGCGCTGTGAAAGGCCGGCGAt   >  4:365264/1‑149 (MQ=255)
                                                                                                                                                ccAGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGTCGAGCATGTCGAGGAGCGGCAGCAGCAATTGCTCGCACAGGGTGACCGCCGGGTACAGCGCCATGGCCTGGTTGAGTTGCTGGTCGAGGGCGCGCTGTGAAAGGCCGGCGATg  <  7:78092/149‑1 (MQ=255)
                                                                                                                                                 cAGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGTCGAGCATGTCGAGGAGCGGCAGCAGCAATTGCTCGCACAGGGTGACCGCCGGGTACAGCGCCATGGCCTGGTTGAGTTGCTGGTCGAGGGCGCGCtgtg                >  6:49603/1‑134 (MQ=255)
                                                                                                                                                 cAGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGTCGAGCATGTCGAGGAGCGGCAGCAGCAATTGCTCGCACAGGGTGACCGCCGGGTACAGCGCCATGGCCTGGTTGAGTTGCTGGTCGAGGGCGCGCTGTGaaa             >  2:384314/1‑137 (MQ=255)
                                                                                                                                                 cAGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGTCGAGCATGTCGAGGAGCGGCAGCAGCAATTGCTCGCACAGGGTGACCGCCGGGTACAGCGCCATGGCCTGGTTGAGTTGCTGGTCGAGGGCGCGCTGTGaaa             <  1:384314/137‑1 (MQ=255)
                                                                                                                                                 cAGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGGCGAGCATGTCGAGGAGCGGCAGCAGCAATTGCTCGCACAGGGTGACCGCCGGGTACAGCGCCATGGCCTGGTTGAGTTGCTGGTCGAGGGCGCGCtgtg                <  5:49603/134‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GCCAGGCGCACAACCATAAGTCCGGGCTGAAGGTGCGTTCGCTGTCTTCGGCCAGCAACACGGCAGGCCCTTGCAGCAGCTGGTTGTCGTGGTACACGCGGGCACCGAGCTTGCTGCGCAGCCAGCTGTGAAAAAACACCTGCTCCAGCCGAGCGTTGAAGTAATTGCGCCAGCGCAGGTCGAGCATGTCGAGGAGCGGCAGCAGCAATTGCTCGCACAGGGTGACCGCCGGGTACAGCGCCATGGCCTGGTTGAGTTGCTGGTCGAGGGCGCGCTGTGAAAGGCCGGCGATG  >  NC_002947/859137‑859429

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: