Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,253,935 G→C 5.2% F249L (TTC→TTG PP_4630 ← MerR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,253,9350GC 5.2% 245.8 / 3.0 77F249L (TTC→TTGPP_4630MerR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (39/33);  new base C (2/2);  total (42/35)
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.64e-01

CAGAAACTGCCGCAAGCGCTGGCGCATCACACTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTTCGTCGCTGCCCA  >  NC_002947/5253791‑5254077
                                                                                                                                                |                                                                                                                                              
cAGAAACTGCCGCAAGCGCTGGCGCATCACACTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTg                                                                                                                                           >  1:44758/1‑150 (MQ=255)
cAGAAACTGCCGCAAGCGCTGGCGCATCACACTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCACTTCAAgcg                                                                                                                                           >  1:470233/1‑148 (MQ=255)
    aaCTGCCGCAAGCGCTGGCGCATCACACTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCACTTCAAGCGCCTg                                                                                                                                       <  2:470233/150‑1 (MQ=255)
     aCTGCCGCAAGCGCTGGCGCATCACACTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACAAACCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTgg                                                                                                                                      >  1:1020773/1‑150 (MQ=255)
     aCTGCCGCAAGCGCTGGCGCATCACACTGCCCTGCGCGCCAAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTgg                                                                                                                                      >  1:828803/1‑150 (MQ=255)
     aCTGCCGCAAGCGCTGGCGCATCACACTGCACTAAACCAACAAGCCCCACTCTCGGATTCCACCGGGGCTTTCCTGGACCGCTTCAGCTGCTTCACCAGCCAGTAAGAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTgg                                                                                                                                      >  1:739028/1‑150 (MQ=255)
      cTGCCGCAAGCGCTGGCGCATCACACTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAGGTGCCTGGc                                                                                                                                     <  2:567102/150‑1 (MQ=255)
            cAAGCGCTGGCGCATCACACTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTc                                                                                                                               <  1:995445/150‑1 (MQ=255)
             aaGCGCTGGCGCATCACACTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTct                                                                                                                              >  1:297351/1‑150 (MQ=255)
                cgcTGGCGCATCACACTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTctgct                                                                                                                           <  1:1145323/150‑1 (MQ=255)
                  cTGGCGCATCACACTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCtgt                                                                                                                         <  2:1079825/150‑1 (MQ=255)
                    ggCGCATCACACTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTgg                                                                                                                       >  2:1164642/1‑150 (MQ=255)
                    ggCGCATCACACTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTgg                                                                                                                       <  2:844618/150‑1 (MQ=255)
                      cgcATCACACTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGc                                                                                                                     <  2:1304928/150‑1 (MQ=255)
                       gcATCACACTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCg                                                                                                                    >  1:643590/1‑150 (MQ=255)
                        cATCACACTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTgg                                                                                                                       <  1:1362285/146‑1 (MQ=255)
                        cATCACACTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTgg                                                                                                                       >  2:1362285/1‑146 (MQ=255)
                         aTCACACTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGAc                                                                                                                  >  1:902208/1‑150 (MQ=255)
                          tCACACTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACt                                                                                                                 >  1:197382/1‑150 (MQ=255)
                              acTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATg                                                                                                             <  1:1397241/150‑1 (MQ=255)
                                tGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATggg                                                                                                           >  1:670930/1‑150 (MQ=255)
                                    cTGCGCGCCCAGGCAGGCTACAGGGCTGCCAACCAGGCAGATTTGAAACACTTCCGTTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGACTGGAGAGTTCTGCTTTGGTCGACTGATTTGaaaa                                                                                                       <  2:1020773/150‑1 (MQ=255)
                                    cTGCGCGCCAAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGaaaa                                                                                                       >  2:280394/1‑150 (MQ=255)
                                          gcCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAAAAACGGGCATTCCAGCATCGACTCCAGGCGCGTTAAGGGCCTGGCTAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCAc                                                                                                 >  2:871211/1‑150 (MQ=255)
                                            ccAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCAccc                                                                                               >  2:417615/1‑149 (MQ=255)
                                            ccAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCa                                                                                               >  1:476775/1‑150 (MQ=255)
                                              aggcaggcTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAgg                                                                                             <  2:670930/150‑1 (MQ=255)
                                              aggcaggcTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAgg                                                                                             <  2:643590/150‑1 (MQ=255)
                                              aggcaggcTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAgg                                                                                             <  2:828803/150‑1 (MQ=255)
                                              aggcaggcTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAgg                                                                                             <  2:902208/150‑1 (MQ=255)
                                                    gcTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTcc                                                                                       >  2:609097/1‑150 (MQ=255)
                                                      tACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGc                                                                                                                                          <  2:792229/97‑1 (MQ=255)
                                                      tACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGc                                                                                                                                          <  2:792818/97‑1 (MQ=255)
                                                      tACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGc                                                                                                                                          >  1:792229/1‑97 (MQ=255)
                                                      tACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGc                                                                                                                                          >  1:792818/1‑97 (MQ=255)
                                                      tACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCgg                                                                                     >  2:302641/1‑150 (MQ=255)
                                                      tACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCgg                                                                                     >  2:35850/1‑150 (MQ=255)
                                                           ggCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCa                                                                                <  1:117375/150‑1 (MQ=255)
                                                              tGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGATTAACCCGCCAGTAACAACGGGCATTCCAGCATCAACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGTCTTCAGcc                                                                             >  2:83220/1‑150 (MQ=255)
                                                                     cagGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCa                                                                                >  1:849800/1‑140 (MQ=255)
                                                                     cagGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCa                                                                                <  2:849800/140‑1 (MQ=255)
                                                                     cagGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCa                                                                      <  1:871211/150‑1 (MQ=255)
                                                                          tGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCa                                                                 >  2:908422/1‑150 (MQ=255)
                                                                                 tACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAg                                                          <  2:476775/150‑1 (MQ=255)
                                                                                   ccGCTTCCGCTGCTTCACCCGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGCCGAGCCGTGCTTTGGGCGACTGATGGGAAAACAGCCCCAGGGCGCCCGGCTTCGGCCGTTCGCAGACCAGCGCAAGcc                                                        >  1:1026421/1‑149 (MQ=255)
                                                                                    cGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTc                                                       >  1:390741/1‑150 (MQ=255)
                                                                                       ttCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCtgtg                                                                                                                        >  1:508991/1‑82 (MQ=255)
                                                                                       ttCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCtgtg                                                                                                                        <  2:508991/82‑1 (MQ=255)
                                                                                           gctgctTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACt                                                                                                                 >  2:115279/1‑85 (MQ=255)
                                                                                           gctgctTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACt                                                                                                                 <  1:115279/85‑1 (MQ=255)
                                                                                            ctgctTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTggg                                                                                                                      <  1:646451/79‑1 (MQ=255)
                                                                                            ctgctTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTggg                                                                                                                      >  2:646451/1‑79 (MQ=255)
                                                                                                    ccagccagTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAg                                                                                                                                             <  2:751399/48‑1 (MQ=255)
                                                                                                    ccagccagTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAg                                                                                                                                             >  1:751399/1‑48 (MQ=255)
                                                                                                     aagccagTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGt                                      >  2:975350/2‑150 (MQ=255)
                                                                                                       gccagTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGTACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGGTGCCCGTTc                                    <  2:309900/150‑1 (MQ=255)
                                                                                                       gccagTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTc                                    <  1:302641/150‑1 (MQ=255)
                                                                                                          agTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGcc                                 <  2:390741/150‑1 (MQ=255)
                                                                                                                aacGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGCTGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTc                                    <  1:648887/141‑1 (MQ=255)
                                                                                                                aacGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAgg                                                                                             >  2:660506/1‑84 (MQ=255)
                                                                                                                aacGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAgg                                                                                             <  1:660506/84‑1 (MQ=255)
                                                                                                                aacGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTc                                    >  2:648887/1‑141 (MQ=255)
                                                                                                                  cGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGAt                         >  2:602347/1‑150 (MQ=255)
                                                                                                                   gggCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATc                        >  1:1263334/1‑150 (MQ=255)
                                                                                                                    ggaaTTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCa                       >  1:1426119/4‑150 (MQ=255)
                                                                                                                     gCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCAc                      >  1:1238047/1‑150 (MQ=255)
                                                                                                                          ccAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCgtg                 <  2:371765/150‑1 (MQ=255)
                                                                                                                                  gACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAg                                                                                              >  1:1364957/1‑65 (MQ=255)
                                                                                                                                  gACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAg                                                                                              <  2:1364957/65‑1 (MQ=255)
                                                                                                                                  gACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTc                                                                                 <  1:497213/78‑1 (MQ=255)
                                                                                                                                  gACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTc                                                                                 >  2:497213/1‑78 (MQ=255)
                                                                                                                                     gCCAGGCACTTCAAGCGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGc                                                                                    >  2:133813/1‑72 (MQ=38)
                                                                                                                                     gCCAGGCACTTCAAGCGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGc                                                                                    <  1:133813/72‑1 (MQ=37)
                                                                                                                                        aGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTTCGTCGCTGccc   <  2:1238047/150‑1 (MQ=255)
                                                                                                                                         ggCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTc                                    >  2:1176513/1‑116 (MQ=255)
                                                                                                                                         ggCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTc                                    <  1:1176513/116‑1 (MQ=255)
                                                                                                                                         ggCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTTCGTCGCTGCCCa  >  2:839992/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                              
CAGAAACTGCCGCAAGCGCTGGCGCATCACACTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTTCGTCGCTGCCCA  >  NC_002947/5253791‑5254077

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