Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,646,886 G→A 7.2% S157S (AGC→AGT)  gcd ← quinoprotein glucose dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,646,8860GA 7.2% 164.4 / 2.1 55S157S (AGC→AGT) gcdquinoprotein glucose dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (25/26);  new base A (2/2);  total (27/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

CACCGGTGCGAATGCGCCAGGCTTCTTCCAGACGGTAGGCGTTCTGCGGGGTGATCTGGCGCAGCGGCGAGTAGCGGTCGCCATGCTCGGTACGGCCGTAGGCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAACTGGCTGGCCAGCGCGCAAGCGCCCGAAGCGACCACTGCCACGCCGAGCAGTGCAGTGTTGGCTTTGCTGGCTGGGCCGATCAGCGGGCGACGTGC  >  NC_002947/1646738‑1647029
                                                                                                                                                    |                                                                                                                                               
cACCGGTGCGAATGCGCCAGGCTTCTTCCAGACGGTAGGCGTTCTGCGGGGTGATCTGGCGCAGCGGCGAGTAGCGGTCGCCATGCTCGGTACGGCCGTAGGCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCg                                                                                                                                                 >  4:31954/1‑149 (MQ=255)
  ccGGTGCGAATGCGCCAGGCTTCTTCCAGACGGTAGGCGTTCTGCGGGGTGATCTGGCGCAGCGGCGAGTAGCGGTCGCCATGCTCGGTACGGCCGTAGGCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCgct                                                                                                                                               >  1:133777/1‑149 (MQ=255)
  ccGGTGCGAATGCGCCAGGCTTCTTCCAGACGGTAGGCGTTCTGCGGGGTGATCTGGCGCAGCGGCGAGTAGCGGTCGCCATGCTCGGTACGGCCGTAGGCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCgct                                                                                                                                               <  3:31954/149‑1 (MQ=255)
     gTGCGAATGCGCCAGGCTTCTTCCAGACGGTAGGCGTTCTGCGGGGTGATCTGGCGCAGCGGCGAGTAGCGGTCGCCATGCTCGGTACGGCCGTAGGCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCgctgct                                                                                                                                            >  6:103689/1‑149 (MQ=255)
         gAATGCGCCAGGCTTCTTCCAGACGGTAGGCGTTCTGCGGGGTGATCTGGCGCAGCGGCGAGTAGCGGTCGCCATGCTCGGTACGGCCGTAGGCCTGCCATTCGCCGTCGGGCATGTCCGGGGCGGCGCTGGCCATTTCGCTGCTGTcg                                                                                                                                        >  6:264325/1‑149 (MQ=255)
          aaTGCGCCAGGCTTCTTCCAGACGGTAGGCGTTCTGCGGGGTGATCTGGCGCAGCGGCGAGTAGCGGTCGCCATGCTCGGTACGGCCGTAGGCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTcgc                                                                                                                                       <  5:24472/149‑1 (MQ=255)
                 cAGGCTTCTTCCAGACGGTAGGCGTTCTGCGGGGTGATCTGGCGCAGCGGCGAGTAGCGGTCGCCATGCTCGGTACGGCCGTAGGCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAg                                                                                                                                >  4:283201/1‑149 (MQ=255)
                 cAGGCTTCTTCCAGACGGTAGGCGTTCTGCGGGGTGATCTGGCGCAGCGGCGAGTAGCGGTCGCCATGCTCGGTACGGCCGTAGGCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAg                                                                                                                                >  2:214898/1‑149 (MQ=255)
                   ggCTTCTTCCAGACGGTAGGCGTTCTGCGGGGTGATCTGGCGCAGCGGCGAGTAGCGGTCGCCATGCTCGGTACGGCCGTAGGCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGtt                                                                                                                              >  6:331412/1‑149 (MQ=255)
                     cttcttCCAGACGGTAGGCGTTCTGCGGGGTGATCTGGCGCAGCGGCGAGTAGCGGTCGCCATGCTCGGTACGGCCGTAGGCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCg                                                                                                                            >  7:760/1‑149 (MQ=255)
                             aGACGGTAGGCGTTCTGCGGGGTGATCTGGCGCAGCGGCGAGTAGCGGTCGCCATGCTCGGTACGGCCGTAGGCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAacac                                                                                                                    >  3:91891/1‑149 (MQ=255)
                                    aGGCGTTCTGCGGGGTGATCTGGCGCAGCGGCGAGTAGCGGTCGCCATGCTCGGTACGGCCGTAGGCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAg                                                                                                                                >  8:107894/1‑130 (MQ=255)
                                    aGGCGTTCTGCGGGGTGATCTGGCGCAGCGGCGAGTAGCGGTCGCCATGCTCGGTACGGCCGTAGGCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAg                                                                                                                                <  7:107894/130‑1 (MQ=255)
                                            tGCGGGGTGATCTGGCGCAGCGGCGAGTAGCGGTCGCCATGCTCGGTACGGCCGTAGGCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACACGGATGGGt                                                                                                     >  6:106642/1‑149 (MQ=255)
                                                gggTGATCTGGCGCAGCGGCGAGTAGCGGTCGCCATGCTCGGTACGGCCGTAGGCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAAc                                                                                                 >  2:259631/1‑149 (MQ=255)
                                                  gTGATCTGGCGCAGCGGCGAGTAGCGGTCGCCATGCTCGGTACGGCCGTAGGCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAActg                                                                                               >  8:242712/1‑149 (MQ=255)
                                                     aTCTGGCGCAGCGGCGAGTAGCGGTCGCCATGCTCGGTACGGCCGTAGGCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAActggct                                                                                            >  4:310164/1‑149 (MQ=255)
                                                     aTCTGGCGCAGCGGCGAGTAGCGGTCGCCATGCTCGGTACGGCCGTAGGCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAActggc                                                                                             <  8:760/148‑1 (MQ=255)
                                                            gcAGCGGCGAGTAGCGGTCGCCATGCTCGGTACGGCCGTAGGCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAACTGGCTGGCCAgc                                                                                     <  5:103689/149‑1 (MQ=255)
                                                                 ggcgAGTAGCGGTCGCCATGCTCGGTACGGCCGTAGGCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAACTGGCTGGCCAgcgcgc                                                                                 >  6:354241/1‑148 (MQ=255)
                                                                  gcgAGTAGCGGTCGCCATGCTCGGTACGGCCGTAGGCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAACTGGCTGGCCAGCGCGCaa                                                                               <  3:310164/149‑1 (MQ=255)
                                                                  gcgAGTAGCGGTCGCCATGCTCGGTACGGCCGTAGGCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAACTGGCTGGCCAGCGCGCaa                                                                               <  1:214898/149‑1 (MQ=255)
                                                                      gTAGCGGTCGCCATGCTCGGTACGGCCGTAGGCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAACTGGCTGGCCAGCGCGCAAgcgc                                                                           <  1:259631/149‑1 (MQ=255)
                                                                         gCGGTCGCCATGCTCGGTACGGCCGTAGGCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCgg                                                                                                            >  7:76849/1‑113 (MQ=255)
                                                                         gCGGTCGCCATGCTCGGTACGGCCGTAGGCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCgg                                                                                                            <  8:76849/113‑1 (MQ=255)
                                                                            gTCGCCATGCTCGGTACGGCCGTAGGCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAActggct                                                                                            >  4:137104/1‑126 (MQ=255)
                                                                            gTCGCCATGCTCGGTACGGCCGTAGGCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAActggct                                                                                            <  3:137104/126‑1 (MQ=255)
                                                                               gCCATGCTCGGTACGGCCGTAGGCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAACTGGCTGGTCAGCGCGCAAGCGCCCGAAGCGa                                                                  <  1:46965/149‑1 (MQ=255)
                                                                                ccATGCTCGGTACGGCCGTAGGCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAACTGGCTGGCCAGCGCGCAAGCGCCCGAAGCGAc                                                                 <  5:354241/149‑1 (MQ=255)
                                                                                ccATGCTCGGTACGGCCGTAGGCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAACTGGCTGGCCAGCGCGCAAGCGCCCGAAGCGAc                                                                 <  4:71541/149‑1 (MQ=255)
                                                                                 cATGCTCGGTACGGCCGTAGGCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAACTGGCTGGCCAGCGCGCAAGCGCCCGAAGCGAct                                                                <  4:91891/149‑2 (MQ=255)
                                                                                     cTCGGTACGGCCGTAGGCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAACTGGCTGGCCAgc                                                                                     <  2:24537/124‑1 (MQ=255)
                                                                                     cTCGGTACGGCCGTAGGCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAACTGGCTGGCCAgc                                                                                     >  1:24537/1‑124 (MQ=255)
                                                                                     cTCGGTACGGCCGTAGGCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAACTGGCTGGCCAGCGCGCAAGCGCCCGAAGCGACCACTg                                                            <  7:242712/149‑1 (MQ=255)
                                                                                      tCGGTACGGCCGTAGGCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAACTGGCTGGCCAGCGCGCAAGCGCCCGAAGCGACCACTGc                                                           <  5:106642/149‑1 (MQ=255)
                                                                                         gTACGGCCGTAGGCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCACTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAACTGGCTGGCCAGCGCGCAAGCGCCCGAAGCGACCACTGCCAc                                                        >  7:85818/1‑149 (MQ=255)
                                                                                          tACGGCCGTAGGCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGGCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAACTGGCTGGCCAGCGCGCAAGCGCCCGAAGCGACCACTGCCACg                                                       <  8:214261/149‑1 (MQ=255)
                                                                                            cGGCCGTAGGCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAACTGGCTGGCCAGCGCGCAAGCGCCCGAAGCGACCACTGCCACGcc                                                     <  4:363634/149‑1 (MQ=255)
                                                                                              gCCGTAGGCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAACTGGCTGGCCAGCGCGCAAGCGCCCGAAGCGACCACTGCCACGCCGa                                                   <  5:264325/149‑1 (MQ=255)
                                                                                               ccGTAGGCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCACTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAACTGGCTGGCCAGCGCGCAAGCGCCCGAAGCGACCACTGCCACGCCGAg                                                  <  8:85818/149‑1 (MQ=255)
                                                                                                 gTAGGCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAACTGGCTGGCCAGCGCGCAAGCGCCCGAAGCGACCACTGCCACGCCGAgca                                                <  4:141789/149‑1 (MQ=255)
                                                                                                 gTAGGCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAACTGGCTGGCCAGCGCGCAAGCGCCCGAAGCGACCACTGCCACGCCGAg                                                  >  8:186868/1‑147 (MQ=255)
                                                                                                 gTAGGCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGGCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAACTGGCTGGCCAGCGCGCAAGCGCCCGAAGCGACCACTGCCACGCCGAgca                                                <  2:133777/149‑1 (MQ=255)
                                                                                                    ggCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAACTGGCTGGCCAGCGCGCAAGCGCCCGAAGCGACCACTGCCACGCCGAgcagtg                                             <  7:186868/149‑1 (MQ=255)
                                                                                                     gCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCACTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAACTGGCTGGCCAGCGCGCAAGCGCCCGAAg                                                                     <  2:100443/124‑1 (MQ=255)
                                                                                                     gCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCACTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAACTGGCTGGCCAGCGCGCAAGCGCCCGAAg                                                                     >  1:100443/1‑124 (MQ=255)
                                                                                                       cTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTg                                                                                                    >  2:29667/1‑91 (MQ=255)
                                                                                                       cTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAACTGGCTGGCCAGCGCGCAAGCGCCCGAAGCGACCACTGCCACGCCGAgcagtgcag                                          >  5:173546/1‑149 (MQ=255)
                                                                                                        tGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAACTGGCTGGCCAGCGCGCAAGCGCCCGAAGCGACCACTGCCACGCCGAgcagtgcagt                                         >  3:228961/1‑149 (MQ=255)
                                                                                                        tGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAACTGGCTGGCCAGCGCGCAAGCGCCCGAAGCGACCACTGCCACGCCGAgcagtgcagt                                         >  5:275341/1‑149 (MQ=255)
                                                                                                         gCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAACTGGCTGGCCAGCGCGCAAGCGCCCGAAGCGACCACTGCCACGCCGAGcagtgcagtg                                        >  7:218952/1‑149 (MQ=255)
                                                                                                           cATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAACTGGCTGGCCAGCGCGCAAGCGCCCGAAGCGACCACTGCCACGCCGAGCAGTGCAGTGtt                                      <  4:228961/149‑1 (MQ=255)
                                                                                                                         cATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAACTGGCTGGCCAGCGCGCAAGCGCCCGAAGCGACCACTGCCACGCCGAGCAGTGTAGTGTTGGCTTTgctggctg                        >  3:222579/1‑149 (MQ=255)
                                                                                                                                ggggCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAACTGGCTGGCCAGCGCGCAAGCGCCCGAAGCGACCACTGCCACGCCGAGCAGTGTAGTGTTGGCTTTGCTGGCTGGGCCGAt                 <  4:222579/149‑1 (MQ=255)
                                                                                                                                ggggCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAACTGGCTGGCCAGCGCGCAAGCGCCCGAAGCGACCACTGCCACGCCGAGCAGTGCAGTGTTGGCTTTGCTGGCTGGGCCGAt                 <  6:275341/149‑1 (MQ=255)
                                                                                                                                ggggCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAACTGGCTGGCCAGCGCGCAAGCGCCCGAAGCGACCACTGCCACGCCGAGCAGTGCAGTGTTGGCTTTGCTGGCTGGGCCGAt                 >  4:257690/1‑149 (MQ=255)
                                                                                                                                   gcggcgCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAACTGGCTGGCCAGCGCGCAAGCGCCCGAAGCGACCACTGCCACGCCGAGCAGTGCAGTGTTGGCTTTGCTGGCTGGGCCGATCAg              >  4:292476/1‑149 (MQ=255)
                                                                                                                                     ggcgCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAACTGGCTGGCCAGCGCGCAAGCGCCCGAAGCGACCACTGCCACGCCGAGCAGTGCAGTGTTGGCTTTGCTGGCTGGGCCGATCAGCg            <  5:331412/149‑1 (MQ=255)
                                                                                                                                               aTTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAACTGGCTGGCCAGCGCGCAAGCGCCCGAAGCGACCACTGCCACGCCGAGCAGTGCAGTGTTGGCTTTGCTGGCTGGGCCGATCAGCGGGCGACGTGc  >  5:107675/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
CACCGGTGCGAATGCGCCAGGCTTCTTCCAGACGGTAGGCGTTCTGCGGGGTGATCTGGCGCAGCGGCGAGTAGCGGTCGCCATGCTCGGTACGGCCGTAGGCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAACTGGCTGGCCAGCGCGCAAGCGCCCGAAGCGACCACTGCCACGCCGAGCAGTGCAGTGTTGGCTTTGCTGGCTGGGCCGATCAGCGGGCGACGTGC  >  NC_002947/1646738‑1647029

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