Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,669,674 C→A 9.7% intergenic (+41/+34) prpD → / ← acnB 2‑methylcitrate dehydratase/aconitate hydratase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,669,6740CA 9.7% 186.7 / 8.3 62intergenic (+41/+34)prpD/acnB2‑methylcitrate dehydratase/aconitate hydratase B
Reads supporting (aligned to +/- strand):  ref base C (30/26);  new base A (2/4);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 4.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

CCGCTTTGCCCGGCAGCGGTGTCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCAC‑‑‑GGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTGTCGATGTCTTTC  >  NC_002947/2669528‑2669819
                                                                                                                                                     |                                                                                                                                                 
ccGCTTTGCCCGGCAGCGGTGTCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCAC‑‑‑GGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTg                                                                                                                                                 >  7:330562/1‑149 (MQ=255)
   cTTTGCCCGGCAGCGGTGTCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCAC‑‑‑GGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCg                                                                                                                                              >  5:228257/1‑149 (MQ=255)
         ccGGCAGCGGTGTCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCAC‑‑‑GGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCAGCGATGCCGGGGCtt                                                                                                                                        <  4:112821/149‑1 (MQ=255)
          cGGCAGCGGTGTCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCAC‑‑‑GGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCttt                                                                                                                                       >  5:373687/1‑149 (MQ=255)
          cGGCAGCGGTGTCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCAC‑‑‑GGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCAGCGATGCCGGGGCttt                                                                                                                                       <  7:47815/149‑1 (MQ=255)
             cAGCGGTGTCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCAC‑‑‑GGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGtt                                                                                                                                    <  8:150400/149‑1 (MQ=255)
              aGCGGTGTCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCAC‑‑‑GGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGttt                                                                                                                                   <  6:373687/149‑1 (MQ=255)
              aGCGGTGTCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCAC‑‑‑GGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGttt                                                                                                                                   <  5:63838/149‑1 (MQ=255)
                  gtgtCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCAC‑‑‑GGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACaa                                                                                                                               <  5:324667/149‑1 (MQ=255)
                   tgtCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCAC‑‑‑GGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAAc                                                                                                                              <  8:17638/149‑1 (MQ=255)
                     tCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCCCTCTTAGTGTGCACAGGTTCGTGTACTTGTTGGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACa                                                                                                                                <  1:335163/148‑1 (MQ=255)
                       agcagGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCAC‑‑‑GGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTg                                                                                                                           >  2:112715/1‑148 (MQ=255)
                          agGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCAC‑‑‑GGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTa                                                                                                                                  <  3:368199/138‑1 (MQ=255)
                          agGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCAC‑‑‑GGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCCGCATCGCTGCCGGGGCTTTGTTTa                                                                                                                                  >  4:368199/1‑138 (MQ=255)
                            gTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCAC‑‑‑GGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTc                                                                                                                     <  6:228257/149‑1 (MQ=255)
                               ttCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCAC‑‑‑GGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAg                                                                                                                   <  2:148177/148‑1 (MQ=255)
                                tCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCAC‑‑‑GGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAAc                                                                                                                              <  3:303541/136‑1 (MQ=255)
                                tCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCAC‑‑‑GGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAAc                                                                                                                              >  4:303541/1‑136 (MQ=255)
                                tCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCAC‑‑‑GGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGc                                                                                                                  <  4:192968/148‑1 (MQ=255)
                                 cGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCAC‑‑‑GGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCt                                                                                                                 >  2:45857/1‑148 (MQ=255)
                                 cGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCAC‑‑‑GGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCt                                                                                                                 <  8:280538/148‑1 (MQ=255)
                                  gAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCAC‑‑‑GGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGc                                                                                                                                          >  7:61575/1‑122 (MQ=255)
                                  gAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCAC‑‑‑GGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGc                                                                                                                                          <  8:61575/122‑1 (MQ=255)
                                          cAAGGACCAGCAGAAGCTGGAAGGCAC‑‑‑GGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGaa                                                                                                       <  1:396600/149‑1 (MQ=255)
                                                 cagcagAAGCTGGAAGGCAC‑‑‑GGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCgg                                                                                                <  8:330562/149‑1 (MQ=255)
                                                  agcagAAGCTGGAAGGCAC‑‑‑GGCTGTGCACAGGATCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTATGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCgg                                                                                                >  1:122827/1‑148 (MQ=255)
                                                       aaGCTGGAAGGCAC‑‑‑GGCTGTGCACAGGTTCGTTGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCtt                                                                                          <  3:103382/149‑1 (MQ=255)
                                                        aGCTGGAAGGCAC‑‑‑GGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCACCTTACGCCTGAACCACCGGGATCTTg                                                                                         <  1:216778/149‑1 (MQ=255)
                                                             gAAGGCAC‑‑‑GGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGtt                                                                                    >  3:170897/1‑149 (MQ=255)
                                                                   aC‑‑‑GGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCAGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGc                                                                              >  2:336838/1‑149 (MQ=255)
                                                                        ggCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCg                                                                            <  3:275961/149‑1 (MQ=255)
                                                                                    gTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACt                                                                >  6:234382/1‑149 (MQ=255)
                                                                                           gACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGAt                                                         >  8:242244/1‑149 (MQ=255)
                                                                                                 tttgtGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATc                                                        <  4:383803/144‑1 (MQ=255)
                                                                                                 tttgtGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATc                                                        >  3:383803/1‑144 (MQ=255)
                                                                                                 tttgtGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGt                                                    <  1:45857/148‑1 (MQ=255)
                                                                                                           gAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAgg                                         >  1:67636/1‑149 (MQ=255)
                                                                                                           gAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCCCAgg                                         >  1:409099/1‑149 (MQ=255)
                                                                                                             aGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAgg                                         >  1:156212/1‑147 (MQ=255)
                                                                                                                         aGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTc                           >  8:334919/1‑149 (MQ=255)
                                                                                                                         aGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACATc                           >  1:163104/1‑149 (MQ=255)
                                                                                                                         aGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACCCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTc                           >  5:408950/1‑149 (MQ=255)
                                                                                                                           gCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCTCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCgg                         >  8:194555/1‑149 (MQ=255)
                                                                                                                           gCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCgg                         >  3:328091/1‑149 (MQ=255)
                                                                                                                           gCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCgg                         >  7:386005/1‑149 (MQ=255)
                                                                                                                               aaaCAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCAAAGCTCAGGTAGCGGTAAACGTCGGCAGc                     >  2:111632/1‑149 (MQ=255)
                                                                                                                                   aacaaAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCTGCGATCTGGTCGAAGCTCATGTAGCGGTAAACGTCGGCAGCCATg                 <  2:122827/149‑1 (MQ=255)
                                                                                                                                   aacaaAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTATCGGTAAACGTAGGCAGCCAt                  >  6:285407/1‑148 (MQ=255)
                                                                                                                                     caaAGCCCCGGCAGCGATGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCt               >  5:264609/1‑149 (MQ=255)
                                                                                                                                      aaaGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTg              >  7:249778/1‑149 (MQ=255)
                                                                                                                                      aaaGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTg              >  2:298022/1‑149 (MQ=255)
                                                                                                                                      aaaGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTg              >  4:107351/1‑149 (MQ=255)
                                                                                                                                       aaGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTGt             <  2:156212/149‑1 (MQ=255)
                                                                                                                                        aGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTg              >  2:365282/1‑147 (MQ=255)
                                                                                                                                          ccccGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTGTCGa          <  7:242244/149‑1 (MQ=255)
                                                                                                                                          ccccGGCAGCGATGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTGTCGa          <  6:264609/149‑1 (MQ=255)
                                                                                                                                           cccGGCAGCGATGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTGTc            >  7:298675/1‑146 (MQ=255)
                                                                                                                                           cccGGCAGCGATGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTGTc            <  8:298675/146‑1 (MQ=255)
                                                                                                                                             cGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCCCAGGTAGCGGTAAACGTCGGCAGCCATGCTGTCGATGt       <  2:409099/149‑1 (MQ=255)
                                                                                                                                               gCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTGTCGATGTCt     >  4:146520/1‑149 (MQ=255)
                                                                                                                                                  tCGCTGCCGGTGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTGTCGATGTCTTTc  <  3:107351/149‑1 (MQ=255)
                                                                                                                                                  tCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTGTCGATGTCTTTc  <  5:234382/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
CCGCTTTGCCCGGCAGCGGTGTCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCAC‑‑‑GGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTGTCGATGTCTTTC  >  NC_002947/2669528‑2669819

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