Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,743,394 T→C 8.0% 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,743,3940TC 8.0% 168.8 / 4.7 65intergenic (‑227/+20)fbp/glgPclass 1 fructose‑1,6‑bisphosphatase/glycogen phosphorylase
Reads supporting (aligned to +/- strand):  ref base T (29/28);  new base C (3/2);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.10e-01

TCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTGCGCCACCACTCCTGCGGTGTGTGCCACAGATCCTC  >  NC_002947/5743247‑5743542
                                                                                                                                                   |                                                                                                                                                    
tCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTg                                                                                                                                                     >  7:210765/1‑149 (MQ=255)
    cACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCcg                                                                                                                                                 <  6:151984/149‑1 (MQ=255)
      cGGCTCTCAGACGGCACAATAGGGCGGTAGTGCAGCCTTCGCCAGGCCGGTGGTCTGGTTCTGTGTCCATTCGCGGGCGCGCTCGCTCCCACCGGTACAGCGCAAGCCTGGAAGCTTGCGAGACACCTGTCGGAGCGGGCACGCCcgct                                                                                                                                               <  5:106690/149‑2 (MQ=255)
        gCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCGCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGaa                                                                                                                                             <  5:180102/149‑1 (MQ=255)
         ctctCAGACGGCAAAATAGGGCGGTTGTTCAGCCTTCGCCAGGCGGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGTTACAGCGCAAGCTTTTAAGCTTGCTAGACACGTGTGGGTGGGGGGGGGCCCGCGAAc                                                                                                                                            <  8:121739/149‑1 (MQ=255)
                     acaATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGcctcacctca                                                                                                                                <  8:210765/149‑1 (MQ=255)
                     acaATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGcctcacctc                                                                                                                                 >  4:81152/1‑148 (MQ=255)
                     acaATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGcctcacct                                                                                                                                  <  4:88879/147‑1 (MQ=255)
                     acaATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGcctcacct                                                                                                                                  >  3:88879/1‑147 (MQ=255)
                      caATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCaa                                                                                                                               <  3:81152/149‑1 (MQ=255)
                        aTAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATc                                                                                                                             <  4:70013/149‑1 (MQ=255)
                                gTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCt                                                                                                                     >  3:343136/1‑149 (MQ=255)
                                    ttCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGGTCGCTCCCACAGGAACACCGAAAGCCTTGACGCTTGCTAGACACCTGTGGGAGCGGGCGGGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCa                                                                                                                 <  5:213610/149‑1 (MQ=255)
                                    ttCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCa                                                                                                                 <  5:155940/149‑1 (MQ=255)
                                     tCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTc                                                                                                                   >  8:69072/1‑146 (MQ=255)
                                      cAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGCCACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTc                                                                                                                   >  3:82542/1‑145 (MQ=255)
                                          cTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCt                                                                                                                     >  2:283287/1‑139 (MQ=255)
                                          cTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCt                                                                                                                     <  1:283287/139‑1 (MQ=255)
                                                 aGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCa                                                                                                    >  4:239463/1‑149 (MQ=255)
                                                        tggtGTGGTTGTGTGGCATTTGGGGGGCGGGCTCTCTCCCAGAGGTACCGCGCATGCCTTGATGCTTGGGAGTCACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTcgc                                                                                             <  1:265151/149‑1 (MQ=255)
                                                                ttGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTg                                                                                     <  7:69072/149‑1 (MQ=255)
                                                                     tgtCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGt                                                                                                                                                      >  8:110900/1‑79 (MQ=255)
                                                                      gtCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCCCCTCAAGCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTcc                                                                               >  5:23851/1‑149 (MQ=255)
                                                                         cATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCcga                                                                             <  7:114609/148‑1 (MQ=255)
                                                                              gcgGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTg                                                                                     >  8:179484/1‑135 (MQ=255)
                                                                              gcgGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTg                                                                                     <  7:179484/135‑1 (MQ=255)
                                                                                      gctcgctcCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCa                                                               >  5:10424/1‑149 (MQ=255)
                                                                                      gcccgctcccacaggtgtctcgcaagcttcaaggcttgcgctgtacctgtgggagcgGGCGTGCCCGCGAACGGCCTCACCTCAATCCAAt                                                                                                                         >  8:294192‑M2/58‑91 (MQ=255)
                                                                                      gcccgctcccacaggtgtctcgcaagcttcaaggcttgcgctgtacctgtgggagcgGGCGTGCCCGCGAACGGCCTCACCTCAATCCAAt                                                                                                                         <  7:294192‑M2/34‑1 (MQ=255)
                                                                                         CGCTCCCACAGGTCCAGCGCAAGCCTTGACGCTTGCGTGACCCCTGTGGGAGCGGGCGGGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATcc                                                            >  5:195300‑M1/1‑149 (MQ=255)
                                                                                          gcgcccACTGGTCCCGCGCAAGCCGGGAAGCTGGCGAGCCCCCTGTGGTCGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATcc                                                            <  5:190348/145‑1 (MQ=255)
                                                                                          GCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGAGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATccg                                                           <  1:5623‑M1/149‑1 (MQ=255)
                                                                                            tcccacaggtgtctcgcaagcttcaaggcttgcgctgtacctgtgggagcgAGCGCGCCCGCGAATGGACACACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATccgtg                                                         >  2:332878‑M2/52‑149 (MQ=255)
                                                                                             CCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGt                                                                                 <  5:306923‑M1/124‑1 (MQ=255)
                                                                                             CCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGt                                                                                 >  6:306923‑M1/1‑124 (MQ=255)
                                                                                              ccacaggtgtctcgcaagcttcaaggcttgcgctgtacctgtgggagcgGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATccgtgcc                                                       >  7:309913‑M2/50‑149 (MQ=255)
                                                                                              ccacaggtgtctcgcaagcttcaaggcttgcgctgtacctgtgggagcgAGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATccgtgcc                                                       >  3:234445‑M2/50‑149 (MQ=255)
                                                                                               CACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATccgtgccg                                                      <  8:299771‑M1/149‑1 (MQ=255)
                                                                                                acaggtgtctcgcaagcttcaaggcttgcgctgtacctgtgggagcgAGCGCGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATccgtgccgt                                                     >  1:110659‑M2/48‑149 (MQ=255)
                                                                                                   GGTCCAGCGCAAGCCTTCACGCTGGCGAGACACCTGGGGGCGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGtt                                                  >  6:23851‑M1/1‑149 (MQ=255)
                                                                                                    gtgtctcgcaagcttcaaggcttgcgctgtacctgtgggagcgGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTc                                                 <  7:37569‑M2/106‑1 (MQ=255)
                                                                                                    gtgtctcgcaagcttcaaggcttgcgctgtacctgtgggagcgAGCGCGCCCGCGAATGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGtt                                                  >  6:153578‑M2/44‑148 (MQ=255)
                                                                                                       tctcgcaagcttcaaggcttgcgctgtacctgtgggagcgGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCa                                                                                                    >  6:100099‑M2/41‑95 (MQ=255)
                                                                                                       tctcgcaagcttcaaggcttgcgctgtacctgtgggagcgGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCa                                                                                                    <  5:100099‑M2/55‑1 (MQ=255)
                                                                                                       tctcgcaagcttcaaggcttgcgctgtacctgtgggagcgAGCGCGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAAt                                              <  2:110659‑M2/109‑1 (MQ=255)
                                                                                                          cgcAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATAcc                                           >  7:5650/1‑149 (MQ=255)
                                                                                                          cgcAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGCTGCGGTCCGACGAGAACCAGCCCATCCGTGCCCTGTTCAATAcc                                           >  4:47027/1‑149 (MQ=255)
                                                                                                             aaGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGcc                                        >  8:252169/1‑149 (MQ=255)
                                                                                                             aaGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAAACAGCCCATCCGTGACGTGTTCAATCCAGcc                                        >  8:257093/1‑149 (MQ=255)
                                                                                                                 cTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTgc                                    >  4:302326/1‑149 (MQ=255)
                                                                                                                 cTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTgc                                    >  5:103286/1‑149 (MQ=255)
                                                                                                                 agtGATGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTgc                                    <  4:82542/147‑1 (MQ=255)
                                                                                                                   tGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTgcgc                                  >  1:213275/1‑149 (MQ=255)
                                                                                                                      aGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTGCGccac                               <  7:257093/149‑1 (MQ=255)
                                                                                                                      aGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTGCGccac                               <  7:252169/149‑1 (MQ=255)
                                                                                                                            cGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCa                                                                    >  2:177249/1‑106 (MQ=255)
                                                                                                                            cGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCa                                                                    <  1:177249/106‑1 (MQ=255)
                                                                                                                            cGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTGCGCCACCACTcc                         <  2:213275/149‑1 (MQ=255)
                                                                                                                            cGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTGCGCCACCACTc                          <  8:5650/148‑1 (MQ=255)
                                                                                                                             gagaCACCTGTGGGAGCGGGCGTGACCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATAGTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTGCGACACCACTCCt                        >  5:348338/1‑149 (MQ=255)
                                                                                                                               gaCACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTGCGCCACCACTCCTGc                      >  8:330733/1‑149 (MQ=255)
                                                                                                                                   cctgtgggagcgGGCGTGCCCGCGAACGGCCTCACCTCAATCCAAAGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTACGACGAGAACCAGCACATCCGTGCCGTGTTCAATACCGCCTTGCGCCACCACTCCTGCGgtg                  >  4:223601‑M2/13‑149 (MQ=255)
                                                                                                                                    ctgtgggagcgGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTGCGCCACCACTCCTGCGgtgt                 <  6:10424‑M2/138‑1 (MQ=255)
                                                                                                                                        gggagcgGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTGCGCCACCACTCCTGCGGTTTGTGc             <  7:330733‑M2/142‑1 (MQ=255)
                                                                                                                                          gagcgGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTGCGCCACCACTCCTGCGGTGTGTGCca           >  3:245632‑M2/6‑149 (MQ=255)
                                                                                                                                                 cGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTGCGCCACCACTCCTGCGGTGTGTGCCACAGATcc    <  3:239463/149‑1 (MQ=255)
                                                                                                                                                   tGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTGCGCCACCACTCCTGCGGTGTGTGCCACAGATCCTc  >  4:286819/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTGCGCCACCACTCCTGCGGTGTGTGCCACAGATCCTC  >  NC_002947/5743247‑5743542

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