Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,351,640 C→G 6.3% V33V (GTG→GTC PP_2066 ← GntR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,351,6400CG 6.3% 303.8 / 9.7 95V33V (GTG→GTCPP_2066GntR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (41/48);  new base G (3/3);  total (44/51)
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.50e-01

CGGGCTGCTGGCCTCCACCCGGGTGCCCACGCGCTTGCGCCGCGACACCAGCCCTTGAGTGAGCAACTGGTCGATAGCGGCACGCACCGTGTGCCGGCTCACATTGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCTTTGCCACGGTGGAGTAGCGCGTTGGGTTCATACCTTTGCTCGGGCTCGCCAATTTCATCGGTTGACAGTGTAACAGCGGCCGAATGCCGAACGCTCCGGCCGCCCACAG  >  NC_002947/2351493‑2351789
                                                                                                                                                   |                                                                                                                                                     
cgGGCTGCTGGCCTCCACCCGGGTGCCCACGCGCTTGCGCCGCGACACCAGCCCTTGAGTGAGCAACTGGTCGATAGCGGCACGCACCGTGTGCCGGCTCACATTGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCAc                                                                                                                                                     <  2:780423/150‑1 (MQ=255)
  ggatgctgGCCTCCACCCGGGTGCCCACGCGCTTGCGCCGCGACACCAGCCCTTGAGTGAGCAACTGGTCGATAGCGGCACGCACCGTGTGCCGGCTCACATTGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCg                                                                                                                                                   >  1:2336269/4‑150 (MQ=255)
  ggCTGCTGGCCTCCACCCGGGTGCCCACGCGCTTGCGCCGCGACACCAGCCCTTGAGTGAGCAACTGGTCGATAGCGGCACGCACCGTGTGCCGGCTCACATTGTAAAGCTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCg                                                                                                                                                   <  1:2225416/150‑1 (MQ=255)
   gctgctGGCCTCCACCCGGGTGCCCACGCGCTTGCGCCGCGACACCAGCCCTTGAGTGAGCAACTGGTCGATAGCGGCACGCACCGTGTGCCGGCTCACATTGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCgg                                                                                                                                                  <  2:2173254/150‑1 (MQ=255)
   gctgctGGCCTCCACCCGGGTGCCCACGCGCTTGCGCCGCGACACCAGCCCTTGAGTGAGCAACTGGTCGATAGCGGCACGCACCGTGTGCCGGCTCACATTGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCgg                                                                                                                                                  >  1:849247/1‑150 (MQ=255)
           tttCCAGCCGGGTGCCCACGCGCTTGCGCCGCAACACCAGCCCATAATTGAGCAACTGGTTGATAGCGGCAAGGCCGGTGTGCCGGCTCACATTGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGc                                                                                                                                          <  1:1685398/148‑1 (MQ=255)
           ccTCCACCCGGGTGCCCACGCGCTTGCGCCGCGACACCAGCCCTTGAGTGAGCAACTGGTCGATAGCGGCACGCACCGTGTGCCGGCTCACATTGTAAAGGTCGCACAGCTCGAACTCGGTGGGCAGCCTGCTGCCGACCGAGt                                                                                                                                                <  2:2124101/144‑1 (MQ=255)
           ccTCCACCCGGGTGCCCACGCGCTTGCGCCGCGACACCAGCCCTTGAGTGAGCAACTGGTCGATAGCGGCACGCACCGTGTGCCGGCTCACATTGTAAAGGTCGCACAGCTCGAACTCGGTGGGCAGCCTGCTGCCGACCGAGt                                                                                                                                                >  1:2124101/1‑144 (MQ=255)
            cTCCACCCGGGTGCCCACGCGCTTGCGCCGCGACACCAGCCCTTGAGTGAGCAACTGGTCGATAGCGGCACGCACCGTGTGCCGGCTCACATTGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGcc                                                                                                                                         >  2:1587644/1‑150 (MQ=255)
            cTCCACCCGGGTGCCCACGCGCTTGCGCCGCGACACCAGCCCTTGAGTGAGCAACTGGTCGATAGCGGCACGCACCGTGTGCCGGCTCACATTGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGcc                                                                                                                                         <  2:77148/150‑1 (MQ=255)
             tCCACCCGGGTGCCCACGCGCTTGCGCCGCGACACCAGCCCTTGAGTGAGCAACTGGTCGATAGCGGCACGCACAGTGTGCCGGCTCACATTGTAAAGGTCGCACAGCTCGAACTCGGTGGGCAGCCTGCTGCCGACCga                                                                                                                                                  <  2:469818/140‑2 (MQ=255)
             tCCACCCGGGTGCCCACGCGCTTGCGCCGCGACACCAGCCCTTGAGTGAGCAACTGGTCGATAGCGGCACGCACAGTGTGCCGGCTCACATTGTAAAGGTCGCACAGCTCGAACTCGGTGGGCAGCCTGCTGCCGACCga                                                                                                                                                  >  1:469818/1‑139 (MQ=255)
                aCCCGGGTGCCCACGCGCTTGCGCCGCGACACCAGCCCTTGAGTGAGCAACTGGTCGATAGCGGCACGCACCGTGTGCCGGCTCACATTGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTg                                                                                                                                     >  1:1253218/1‑150 (MQ=255)
                aCCCGGGTGCCCACGCGCTTGCGCCGCGACACCAGCCCTTGAGTGAGCAACTGGTCGATAGCGGCACGCACCGTGTGCCGGCTCACATTGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTg                                                                                                                                     <  2:1253218/150‑1 (MQ=255)
                aCCCGGGTGCCCACGCGCTTGCGCCGCGACACCAGCCCTTGAGTGAGCAACTGGTCGATAGCGGCACGCACCGTGTGCCGGCTCACATTGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTg                                                                                                                                     >  2:2120465/1‑150 (MQ=255)
                  ccGGGTGCCCACGCGCTTGCGCCGCGACACCAGCCCTTGAGTGAGCAACTGGTCGATAGCGGCACGCACCGTGTGCCGGCTCACATTGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCgg                                                                                                                                                  >  1:2126410/1‑135 (MQ=255)
                  ccGGGTGCCCACGCGCTTGCGCCGCGACACCAGCCCTTGAGTGAGCAACTGGTCGATAGCGGCACGCACCGTGTGCCGGCTCACATTGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCgg                                                                                                                                                  <  2:2126410/135‑1 (MQ=255)
                  ccGGGTGCCCACGCGCTTGCGCCGCGACACCAGCCCTTGAGTGAGCAACTGGTCGATAGCGGCACGCACCGTGTGCCGGCTCACATTGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTgg                                                                                                                                    <  2:1937733/149‑1 (MQ=255)
                  ccGGGTGCCCACGCGCTTGCGCCGCGACACCAGCACTTGAGTGAGCAACTGGTCGATAGCGGCACGCACCGTGTGCCGGCTCACATTGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTgg                                                                                                                                    >  1:1937733/1‑149 (MQ=255)
                      gTGCCCACGCGCTTGCGCCGCGACACCAGCCCTTGAGTGAGCAACTGGTCGATAGCGGCACGCACCGTGTGCCGGCTCACATTGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATc                                                                                                                               <  2:2052066/150‑1 (MQ=255)
                          ccACGCGCTTGCGCCGCGACACCAGCCCTTGAGTGAGCAACTGGTCGATAGCGGCACGCACCGTGTGCCGGCTCACATTGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCtt                                                                                                                           >  2:1954448/1‑150 (MQ=255)
                          ccACGCGCTTGCGCCGCGACACCAGCCCTTGAGTGAGCAACTGGTCGATAGCGGCACGCACCGTGTGCCGGCTCACATTGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCtt                                                                                                                           >  2:114181/1‑150 (MQ=255)
                          ccACGCGCTTGCGCCGCGACACCAGCCCTTGAGTGAGCAACTGGTCGATAGCGGCACGCACCGTGTGCCGGCTCACATTGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCtt                                                                                                                           >  1:805017/1‑150 (MQ=255)
                           cgcgcgCTTGCGCCGCGACACCAGCCCTTGAGTGAGCAACAGGTCAATACTGGCACTGATCGTGTGCCGGCTAACATTGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTg                                                                                                                          <  1:523827/148‑1 (MQ=255)
                           cACGCGCTTGCGCCGCGACACCAGCCCTTGAGTGAGCAACTGGTCGATAGCGGCACGCACCGTGTGCCGGCTCACATTGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTg                                                                                                                          <  2:517553/150‑1 (MQ=255)
                           cACGCGCTTGCGCCGCGACACCAGCCCTTGAGTGAGCAACTGGTCGATAGCGGCACGCACCGTGTGCCGGCTCACATTGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTg                                                                                                                          <  2:789090/150‑1 (MQ=255)
                            aCGCGCTTGCGCCGCGACACCAGCCCTTGAGTGAGCAACTGGTCGATAGCGGCACGCACCGTGTGCCGGCTCACATTGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGc                                                                                                                         >  1:296451/1‑150 (MQ=255)
                            aCGCGCTTGCGCCGCGACACCAGCCCTTGAGTGAGCAACTGGTCGATAGCGGCACGCACCGTGTGCCGGCTCACATTGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGc                                                                                                                         <  2:2336269/150‑1 (MQ=255)
                               cgcTTGCGCCGCGACACCAGCCCTTGAGTGAGCAACTGGTCGATAGCGGCACGCACCGTGTGCCGGCTCACATTGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATc                                                                                                                      <  1:1931690/150‑1 (MQ=255)
                               cgcTTGCGCCGCGACACCAGCCCTTGAGTGAGCAACTGGTCGATAGCGGCACGCACCGTGTGCCGGCTCACATTGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATc                                                                                                                      <  1:991281/150‑1 (MQ=255)
                               cgcTTGCGCCGCGACACCAGCCCTTGAGTGAGCAACTGGTCGATAGCGGCACGCACCGTGTGCCGGCTCACATTGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATc                                                                                                                      <  1:593329/150‑1 (MQ=255)
                                gcTTGCGCCGCGACACCAGCCCTTGAGTGAGCAACTGGTCGATAGCGGCACGCACCGTGTGCCGGCTCACATTGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCa                                                                                                                     >  1:1161285/1‑150 (MQ=255)
                                gcTTGCGCCGCGACACCAGCCCTTGAGTGAGCAACTGGTCGATAGCGGCACGCACCGTGTGCCGGCTCACATTGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCa                                                                                                                     >  1:1164001/1‑150 (MQ=255)
                                      gccgcGACACCAGCCCTTGAGTGAGCAACTGGTCGATAGCGGCACGCACCGTGTGCCGGCTCACATTGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCACTCGCTTGCATCAAGTCCt                                                                                                               >  1:1981826/1‑150 (MQ=255)
                                      gccgcGACACCAGCCCTTGAGTGAGCAACTGGTCGATAGCGGCACGCACCGTGTGCCGGCTCACATTGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCaa                                                                                                                    <  1:68088/145‑1 (MQ=255)
                                      gccgcGACACCAGCCCTTGAGTGAGCAACTGGTCGATAGCGGCACGCACCGTGTGCCGGCTCACATTGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCaa                                                                                                                    >  2:68088/1‑145 (MQ=255)
                                        cgcgACACCAGCCCTTGAGTGCGCAACTGGTCGATAGCGGCACGCACCGTGTGCCGGCTCACATTGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCtgt                                                                                                             <  2:2057118/150‑3 (MQ=255)
                                        cgcgACACCAGCCCTTGAGTGAGCAACTGGTCGATAGCGGCACGCACCGTGTGCCGGCTCACATTGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCttt                                                                                                             <  2:403012/150‑1 (MQ=255)
                                                cAGCCCTTGAGTGAGCAACTGGTCGATAGCGGCACGCACCGTGTGCCGGCTCACATTGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCTTTGCCACGGt                                                                                                     >  2:2202816/1‑150 (MQ=255)
                                                cAGCCCTTGAGTGAGCAACTGGTCGATAGCGGCACGCACCGTGTGCCGGCTCACATTGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCTTTGCCACGGt                                                                                                     >  2:1017018/1‑150 (MQ=255)
                                                 aGCCCTTGAGTGAGCAACTGGTCGATAGCGGCACGCACCGTGTGCCGGCTCACATTGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCTTTGCCACGGTg                                                                                                    <  1:2120465/150‑1 (MQ=255)
                                                 aGCCCTTGAGTGAGCAACTGGTCGATAGCGGCACGCACCGTGTGCCGGCTCACATTGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCTTTGCCACGGTg                                                                                                    <  1:2202816/150‑1 (MQ=255)
                                                    ccTTGAGTGAGCAACTGGTCGATAGCGGCACGCACCGTGTGCCGGCTCACATTGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCTTTGCCACGGTGGAg                                                                                                 <  1:602412/150‑1 (MQ=255)
                                                      ttGAGTGAGCAACTGGTCGATAGCGGCACGCACCGTGTGCCGGCTCACATTGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCTTTGCCACGGTGGAGTa                                                                                               <  1:2101997/150‑1 (MQ=255)
                                                        gagtgagCAACTGGTCGCTAGCGGCTCGCACCGTGTGCCGGCTCACATTGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCTTTGCCACGGTGGAGTAgc                                                                                             <  2:1705286/150‑1 (MQ=255)
                                                        gagtgagCAACTGGTCGATAGCGGCACGCACCGTGTGCCGGCTCACATTGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCTTTGCCACGGTGGAGTAgc                                                                                             <  1:2020320/150‑1 (MQ=255)
                                                           tgagCAACTGGTCGATAGCGGCACGCACCGTGTGCCGGCTCACATTGTAAAGGTCGCACACCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCTTTGCCACGGTGGAGTAGcgcg                                                                                          >  2:680283/1‑150 (MQ=255)
                                                             agCAACTGGTCGATAGCGGCACGCACCGTGTGCCGGCTCACATTGTAAAGGTCGCACCACTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCTTTGCCACGGTGGAGTAGCGCGtt                                                                                        <  1:384338/150‑1 (MQ=255)
                                                             agCAACTGGTCGATAGCGGCACGCACCGTGTGCCGGCTCACATTGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCTTTGCCACGGTGGAGTAGCGCGtt                                                                                        <  1:1921998/150‑1 (MQ=255)
                                                              gCAACTGGTCGATAGCGGCACGCACCGTGTGCCGGCTCACATTGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCTTTGCCACGGTGGAGTAGCGCGTTg                                                                                       <  2:1161285/150‑1 (MQ=255)
                                                              gCAACTGGTCGATAGCGGCACGCACCGTGTGCCGGCTCACATTGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCTTTGCCACGGTGGAGTAGCGCGTTg                                                                                       <  2:1164001/150‑1 (MQ=255)
                                                              gCAACTGGTCGATAGCGGCACGCACCGTGTGCCGGCTCACATTGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCTTTGCCACGGTGGAGTAGCGCGTTg                                                                                       <  2:594447/150‑1 (MQ=255)
                                                                        gATAGCGGCACGCACCGTGTGCCGGCTCACATTGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCTTTGCCACGGTGGAGTAGCGCGTTGGGTTCATAcc                                                                             >  2:2091200/1‑150 (MQ=255)
                                                                               gcacgcacCGTGTGCCGGCTCACATTGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCTTTGCCACGGTGGAGTAGCGCGTTGGGTTCATACCTTTGCTc                                                                      <  1:2091200/150‑1 (MQ=255)
                                                                               gcacgcacCGTGTGCCGGCTCACATTGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCTTTGCCACGGTGGAGTAGCGCGTTGGGTTCATACCTTTGCTc                                                                      <  2:296451/150‑1 (MQ=255)
                                                                                   gcacCGTGTGCCGGCTCACATTGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCTTTGCCACGGTGGAGTAGCGCGTTGGGTTCATACCTTTGCTCGGGc                                                                  <  2:805017/150‑1 (MQ=255)
                                                                                     acCGTGTGCCGGCTCACATTGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCTTTGCCACGGTGGAGTAGCGCGTTgg                                                                                      <  2:2211898/128‑1 (MQ=255)
                                                                                     acCGTGTGCCGGCTCACATTGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCTTTGCCACGGTGGAGTAGCGCGTTgg                                                                                      >  1:2211898/1‑128 (MQ=255)
                                                                                       cGTGTGCCGGCTCACATTGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCaa                                                                                                                                 <  1:2158587/83‑1 (MQ=255)
                                                                                       cGTGTGCCGGCTCACATTGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCaa                                                                                                                                 >  2:2158587/1‑83 (MQ=255)
                                                                                         tgtgGCGGCGAACATGGTAAAGGTCGCACAGCTCGAACTCGGTGGGCTGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCTTTGCCACGGTGGAGTAGCGCGTTGGGTTCATACCTTTGCTCGGGCTCGCCa                                                            >  2:1256298/1‑150 (MQ=255)
                                                                                            gCCGGCTCACATTGTAAAGGTCGCACACCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATTAAGTCCTTTGCCACGGTGGAGTAGCGCGTTGGGTTCATACCTTTGCTCGGGCTCGCCAAtt                                                         <  1:680283/150‑1 (MQ=255)
                                                                                             ccGGCTCACATTGTAAAGGTCGCACAGCTCGGACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCTTTGCCACGGTGGAGTAGCGCGTTGGGTTCATACCTTTGCTCGGGCTCGCCAAttt                                                        <  2:772603/150‑1 (MQ=255)
                                                                                                    acaTTGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCTTTGCCACGGTGGAGTAGCGCGTTGGGTTCATACCTTTGCTCGGGCTCGCCAATTTCATCGGt                                                 >  1:677987/1‑150 (MQ=255)
                                                                                                     caTTGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCCCCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCTTTGCCACGGTGGAGTAGCGCGTTGGGTTCATACCTTTGCTCGGGCTCGCCAATTTCATCGGtt                                                <  2:849247/150‑1 (MQ=255)
                                                                                                      aTTGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGATGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCTTTGCCACGGTGGAGTAGCGCGTTGGGTTCATACCTTTGCTCGGGCTCGCCAATTTCATCGGTTg                                               >  2:894407/1‑150 (MQ=255)
                                                                                                       ttGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGACGCTGGCAATCGCTTGCATCAAGTCCTTTGCCACGGTGGAGTAGCGCGTTGGGTTCATACCTTTGCTCGGGCTCGCCAATTTCATCGGTTGa                                              >  2:904507/1‑150 (MQ=255)
                                                                                                         gTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCTTTGCCACGGTGGAGTAGCGCGTTGGGTTCATACCTTTGCTCGGGCTCGCCAATTTCATCGGTTGACa                                            <  1:1017018/150‑1 (MQ=255)
                                                                                                         gTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGATGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCTTTGCCACGGTGGAGTAGCGCGTTGGGTTCATACCTTTGCTCGGGCTCGCCAATTTCATCGGTTGACa                                            <  1:894407/150‑1 (MQ=255)
                                                                                                              ggTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCTTTGCCACGGTGGAGTAGCGCGTTGGGTTCATACCTTTGCTCGGGCTCGCCAATTTCATCGGTTGACAGTGTa                                       >  2:2110216/1‑150 (MQ=255)
                                                                                                                  gCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCTTTGCCACGGTGGAg                                                                                                 <  1:1077439/88‑1 (MQ=255)
                                                                                                                  gCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCTTTGCCACGGTGGAg                                                                                                 >  2:1077439/1‑88 (MQ=255)
                                                                                                                  gCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCTTTGCCACGGTGGAGTAGCGCGTTGGGTTCATACCTTTGCTCGGGCTCGCCAATTTCATCGGTTGACAGTGTAACAg                                   >  2:835482/1‑150 (MQ=255)
                                                                                                                    acaGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCTTTGCCACGGTGGAGTAGCGCGTTGGGTTCATACCt                                                                            >  2:370752/1‑107 (MQ=255)
                                                                                                                    acaGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCTTTGCCACGGTGGAGTAGCGCGTTGGGTTCATACCt                                                                            <  1:370752/107‑1 (MQ=255)
                                                                                                                    acaGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCTTTGCCACGGTGGAGTAGCGCGTTGGGTTCATACCTTTGCTCGGGCTCGCCAATTTCATCGGTTGACAGTGTAACAGCg                                 >  2:745411/1‑150 (MQ=255)
                                                                                                                    acaGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCTTTGCCACGGTGGAGTAGCGCGTTGGGTTCATACCTTTGCTCGGGCTCGCCAATTTCATCGGTTGACAGTGTAACAGCg                                 >  2:1025997/1‑150 (MQ=255)
                                                                                                                       gCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCTTTGCCACGGTGGAGTAGCGCGTTGGGTTCATACCTTTGCTCGGGCTCGCCAATTTCATCGGTTGACAGTGTAACAGCGGcc                              <  2:216372/150‑1 (MQ=255)
                                                                                                                           gAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCTTTGCCACGGTGGAGTAGCGCGTTGGGTTCATACCTTTGCTCGGGCTCGCCAATTTCATCGGtt                                                <  2:1754127/128‑1 (MQ=255)
                                                                                                                           gAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCTTTGCCACGGTGGAGTAGCGCGTTGGGTTCATACCTTTGCTCGGGCTCGCCAATTTCATCGGtt                                                >  1:1754127/1‑128 (MQ=255)
                                                                                                                           gAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCTTTGCCACGGTGGAGTAGCGCGTTGGGTTCATACCTTTGCTCGGGCTCGCCAATTTCATCGGTTGACAGTGTAACAGCGGCCGAAt                          <  2:1070828/150‑1 (MQ=255)
                                                                                                                           gAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCTTTGCCACGGTGGAGTAGCGCGTTGGGTTCATACCTTTGCTCGGGCTCGCCAATTTCATCGGTTGACAGTGTAACAGCGGCCGAAt                          >  2:1020178/1‑150 (MQ=255)
                                                                                                                             aCTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCTTTGCCACGGTGGAGTAGCGCGTTGGGTTCATACCTTTGCTc                                                                      <  2:731850/104‑1 (MQ=255)
                                                                                                                             aCTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCTTTGCCACGGTGGAGTAGCGCGTTGGGTTCATACCTTTGCTc                                                                      >  1:731850/1‑104 (MQ=255)
                                                                                                                             aCTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCTTTGCCACGGTGGAGTAGCGCGTTGGGTTCATACCTTTGCTCGGGCTCGCCAATTTCATCGGTTGACAGTGTAACAGCGGCCGAATGc                        >  1:1448496/1‑150 (MQ=255)
                                                                                                                               tcggtcggCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCTTTGCCACGGTGGAGTAGCGCGTTGGGTTCATACCTTTGCTCGGGCTCGCCAATTTCATCgg                                                  <  2:1046627/122‑1 (MQ=255)
                                                                                                                               tcggtcggCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCTTTGCCACGGTGGAGTAGCGCGTTGGGTTCATACCTTTGCTCGGGCTCGCCAATTTCATCgg                                                  >  1:1046627/1‑122 (MQ=255)
                                                                                                                                cggtgggCAGCCTGCTGCCGACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCTTTGCCACGGTGGAGTAgc                                                                                             >  2:1888812/1‑78 (MQ=37)
                                                                                                                                cggtgggCAGCCTGCTGCCGACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCTTTGCCACGGTGGAGTAgc                                                                                             <  1:1888812/78‑1 (MQ=37)
                                                                                                                                      gcagcaGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCTTTGCCACGGTGGAGTAGCGCGTTgg                                                                                      >  2:1173222/1‑79 (MQ=255)
                                                                                                                                      gcagcaGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCTTTGCCACGGTGGAGTAGCGCGTTgg                                                                                      <  1:1173222/79‑1 (MQ=255)
                                                                                                                                        agcagGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCTTTGCCACGGTGGAGTAGCGCGTTGGGTTCATACCTTTGCTCGGGCTCGCCAATTTCATCGGTTGACAGTGTAACAGCGGCCGAATGCCGAACGCTCCg             >  1:864243/1‑150 (MQ=255)
                                                                                                                                            ggCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCTTTGCCACGGTGGAGTAGCGCGTTGGGTTCATACCTTTGCTCGGGCTCGCCAATTTCATCGGTTGACAGTGTAACAGCGGCCGAATGCCGAACGCTCCGgccg         >  2:1216251/1‑150 (MQ=255)
                                                                                                                                             gCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCTTTGCCACGGTGGAGTAGCGCGTTGGGTTCATACCTTTGCTCGGGCTCGCCAATTTCATCGGTTGACAGTGTAACAGCGGCCGAATGCCGAACGCTCCGgccgc        >  2:750645/1‑150 (MQ=255)
                                                                                                                                             gCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCTTTGCCACGGTGGAGTAGCGCGTTGGGTTCATACCTTTGCTCGGGCTCGCCAATTTCATCGGTTGACAGTGTAACAGCGGCCGAATGCCGAACGCTCCGgccgc        >  2:750604/1‑150 (MQ=255)
                                                                                                                                                  ccACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCTTTGCCACGGTGGAGTAGCGCGTTGGGTTCATACCTTTGCTCGGGCTCGCCAATTTCATCGGTTGACAGTGTAACAGCGGCCGAATGCCGAACGCTCCGGCCGCCcaca   >  1:1241371/1‑150 (MQ=255)
                                                                                                                                                   cACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTc                                                                                                                 >  1:72180/1‑39 (MQ=255)
                                                                                                                                                   cACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTc                                                                                                                 <  2:72180/39‑1 (MQ=255)
                                                                                                                                                   cACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCTTTGCCACGGTGGAGAAGCCAGTTGGGTTCATACCGTTACTCGGGCTCGCCAATTTCATCGGTTGACAGTGTAACAGCGGCCGAATGCCGAACGCTCCGGCCGCCCACAg  >  1:870894/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CGGGCTGCTGGCCTCCACCCGGGTGCCCACGCGCTTGCGCCGCGACACCAGCCCTTGAGTGAGCAACTGGTCGATAGCGGCACGCACCGTGTGCCGGCTCACATTGTAAAGGTCGCACAGCTCGAACTCGGTCGGCAGCAGGCTGCCCACCGGGTAGCGGCCGCTGGCAATCGCTTGCATCAAGTCCTTTGCCACGGTGGAGTAGCGCGTTGGGTTCATACCTTTGCTCGGGCTCGCCAATTTCATCGGTTGACAGTGTAACAGCGGCCGAATGCCGAACGCTCCGGCCGCCCACAG  >  NC_002947/2351493‑2351789

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