Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,748,340 T→C 9.6% intergenic (‑168/+79) 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,3400TC 9.6% 729.9 / 43.9 251intergenic (‑168/+79)fbp/glgPclass 1 fructose‑1,6‑bisphosphatase/glycogen phosphorylase
Reads supporting (aligned to +/- strand):  ref base T (121/102);  new base C (15/9);  total (137/114)
Fisher's exact test for biased strand distribution p-value = 5.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 267 total aligned reads displayed.

TGGGGGGATGAAAAACCCACGCAGTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCC  >  NC_002947/5748191‑5748489
                                                                                                                                                     |                                                                                                                                                     
tGGGGGGATGAAAAACCCACGCAGTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGcgct                                                                                                                                                       <  1:6953623/150‑1 (MQ=255)
  gggggATGAACAACCCACGCAGTTTACAGGGAGAGTGACGGGCTAGCGAGTGTTGAGACTCGCCACGGCTCTCAGACGGCACCATAGGTCGGTGGTTCGGCCTTCGCCAGGCCGGGGGTGTGGTTGTGTGTCCAGTCGCGGGCGcgctcg                                                                                                                                                     >  1:5009963/1‑150 (MQ=255)
  gggggATGAAAAACCCACGCAGTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGGAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGcgctcg                                                                                                                                                     >  2:5401231/1‑150 (MQ=255)
   ggggATGAAAAACCCACGCAGTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGcgctcgc                                                                                                                                                    <  1:5108916/150‑1 (MQ=255)
    gggATGAAAAACCCACGCAGTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGcgctcgct                                                                                                                                                   <  2:599978/150‑1 (MQ=255)
        tGAAAAACCCACGCAGTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCca                                                                                                                                               <  1:2272556/150‑1 (MQ=255)
           aaaaCCCACGCAGTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAg                                                                                                                                            >  2:5693852/1‑150 (MQ=255)
            aaaCCCACGCAGTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAgg                                                                                                                                           >  2:4449697/1‑150 (MQ=255)
               cccACGCAGTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTAc                                                                                                                                        <  1:759954/150‑1 (MQ=255)
                ccACGCAGTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACa                                                                                                                                       <  2:3368387/150‑1 (MQ=255)
                 taCGCAGTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAg                                                                                                                                      <  1:4670683/149‑1 (MQ=255)
                 cACGCAGTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAg                                                                                                                                      >  2:2607184/1‑150 (MQ=255)
                      aGTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCaa                                                                                                                                 <  1:724349/150‑1 (MQ=255)
                        tttACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGc                                                                                                                               >  1:2914832/1‑150 (MQ=255)
                         ttACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGcc                                                                                                                              >  2:1025289/1‑150 (MQ=255)
                              gggAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGaa                                                                                                                         <  2:6902453/150‑1 (MQ=255)
                              gggAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGaa                                                                                                                         <  1:5795467/150‑1 (MQ=255)
                              gggAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGaa                                                                                                                         <  2:43456/150‑1 (MQ=255)
                               ggAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAg                                                                                                                        >  1:3145519/1‑150 (MQ=255)
                                  gagTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCATGAAGCtt                                                                                                                     <  1:534123/150‑1 (MQ=255)
                                   agTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTg                                                                                                                    <  1:1375910/150‑1 (MQ=255)
                                   agTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTg                                                                                                                    <  1:5969737/150‑1 (MQ=255)
                                   agTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTg                                                                                                                    <  1:1376183/150‑1 (MQ=255)
                                      ggtattgaacacggcacggatgggctggttctcgtcggaccgcaccattcgcgaatatgccaccgagatctggaaggcattggattgaggtgaggccgTTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCga                                                                                                                 >  2:1324040‑M1/99‑150 (MQ=255)
                                      gACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCga                                                                                                                 >  2:3424740/1‑150 (MQ=255)
                                      gACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCga                                                                                                                 <  1:6172963/150‑1 (MQ=255)
                                      gACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCga                                                                                                                 <  1:4936595/150‑1 (MQ=255)
                                       aCGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCTTCAAGGCTTGCgct                                                                                                                <  2:641161/150‑3 (MQ=255)
                                       aCGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCgag                                                                                                                >  1:4029104/1‑150 (MQ=255)
                                        cGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCgaga                                                                                                               <  1:2682496/150‑1 (MQ=255)
                                         gggCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGac                                                                                                              >  1:1368385/1‑150 (MQ=255)
                                           gCTAGCGCGAGAGGCGACCCGCCCCCCCCCGCAGACGTCACACTCGCACGGGAGTACAGCCTTCGCCAGGCCGGTGGTCGGCTTGGGCGCCCATTCGCGGGCGCGCGCGCTCCCACAGGTACAGCGCAAGCCATGAAGCTTGCGAGacac                                                                                                            <  2:5009963/150‑1 (MQ=255)
                                           gCTAGCGAGTGATGAGACTCGCCAGGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGacac                                                                                                            <  2:4513743/150‑1 (MQ=255)
                                           gCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCATCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTGTCTCGCAAGCTTCAAGGCTTGCgctgtac                                                                                                            <  2:4665486/150‑7 (MQ=25)
                                             tAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGGGAGAGACCt                                                                                                          >  1:3919733/1‑150 (MQ=255)
                                               gCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCtgt                                                                                                        <  1:3107187/150‑1 (MQ=255)
                                               gCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCtgt                                                                                                        <  2:7094683/150‑1 (MQ=255)
                                                cGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCtgtg                                                                                                       <  2:3982589/150‑1 (MQ=255)
                                                acggcacggatgggctggttctcgtcggaccgcaccattcgcgaatatgccaccgagatctggaaggcattggattgaggtgaggccgTTCGCGGGCACGCCCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCtgtg                                                                                                       >  1:6036635‑M1/89‑150 (MQ=255)
                                                 gAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTgg                                                                                                      <  1:2187252/150‑1 (MQ=255)
                                                 cggcacggatgggctggttctcgtcggaccgcaccattcgcgaatatgccaccgagatctggaaggcattggattgaggtgaggccgTTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTgg                                                                                                      >  2:5905828‑M1/88‑150 (MQ=255)
                                                  ggcacggatgggctggttctcgtcggaccgcaccattcgcgaatatgccaccgagatctggaaggcattggattgaggtgaggccgTTCGCGGGCACGCCCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCga                                                                                                                 <  2:4192965‑M1/52‑1 (MQ=255)
                                                  aGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTggg                                                                                                     <  1:363374/150‑1 (MQ=255)
                                                   gcacggatgggctggttctcgtcggaccgcaccattcgcgaatatgccaccgagatctggaaggcattggattgaggtgaggccgTTCGCGGGCACGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGa                                                                                                    >  1:1242115‑M1/86‑150 (MQ=255)
                                                   gcacggatgggctggttctcgtcggaccgcaccattcgcgaatatgccaccgagatctggaaggcattggattgaggtgaggccgTTCGCGGGCACGCCCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGa                                                                                                    >  2:3005535‑M1/86‑150 (MQ=255)
                                                   gTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTGTCTCGCAAGCTTCAAGGCTTGCGCTGTACCTGTGGGa                                                                                                    <  1:4176658/150‑1 (MQ=21)
                                                   gTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGa                                                                                                    <  1:4369042/150‑1 (MQ=255)
                                                    tgatgaGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAg                                                                                                   >  2:1428566/1‑150 (MQ=255)
                                                      atgaGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACACGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCg                                                                                                 >  1:5517535/1‑150 (MQ=255)
                                                       tgaGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTGTCTCGCAAGCTTCAAGGCTTGCGCTGTACCTGTGGGAGCga                                                                                                <  1:4037925/150‑2 (MQ=21)
                                                        gagaCTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCCTGAAGCTTGCGAGACACCTGTGGGAGCggg                                                                                               <  2:2010442/150‑1 (MQ=255)
                                                         agaCTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCACGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGc                                                                                              <  2:2949185/150‑1 (MQ=255)
                                                         agaCTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGc                                                                                              <  1:2545782/150‑1 (MQ=255)
                                                         agaCTCGCCACGGCTCTCAGACGGCACAAGAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCGTGAAGCTTGCGAGACACCTGTGGGAGCGGGc                                                                                              >  2:7105969/1‑150 (MQ=255)
                                                         agaCTCGCCACGGCTCTCAGACGGCACAACAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGc                                                                                              <  2:618231/150‑1 (MQ=255)
                                                           aCTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGGCCGTTCGCGGGCACGCCCGCTCCCACAGGTGTCTCGCAAGCTTCAAGGCTTGCGCTGTACCTGTGGGagcgagcgc                                                                                            <  1:1000586/150‑2 (MQ=14)
                                                           aCTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGGCCGTTCGCGGGCACGCCCGCTCCCACAGGTGTCTCGCAAGCTTCAAGGCTTGCGCTGTACCTGTGGGagcgagcgc                                                                                            <  1:1868919/150‑2 (MQ=14)
                                                             tCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGc                                                                                          <  1:6739233/150‑1 (MQ=255)
                                                             tCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGc                                                                                          >  1:1817756/1‑150 (MQ=255)
                                                             gctgggtcgcggcggaccgcaccattcgcgccgcggcccccgagatctggaaggcattggattgaggtgtgtccaTTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGACGCTTGCGAGACACCTGTGGGAGCGGGCGTGc                                                                                          >  2:7364243‑M1/76‑150 (MQ=255)
                                                                ccACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAgcgc                                                                                                                                   <  2:3272066/106‑1 (MQ=255)
                                                                ccACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCcg                                                                                       >  1:1576048/1‑150 (MQ=255)
                                                                 cACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTg                                                                                                                           >  2:6697482/1‑113 (MQ=255)
                                                                 cACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTg                                                                                                                           <  1:6697482/113‑1 (MQ=255)
                                                                   cGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCACGCCTTGACGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGa                                                                                    <  2:3148319/150‑1 (MQ=255)
                                                                         tcAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTGTCTCGCAAGCTTCAAGGCTTGCGCTGTACCTGTgg                                                                                                      <  1:45107/126‑1 (MQ=18)
                                                                         tcAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGcc                                                                              >  1:5267311/1‑150 (MQ=255)
                                                                            gACGGCACATTTGGTCGGTAGTTCAGCCTCCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGcctca                                                                           <  2:3720356/150‑1 (MQ=255)
                                                                             ccgccccattcgcgaatagcccaccgagatctggaaggctttgttttgcggtgtggccaTTGGCGGGCGCGCTCGCTCCCACAGGTACCGCGCACGCCTTGAAGCTGGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGcctcac                                                                          >  2:2530908‑M1/60‑150 (MQ=255)
                                                                                        ggTCGGTAGTTCCGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCAGTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAAt                                                               <  2:2256424/150‑1 (MQ=255)
                                                                                          tCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTTCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGc                                                             >  2:2844331/1‑150 (MQ=255)
                                                                                          tCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGc                                                             >  2:2736663/1‑150 (MQ=255)
                                                                                              tAGTTCAGCCTTCGCCTGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACTGGTACCGCGCAAGCCTTGACGCTTGCGAGCCACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTc                                                         <  2:788757/150‑1 (MQ=255)
                                                                                                    aGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGAtc                                                   >  1:6174485/1‑150 (MQ=255)
                                                                                                      ccTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGAtctc                                                 >  1:2005316/1‑150 (MQ=255)
                                                                                                        ttCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCgg                                               >  1:6634782/1‑150 (MQ=255)
                                                                                                        ttCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCgg                                               >  2:2213356/1‑150 (MQ=255)
                                                                                                         tCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGt                                              >  1:4396183/1‑150 (MQ=255)
                                                                                                            ccAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGc                                           <  2:4396183/150‑1 (MQ=255)
                                                                                                             cAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCACTGCCTTCCAGATCTCGGGGGCa                                          >  2:6401356/1‑150 (MQ=255)
                                                                                                             cAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGcc                                          >  2:4064859/1‑149 (MQ=255)
                                                                                                                   ggtggtCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTcg                                    >  1:6437744/1‑150 (MQ=255)
                                                                                                                       gtCTGGTTGTGTGTCCAGTCGCGGGCGCGCTCGCTCCCACAGGGACAGCGCACGCCTGGAAGCTGGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGaa                                <  2:6877671/150‑1 (MQ=39)
                                                                                                                        gattgaggtgaggccgTTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGacac                                                                                                            <  2:6137642‑M1/57‑1 (MQ=255)
                                                                                                                            gctgtgtgTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTg                           >  2:6064262/3‑150 (MQ=255)
                                                                                                                            gTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTg                           <  2:6174025/150‑1 (MQ=255)
                                                                                                                                tgtgtCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCGCGGTGGCAGAGGCGCGAGGGGTGCGGt                       >  2:3843932/1‑150 (MQ=38)
                                                                                                                                 gtgtCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTc                      <  1:2330271/150‑1 (MQ=255)
                                                                                                                                    tCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGGATAGTCGAGAATGGTGCGGTCcga                   >  2:7037180/1‑150 (MQ=255)
                                                                                                                                    tCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCcga                   >  2:575627/1‑150 (MQ=255)
                                                                                                                                    tCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCcga                   >  1:3282196/1‑150 (MQ=255)
                                                                                                                                    gccATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGACCGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCCTCTTCGCGAATGGTGCGGTCcgc                   >  2:1622248/2‑149 (MQ=255)
                                                                                                                                     ccaggCGCGGGCGCGCTCGCCCCCCCCGGTACCGCGCACGCCTTGACGCTTGCGAGACACCTGGGGGCGCGGGCGGGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCcgac                  >  1:4752032‑M1/6‑150 (MQ=255)
                                                                                                                                      cATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGTTCcgacg                 >  1:3430174/1‑150 (MQ=255)
                                                                                                                                          CGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGaa             >  1:6064262‑M1/1‑150 (MQ=255)
                                                                                                                                          CGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGaa             >  2:1368385‑M1/1‑150 (MQ=255)
                                                                                                                                               GCGCGCGCGCGCCCACCGGGACAGCGCACGCTTTGACGGTTGCGCGACACGGGTGGGGGGGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGc        >  1:7105969‑M1/1‑150 (MQ=255)
                                                                                                                                                     tcgctcCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATcc  >  2:5924465/1‑150 (MQ=255)
                                                                                                                                                     tcgctcCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGGCCGACGAGAACCAGCCCATcc  >  1:649616/1‑150 (MQ=255)
                                                                                                                                                     tcgctcCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGACTTCCAGATCTCGGTGGAATATTATCGAATGGTGCGGTCCGACGAGAACCAGCCCATcc  >  2:1295320/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGGGGGGATGAAAAACCCACGCAGTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCC  >  NC_002947/5748191‑5748489

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: