Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,748,399 T→C 5.8% intergenic (‑227/+20) fbp ← / ← glgP class 1 fructose‑1,6‑bisphosphatase/glycogen phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,748,3990TC 5.8% 804.6 / 24.3 259intergenic (‑227/+20)fbp/glgPclass 1 fructose‑1,6‑bisphosphatase/glycogen phosphorylase
Reads supporting (aligned to +/- strand):  ref base T (139/101);  new base C (12/3);  total (151/108)
Fisher's exact test for biased strand distribution p-value = 1.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.86e-01

Only 100 of 267 total aligned reads displayed.

ACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTGCGCCACCACTCCTGCGGTGTGTGCCACAGATCCTCG  >  NC_002947/5748250‑5748548
                                                                                                                                                     |                                                                                                                                                     
aCTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGt                                                                                                                                                       >  2:3398806/1‑150 (MQ=255)
 ggctggttctcgtcggaccgcaccattcgcgaatatgccaccgagatctggaaggcattggattgaggtgaggccattcgcgggcgcgctCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTg                                                                                                                                                      >  1:1558451‑M1/91‑150 (MQ=255)
  tCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGc                                                                                                                                                     >  2:307589/1‑150 (MQ=255)
  tCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGAGGGCGTGc                                                                                                                                                     >  1:6180971/1‑150 (MQ=255)
  gctggttctcgtcggaccgcaccattcgcgaatatgccaccgagatctggaaggcattggattgaggtgaggccgttcgcgggcacgcccgcgcccccacgtgtctcgcaagcgtcaaggcttgcgctgttccggtgggggcgGgcgcgc                                                                                                                                                     >  1:1822673‑M2/144‑150 (MQ=255)
    gCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGccc                                                                                                                                                   <  2:5417681/150‑1 (MQ=255)
       aCGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCcgcg                                                                                                                                                >  1:520761/1‑150 (MQ=255)
              tcAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGcc                                                                                                                                         >  2:4231012/1‑150 (MQ=255)
                 gACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGcctca                                                                                                                                      >  1:2694888/1‑150 (MQ=255)
                 gACGGCACAATAGGTCGGTAGTTCAGCCTGCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGcctca                                                                                                                                      <  2:4707757/150‑1 (MQ=255)
                    ggCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGcctcacct                                                                                                                                   >  1:4864278/1‑150 (MQ=255)
                       acaATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCaa                                                                                                                                >  2:4377560/1‑150 (MQ=255)
                        caATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAAt                                                                                                                               >  2:576604/1‑150 (MQ=255)
                        caATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAAt                                                                                                                               >  2:1157270/1‑150 (MQ=255)
                        caATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAAt                                                                                                                               <  2:520761/150‑1 (MQ=255)
                                  gTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGACCGGCCTCACCTCAATCCAATGCCtt                                                                                                                     >  2:2837966/1‑150 (MQ=255)
                                       caccgagatctggaaggcattggattgaggtgtgtccattcgcgggcgcgctCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGa                                                                                                                <  2:5086508‑M1/98‑1 (MQ=255)
                                        accgagatctggaaggcattggattgaggtgaggccattcgcgggcgcgctCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGt                                                                                                                                                       <  1:746199‑M1/59‑1 (MQ=255)
                                           ccTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGAtctc                                                                                                            >  1:6703018/1‑150 (MQ=255)
                                           ccTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGAtctc                                                                                                            >  1:5152860/1‑150 (MQ=255)
                                            ctccGCCAGGCCGGTGGTCTGGTTGTGTGACCTTTCGCGGGCGCGCGCGCCCCCCCCGGTACCCCGAAAGCCTTGACGTTTGCGCTACACCTGTGGGAGCGGGCGTCCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCg                                                                                                           <  1:1942105/147‑1 (MQ=255)
                                            cTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCCCAGGTACAGCGCAAGCCTTGACGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCg                                                                                                           <  2:4075103/150‑1 (MQ=255)
                                            cTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCg                                                                                                           <  2:2373327/150‑1 (MQ=255)
                                              tCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCGCGCTCCCACTGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGTGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGt                                                                                                         <  1:1067020/150‑1 (MQ=255)
                                               cGCCCGGCGGGGGGTGGGGTTGTGGGCCCATCCGGGGGCGCGCCCGCCCCCACCGGCACCGCCCAAGCCTTGATGGTGGCGAGACCCGTGTGGGCGGGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTg                                                                                                        <  1:4231012/150‑1 (MQ=255)
                                                gCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTgg                                                                                                       >  2:5185867/1‑150 (MQ=255)
                                                 ccAGGCCGGTGGTCTGGTTGTGTTTCCTTTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCACGCCTTGAAGCTTGCGAGACACCTGTGGTAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGc                                                                                                      <  1:3273985/150‑1 (MQ=255)
                                                           ggtCTGGTTGGGTGTCCATTGGCGGGCGCGCTCGCTCCCACAGGTACAGCGCACGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGa                                                                                            <  2:357236/150‑1 (MQ=255)
                                                             tCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATAGTCGCGAAt                                                                                          >  2:4361425/1‑150 (MQ=255)
                                                              gtGGTTGTGTGTCCTTTCGGGGGCGCGCTCGCTCCCACAGGTACCGCGAACGCATTCATGCTGGCGAGACACGTGTGGGAGGGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATg                                                                                         <  1:1378385/149‑1 (MQ=255)
                                                              cTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATg                                                                                         <  2:1308103/150‑1 (MQ=255)
                                                                 gTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTg                                                                                      >  1:1049107/1‑150 (MQ=255)
                                                                   tgtgtgTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCg                                                                                    <  1:2235103/150‑1 (MQ=255)
                                                                      gtgtCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTc                                                                                 >  1:3884987/1‑150 (MQ=255)
                                                                          ccATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCcgac                                                                             >  2:179699/1‑150 (MQ=255)
                                                                           cATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCcgacg                                                                            >  2:3143387/1‑150 (MQ=255)
                                                                                gcgGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAAc                                                                       >  1:4991728/1‑150 (MQ=255)
                                                                                  gggCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCa                                                                     <  1:6773866/150‑1 (MQ=255)
                                                                                  gggCGCGCGCGCTCCCACCGGTACAGCGCACGCCTGGACGCGTGCGCGCCCCGGGGGGGAGCGGGCGGGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCa                                                                     <  1:5339914/150‑1 (MQ=255)
                                                                                   ggCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGACCCAg                                                                    >  2:4941828/1‑150 (MQ=255)
                                                                                   ggCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAg                                                                    <  2:6851140/150‑1 (MQ=255)
                                                                                      acgcccgctcccacaggtgtctcgcaagcttcaaggcttgcgctgtacctgtgggagcgAGCGCGCCCGCGAATGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGt                                                                                                         >  2:746199‑M2/60‑110 (MQ=255)
                                                                                        gctcgctcCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCAt                                                               >  2:1533361/1‑150 (MQ=255)
                                                                                              tcccacaggtgtctcgcaagcttcaaggcttgcgctgtacctgtgggagcgGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAAt                                                                                          <  1:6368161‑M2/66‑1 (MQ=255)
                                                                                              TCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATccgtgc                                                         >  2:3998527‑M1/1‑150 (MQ=255)
                                                                                              TCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCACGAATGGTGCGGTCCGACGAGAACCAGCCCATccgtgc                                                         <  1:6022492‑M1/150‑1 (MQ=255)
                                                                                                ccacaggtgtctcgcaagcttcaaggcttgcgctgtacctgtgggagcgAGCGCGCCCGCGAATGGACACACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATccgtgccg                                                       >  1:253852‑M2/50‑150 (MQ=255)
                                                                                                CCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATccgtgccg                                                       <  2:4823102‑M1/150‑1 (MQ=255)
                                                                                                CCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATccgtgccg                                                       <  2:4534904‑M1/150‑1 (MQ=255)
                                                                                                CCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATccgtgccg                                                       <  2:6018992‑M1/150‑1 (MQ=255)
                                                                                                CCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATccgtgccg                                                       <  1:3933944‑M1/150‑1 (MQ=255)
                                                                                                 cacaggtgtctcgcaagcttcaaggcttgcgctgtacctgtgggagcgAGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCAgccgtgccgg                                                      >  2:2560193‑M2/49‑149 (MQ=255)
                                                                                                 CACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATccgtgccgt                                                      <  1:3550090‑M1/150‑1 (MQ=255)
                                                                                                 CACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATccgtgccgt                                                      <  1:1233580‑M1/150‑1 (MQ=255)
                                                                                                 CACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATccgtgccgt                                                      <  1:1717785‑M1/150‑1 (MQ=255)
                                                                                                  ACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATccgtgccgtg                                                     <  2:1650945‑M1/150‑1 (MQ=255)
                                                                                                  ACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATccgtgccgtg                                                     <  1:1171736‑M1/150‑1 (MQ=255)
                                                                                                  ACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATccgtgccgtg                                                     <  1:4590279‑M1/150‑1 (MQ=255)
                                                                                                    aggtgtctcgcaagcttcaaggcttgcgctgtacctgtgggagcgAGCGCGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGtt                                                   >  2:6600590‑M2/46‑150 (MQ=255)
                                                                                                    aggtgtctcgcaagcttcaaggcttgcgctgtacctgtgggagcgAGCGCGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGACCGACGAGAACCAGCCCATCCGTGCCGTGtt                                                   >  2:1915268‑M2/46‑150 (MQ=255)
                                                                                                    AGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGtt                                                   <  2:4575805‑M1/150‑1 (MQ=255)
                                                                                                      gtgtctcgcaagcttcaaggcttgcgctgtacctgtgggagcgAGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCa                                                 >  1:3687322‑M2/44‑150 (MQ=255)
                                                                                                      GTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCCTATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGGGTTCa                                                 <  2:5545190‑M1/150‑1 (MQ=255)
                                                                                                       TACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCaa                                                <  2:6518762‑M1/150‑1 (MQ=255)
                                                                                                       TACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCaa                                                <  2:6173566‑M1/150‑1 (MQ=255)
                                                                                                        aCAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAAt                                               >  1:6454133/1‑150 (MQ=255)
                                                                                                         tctcgcaagcttcaaggcttgcgcgacacctgtgggagcgGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATa                                              >  1:6653969‑M2/41‑150 (MQ=255)
                                                                                                           tcgcaagcttcaaggcttgcgctgcacctgtgggagcgGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATAcc                                            >  2:1095110‑M2/39‑150 (MQ=255)
                                                                                                           tcgcaagcttcaaggcttgcgagacacctgtgggagcgGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATAcc                                            >  1:2404958‑M2/39‑150 (MQ=255)
                                                                                                           gcgcAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCAGCTTCGCGACTGGTGCGGTCCGACGAGAACCAGCCCATCCGGGCCGGGTTCAATAcc                                            >  2:4665470/1‑150 (MQ=255)
                                                                                                            cgcAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCg                                           >  1:5489226/1‑150 (MQ=255)
                                                                                                              cAAGCCTTTACGCTTGCGAGACACGTGTGGGCGGGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGcc                                         <  1:3089475/150‑1 (MQ=255)
                                                                                                               aaGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCt                                        >  2:5232922/1‑150 (MQ=255)
                                                                                                               aaGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCt                                        >  1:5610663/1‑150 (MQ=255)
                                                                                                               aaGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCt                                        >  2:1598842/1‑150 (MQ=255)
                                                                                                                   cTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGaa                                                                        >  1:1641070/1‑114 (MQ=255)
                                                                                                                   cTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTgcg                                    >  1:2864644/1‑150 (MQ=255)
                                                                                                                    ttGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTgcgc                                   >  1:4099239/1‑150 (MQ=255)
                                                                                                                    ttGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTgcgc                                   >  2:675422/1‑150 (MQ=255)
                                                                                                                      gAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTGCGcca                                 <  1:4381613/150‑1 (MQ=255)
                                                                                                                       taGCTTGCGAGACTCCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTGCGccac                                <  2:3893960/149‑1 (MQ=255)
                                                                                                                       aaGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTGCGccac                                <  2:1851434/150‑1 (MQ=255)
                                                                                                                       aaGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTGCGccac                                <  1:4720487/150‑1 (MQ=255)
                                                                                                                            ggCGAGACACCTGTGGCCGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTGCGCCACCACTc                           <  1:3319631/149‑1 (MQ=255)
                                                                                                                               gagaCACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTGCGCCACCACTCCTg                        >  1:4166253/1‑150 (MQ=255)
                                                                                                                               gagaCACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTGCGCCACCACTCCTg                        >  1:4281295/1‑150 (MQ=255)
                                                                                                                                    acctgtgggagcgGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTGCGCCACCACTCCTGCGgtg                   >  2:4102665‑M2/14‑150 (MQ=255)
                                                                                                                                     cctgtgggagcgGGCGTGCCCGCGAACGGCCTCACCTCACTCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTGCGCCACCACTCCTGCGgtgt                  <  1:2252898‑M2/138‑1 (MQ=255)
                                                                                                                                     cctgtgggagcgGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTGCGCCACCACTCCTGCGgtgt                  >  2:2143739‑M2/13‑150 (MQ=255)
                                                                                                                                      ctgtgggagcgGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTGCGCCACCACTCCTGCGgtgtg                 >  2:984775‑M2/12‑150 (MQ=255)
                                                                                                                                      ctgtgggagcgGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCACCTTGCGCCACCACTCCTGCGgtgtg                 >  1:79236‑M2/12‑150 (MQ=255)
                                                                                                                                       tgtgggagcgGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTGCGCCACCACTCCTGCGgtgtgt                <  2:6703018‑M2/140‑1 (MQ=255)
                                                                                                                                         tgggagcgGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTGCGCCACCACTCCTGCGGTGTGTGc              <  1:6628277‑M2/142‑1 (MQ=255)
                                                                                                                                          gttcgcgGGCACGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTGCGCCAGCCCTCGTGCGGTGTGTgcg             >  2:6423211‑M2/8‑149 (MQ=255)
                                                                                                                                           ttcgcgGGCACGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATAATCGCGAATGGTGCGGTCCGACGAGACCCAGCCCATCCGTGCCGTGTTCAATACCGCCTGGCGCCACCACTCCTGCGGTGTGTGCca            >  2:2548113‑M2/7‑150 (MQ=255)
                                                                                                                                             agcgGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTGCGCCACCACTCCTGCGGTGTGTGCcaca          >  1:2203449‑M2/5‑150 (MQ=255)
                                                                                                                                                gGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTGCGCCACCACTCCTGCGGTGTGTGCCACAGAt       <  2:5941956‑M2/149‑1 (MQ=255)
                                                                                                                                                gGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTGCGCCACCACTCCTGCGGTGTGTGCCACAGAt       >  2:2208677‑M2/2‑150 (MQ=255)
                                                                                                                                                 ggCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTGCGCCACCACTCCTGCGGTGTGTGCCACAGATc      <  2:5610663/150‑1 (MQ=255)
                                                                                                                                                  gCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATCTTCGCGCATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTGCGCCACCACTCCTGCGGTGTGTGCCCCAGATcc     >  2:5791735/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTGCGCCACCACTCCTGCGGTGTGTGCCACAGATCCTCG  >  NC_002947/5748250‑5748548

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