Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 4,740,811 (T)7→6 intergenic (+48/+705) gltA → / ← yeiW citrate synthase/putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,740,8050T.92.5% 196.0 / 13.8 54intergenic (+42/+711)gltA/yeiWcitrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (0/4);  new base . (24/25);  total (24/30)
Fisher's exact test for biased strand distribution p-value = 1.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTT  >  NC_002947/4740659‑4740946
                                                                                                                                                  |                                                                                                                                             
cACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttt                                                                                                                                             <  1:285139/149‑1 (MQ=255)
 aCTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCtttt                                                                                                                                            >  6:11827/1‑149 (MQ=255)
 aCTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCtttt                                                                                                                                            <  8:7720/149‑1 (MQ=255)
  cTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGACGCGCGACATCGTTGCCCTGAAAGACCGCCAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCGAGGGCAGCcttctt                                                                                                                                           <  3:5789/149‑1 (MQ=255)
  cTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttt                                                                                                                                           <  1:393332/149‑1 (MQ=255)
  cTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAAAAGGCTGCCCGAGGGCAGCcttctt                                                                                                                                           <  1:239744/149‑1 (MQ=255)
  cTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAAAAGGCTGCCCGAGGGCAGCcttctt                                                                                                                                           >  2:239744/1‑149 (MQ=255)
   tGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttt                                                                                                                                           >  8:138878/1‑148 (MQ=255)
      aaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttatt                                                                                                                                       <  1:270601/149‑4 (MQ=255)
       aagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttattt                                                                                                                                      <  2:251413/149‑5 (MQ=255)
                    cAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCg                                                                                                                          <  8:263903/147‑1 (MQ=255)
                      gCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCGAGGGCAGCcttctt                                                                                                                                           >  2:49606/1‑129 (MQ=255)
                      ccccGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCGAGGGCAGCcttctt                                                                                                                                           <  1:49606/128‑1 (MQ=255)
                          gTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCcact                                                                                                                    >  2:179401/1‑147 (MQ=255)
                            aCAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGACCTGTGCCGGCCTCTTcgc                                                                                                                >  1:128698/1‑149 (MQ=255)
                              aaGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCgg                                                                                                              >  8:370245/1‑149 (MQ=255)
                              aaGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTAGCgg                                                                                                              >  6:29869/1‑149 (MQ=255)
                                gATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGcc                                                                                                            <  2:324611/149‑1 (MQ=255)
                                    ggCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGc                                                                                                        <  5:11827/149‑1 (MQ=255)
                                     gCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGcc                                                                                                       >  8:208917/1‑149 (MQ=255)
                                       cGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAACAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCACTTCGCGGCCTTGCCcg                                                                                                     >  5:299895/1‑149 (MQ=255)
                                           cgcgCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGa                                                                                                  <  7:138878/148‑1 (MQ=255)
                                             cgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAagag                                                                                              <  7:345587/149‑1 (MQ=38)
                                             cgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                               <  7:370245/149‑1 (MQ=38)
                                             cgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                               <  4:344595/149‑1 (MQ=38)
                                             cgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                               <  7:208917/149‑1 (MQ=255)
                                                  gCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCTTGCCcgcg                                                                                                   >  2:199819/1‑140 (MQ=255)
                                                   cTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAAGAGGCCgg                                                                                         <  5:19442/149‑1 (MQ=38)
                                                   cTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTCAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCgg                                                                                         <  7:368193/149‑1 (MQ=255)
                                                         acCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCCCGGCCTTGCCCGCGAAGCGggcggctccga                                                                                   >  7:232047/1‑144 (MQ=255)
                                                          ccGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGac                                                                                  >  2:301826/1‑149 (MQ=255)
                                                           cGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGaca                                                                                 >  6:278558/1‑149 (MQ=255)
                                                           cGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCAGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCATCGCGGCCTTGCCCGCGAAGAGTCAGGCCAAGaca                                                                                 >  1:92281/1‑149 (MQ=255)
                                                               gAGCAGATGCGCTACATCGTTCCCCTCAAAGACCGCAATTAAGCTCCAGTAGCCCAACGAAAGACGGCGGCCCGCGGGCAGCTTTTTTTATTTTTCCGGTGCCGCCGTCTTCGCGCCCTTGCCCGCTAAGAGGCCGGCACAGAcaccac                                                                             <  4:369279/149‑1 (MQ=255)
                                                                   tgatgcgcGACATCGTTGCACTTAACTACCGCAATTAAGCTGCAGTTGCCGAACGCAATACGCCTGCCCTCGGGCAGCCTTTTTTTTTTGTCCTGTGCGGGCCTCTTCGCGGCCTTGCCCGCGAAGTGGCCGGCACAGACACCACaaaa                                                                         <  6:299895/145‑1 (MQ=255)
                                                                        cgcgACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGTCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCa                                                                    <  1:102122/149‑1 (MQ=255)
                                                                           gACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGtt                                                                 <  1:83264/149‑1 (MQ=255)
                                                                                ctttGCCCTGAAAGCCCCAAAGTAACCTGCAGTAGCCGACCCAAACACGGCTGCCCTCGGCCAGCCTTTTTTATTTGTCCTGTCCCGGCCTCTTCGCGGCCTTGCCCGCTAAGTGGCCGGCACAGACACCACaaaa                                                                         <  4:331135/134‑1 (MQ=255)
                                                                                cGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTcgcg                                                                                                               >  7:114780/1‑98 (MQ=255)
                                                                                cGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTcgcg                                                                                                               <  8:114780/98‑1 (MQ=255)
                                                                                 gTTCCCCTCAAACCCCCCAAGAAAGCGGCAGTCGCCCACCGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGGGCCGCCCTCTTCGCGGCCTTCCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTtc                                                             <  8:87677/147‑1 (MQ=255)
                                                                                   tGCCATGAAATACCGAAAGTAAGCTGAAGTATCCGAACGCATGTTGTCTTCCCTCGGGCAGCATTTTTTATTTGTCCTGTGCCGGCTTTTTGGCGGCCTTTCCCGGTATGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTcc                                                          <  2:92281/148‑1 (MQ=255)
                                                                                     cccTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGttt                                                                <  6:185150/140‑1 (MQ=255)
                                                                                     cccTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGttt                                                                >  5:185150/1‑140 (MQ=255)
                                                                                     cccTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTccgc                                                        <  1:301826/148‑1 (MQ=255)
                                                                                      ccTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCCGACTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTATTTCTccgcac                                                      >  2:61649/1‑147 (MQ=255)
                                                                                          aaaGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCATGCCCGCGAAGCGGCCGGCACAGACCCCCCAAAACCTCAGTGTTTCTccgcccgcg                                                   >  7:5551/1‑147 (MQ=255)
                                                                                             gACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCTCTCACg                                               <  5:278558/149‑1 (MQ=255)
                                                                                              aCCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCCCAGCCACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACg                                               >  5:87067/1‑148 (MQ=255)
                                                                                                cGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAg                                            >  3:324283/1‑149 (MQ=255)
                                                                                                  cAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGcc                                          >  2:77372/1‑149 (MQ=255)
                                                                                                     gTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCGCAGGTTTTCTCCGCCCGCTCACGCAGCCCtt                                       >  5:138792/1‑149 (MQ=255)
                                                                                                        aGCTGCAGTAGCGGAACGCACGACGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCGGTGCCGGCCTCTTCGGGGCCTTGCCCGCTAAGTGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTc                                     <  8:5551/148‑1 (MQ=255)
                                                                                                          cTGCAGTAGCGGAACGCACCACGGCGGCCGTCGGGCAGCCTTTTTTTTTTGGCCTGTGGCGGCCTCTTCGCGGCCTTGCCCGCTAATAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                                  <  8:232047/149‑1 (MQ=255)
                                                                                                          cTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                                  <  1:77372/149‑1 (MQ=255)
                                                                                                                      aaCGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTATTCGCCGCCTTGCCACCGAAGAAGCCGGCACACACACCAAAAAAACAAAGTTTTTCTACACCACATCACCCAGCACTTTAAGGGTAATGaaagca                      >  8:22360/1‑145 (MQ=255)
                                                                                                                           aagaagGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACCCCACAAAACCGCAGTTTTTCTCCGCCCGCGCACGCAGCCCGTTCAGGGTAAATAACGGCGCATc                 >  1:332722/1‑149 (MQ=255)
                                                                                                                           aagaagGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCa                                                                    <  5:88727/98‑1 (MQ=37)
                                                                                                                           aagaagGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCa                                                                    >  6:88727/1‑98 (MQ=37)
                                                                                                                                ggCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAgg                                                                                             <  7:328721/68‑1 (MQ=25)
                                                                                                                                 tcTGCCCTCGGGAAGCCTTTTTTATTTGTGCTGTGCCGGCCTCTTCGCGGCATTGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCcaccac           <  7:22360/148‑1 (MQ=255)
                                                                                                                                    gcccgggggCTGCCTTTTTTTTTTTACCTGTCCCGGCCTCTTCGGGGCCTTGCCCGCTACGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCa                                    <  5:272657/115‑1 (MQ=255)
                                                                                                                                     cccTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGaa        <  1:61649/148‑1 (MQ=255)
                                                                                                                                     cccTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTGCGCGGCCTGCCCCGCGAAGCGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTccgccc                                                      <  3:78698/102‑1 (MQ=255)
                                                                                                                                        tCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGGGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTAAGGGTATTGAACGGCGCATCCCCCACGAACttg    >  2:127766/1‑149 (MQ=255)
                                                                                                                                          gggCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCCCGGCCTTGCCCGCGAAGAGGCCGCCACAAACACCACAAAACCAAAGAATTACTACCCCCCCTAACGCAGCCCATTCAGCGTAATGAACCGCCCAAACACCACCAAAatgtt  >  2:122159/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
CACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTT  >  NC_002947/4740659‑4740946

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: