Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,711,994 A→G 5.9% I1037T (ATC→ACC)  metH ← cobalamin‑dependent methionine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,711,9940AG 5.9% 198.0 / 2.1 67I1037T (ATC→ACC) metHcobalamin‑dependent methionine synthase
Reads supporting (aligned to +/- strand):  ref base A (34/29);  new base G (2/2);  total (36/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.53e-01

TCGTCGCCCTTGTCCTGGTAAGCCTTGGCCACTTCCTCGGCACCGATGCCGGCGGTGGTGATGAAACCACCTACGTAGTCGGTAACGCCACTGTCTTTCGGTGCGACGAAGTCGGCCAGCGACCAGTTGGGCTTGCCGTCCGGCTTGATGGTCTGCTGGCGCAGGTGATGCAGGGTGGCCAGTGCCTGGCCGTCTTCGCCGTAGACTTCGATGTCGTCATCGTCCACCTGGTTGGCCGGCCAGAAGCCGAATACCGCGCGAGCGCTGATCAGCTTCTCATCGATCAGCTTGTC  >  NC_002947/2711847‑2712139
                                                                                                                                                   |                                                                                                                                                 
tcgtcgCCCTTGTCCTGGTAAGCCTTGGCCACTTCCTCGGCACCGATGCCGGCGGTGGTGATGAAACCACCTACGTAGTCGGTAACGCCACTGTCTTTCGGTGCGACGAAGTCGGCCAGCGACCAGTTGGGCTTGCCGTCCGGCTTGAt                                                                                                                                                  <  6:295015/149‑1 (MQ=255)
  gtcgCCCTTGTCCTGGTAAGCCTTGGCCACTTCCTCGGCACCGATGCCGGCGGTGGTGATGAAACCACCTACGTAGTCGGTAACGCCACTGTCTTTCGGTGCGACGAAGTCGGCCAGCGACCAGTTGGGCTTGCCGTCCGGCTTGATg                                                                                                                                                 <  8:68596/148‑1 (MQ=255)
  gtcgCCCTTGTCATGGTAAGCCTTGGCCACTTCCTCGGCACCGATGCCGGCGGTGGTGATGAAACCACCTACGTAGTCGGTAACGCCACTGTCTTTCGGTGCGACGAAGTCGGCCAGCGACCAGTTGGGCTTGCCGTCCGGCTTGATgg                                                                                                                                                >  3:15415/1‑149 (MQ=255)
       ccTTGTCCTGGTAAGCCTTGGCCACTTCCTCGGCACCGATGCCGGCGGTGGTGATGAAACCACCTACGTAGTCGGTAACGCCACTGTCTTTCGGTGCGACGAAGTCGGCCAGCGACCAGTTGGGCTTGCCGTCCGGCTTGATGGTctgc                                                                                                                                           >  5:272963/1‑149 (MQ=255)
            tCCTGGTAAGCCTTGGCCACTTCCTCGGCACCGATGCCGGCGGTGGTGATGAAACCACCTACGTAGTCGGTAACGCCACTGTCTTTCGGTGCGACGAAGTCGGCCAGCGACCAGTTGGGCTTGCCGTCCGGCTTGATGGTCTGCTGgcg                                                                                                                                      >  3:229032/1‑149 (MQ=255)
              cTGGTAAGCCTTGGCCACTTCCTCGGCACCGATGCCGGCGGTGGTGATGAAACCACCTACGTAGTCGGTAACGCCACTGTCTTTCGGTGCGACGAAGTCGGCCAGCGACCAGTTGGGCTTGCCGTCCGGCTTGATGGTCTGCTGgcg                                                                                                                                      >  3:112425/1‑147 (MQ=255)
              cTGGTAAGCCTTGGCCACTTCCTCGGCACCGATGCCGGCGGTGGTGATGAAACCACCTACGTAGTCGGTAACGCCACTGTCTTTCGGTGCGACGAAGTCGGCCAGCGACCAGTTGGGCTTGCCGTCCGGCTTGATGGTCTGCTGGCGCa                                                                                                                                    <  6:272963/149‑1 (MQ=255)
                ggTAAGCCTTGGCCACTTCCTCGGCACCGATGCCGGCGGTGGTGATGAAACCACCTACGTAGTCGGTAACGCCACTGTCTTTCGGTGCGACGAAGTCGGCCAGCGACCAGTTGGGCTTGCCGTCCGGCTTGATGGTCTGCTGGCGCAgg                                                                                                                                  >  2:10750/1‑149 (MQ=255)
                 gTAAGCCTTGGCCACTTCCTCGGCACCGATGCCGGCGGTGGTGATGAAACCACCTACGTAGTCGGTAACGCCACTGTCTTTCGGTGCGACGAAGTCGGCCAGCGACCAGTTGGGCTTGCCGTCCGGCTTGATGGTCTGCTGGCGCAg                                                                                                                                   <  4:112425/147‑1 (MQ=255)
                   aaGCCTTGGCCACTTCCTCGGCACCGATGCCGGCGGTGGTGATGAAACCACCTACGTAGTCGGTAACGCCACTGTCTTTCGGTGCGACGAAGTCGGCCAGCGACCAGTTGGGCTTGCAGTCCGGCTTGATGGTCTGCTGGCGCAGGTGa                                                                                                                               >  8:40088/1‑149 (MQ=255)
                      ccTTGGCCACTTCCTCGGCACCGATGCCGGCGGTGGTGATGAAACCACCTACGTAGTCGGTAACGCCACTGTCTTTCGGTGCGACGAAGTCGGCCAGCGACCAGTTGGGCTTGCCGTCCGGCTTGATGGTCTGCTGGCGCAGGTGATg                                                                                                                             >  8:194805/1‑148 (MQ=255)
                       cTTGGCCACTTCCTCGGCACCGATGCCGGCGGTGGTGATGAAACCACCTACGTAGTCGGTAACGCCACTGTCTTTCGGTGCGACGAAGTCGGCCAGCGACCAGTTGGGCTTGCCGTCCGGCTTGATGGTCTGCTGGCGCAGGTGATg                                                                                                                             >  5:347616/1‑147 (MQ=255)
                       cTTGGCCACTTCCTCGGCACCGATGCCGGCGGTGGTGATGAAACCACCTACGTAGTCGGTAACGCCACTGTCTTTCGGTGCGACGAAGTCGGCCAGCGACCAGTTGGGCTTGCCGTCCGGCTTGATGGTCTGCTGGCGCAGGTGATg                                                                                                                             >  1:202351/1‑147 (MQ=255)
                             cACTTCCTCGGCACCGATGCCGGCGGTGGTGATGAAACCACCTACGTAGTCGGTAACGCCACTGTCTTTCGGTGCGACGAAGTCGGCCAGCGACCAGTTGGGCTTGCCGTCCGGCTTTATGGt                                                                                                                                               <  1:310202/123‑1 (MQ=255)
                             cACTTCCTCGGCACCGATGCCGGCGGTGGTGATGAAACCACCTACGTAGTCGGTAACGCCACTGTCTTTCGGTGCGACGAAGTCGGCCAGCGACCAGTTGGGCTTGCCGTCCGGCTTTATGGt                                                                                                                                               >  2:310202/1‑123 (MQ=255)
                                   cTCGGCACCGATGCCGGCGGTGGTGATGAAACCACCTACGTAGTCGGTAACGCCACTGTCTTTCGGTGCGACGAAGTCGGCCAGCGACCAGTTGGGCTTGCCGTCCGGCTTGATGGTCTGCTGGCGCAGGTGAt                                                                                                                              >  5:180882/1‑134 (MQ=255)
                                   cTCGGCACCGATGCCGGCGGTGGTGATGAAACCACCTACGTAGTCGGTAACGCCACTGTCTTTCGGTGCGACGAAGTCGGCCAGCGACCAGTTGGGCTTGCCGTCCGGCTTGATGGTCTGCTGGCGCAGGTGAt                                                                                                                              <  6:180882/134‑1 (MQ=255)
                                    tCGGCACCGATGCCGGCGGTGGTGATGAAACCACCTACGTAGTCGGTAACGCCACTGTCTTTCGGTGCGACGAAGTCGGCCAGCGACCAGTTGGGCTTGCCGTCCGGCTTGATGGTCTGCTGGCGCAGGTg                                                                                                                                >  8:127535/1‑131 (MQ=255)
                                    tCGGCACCGATGCCGGCGGTGGTGATGAAACCACCTACGTAGTCGGTAACGCCACTGTCTTTCGGTGCGACGAAGTCGGCCAGCGACCAGTTGGGCTTGCCGTCCGGCTTGATGGTCTGCTGGCGCAGGTg                                                                                                                                <  7:127535/131‑1 (MQ=255)
                                      ggCACCGATGCCGGCGGTGGTGATGAAACCACCTACGTAGTCGGTAACGCCACTGTCTTTCGGTGCGACGAAGTCGGCCAGCGACCAGTTGGGCTTGCCGTCCGGCTTGATGGTCTGCTGGCGCAGGTGATGCAGGGTGGCCAGTGc                                                                                                              <  7:40088/147‑1 (MQ=255)
                                          ccGATGCCGGCGGTGGTGATGAAACCACCTACGTAGTCGGTAACGCCACTGTCTTTCGGTGCGACGAAGTCTGCCAGCGACCAGTTGGGCTTGCCGTCCGGCTTGATGGTCTGCTGGCGCAGGTGATGCAGGGTGGTCAGTGCCTGGcc                                                                                                        <  4:229032/149‑1 (MQ=255)
                                           cGATGCCGGCGGTGGTGATGAAACCACCTACGTAGTCGGTAACGCCACTGTCTTTCGGTGCGACGAAGTCGGCCAGCGACCAGTTGGGCTTGCCGTCCGGCTTGATGGTCTGCTGGCGCAGGTGATGCAGGGTGGCCAGGCACTGGcca                                                                                                       <  2:202351/149‑2 (MQ=255)
                                                        ggtgATGAAACCACCTACGTAGTCGGTAACGCCACTGTCTTTCGGTGTGACGAAGTCGGCCAGCGACCAGTTGGGCTTGCCGTCCGGCTTGGTGGTCTGCTGGCGCAGGTGATGCAGGGTGGCCAGTGCCTGGcc                                                                                                        >  1:82694/1‑135 (MQ=255)
                                                        ggtgATGAAACCACCTACGTAGTCGGTAACGCCACTGTCTTTCGGTGTGACGAAGTCGGCCAGCGACCAGTTGGGCTTGCCGTCCGGCTTGGTGGTCTGCTGGCGCAGGTGATGCAGGGTGGCCAGTGCCTGGcc                                                                                                        >  1:77878/1‑135 (MQ=255)
                                                        ggtgATGAAACCACCTACGTAGTCGGTAACGCCACTGTCTTTCGGTGCGACGAAGTCGGCCAGCGACCAGTTGGGCTTGCCGTCCGGCTTGGTGGTCTGCTGGCGCAGGTGATGCAGGGTGGCCAGTGCCTGGcc                                                                                                        <  2:82694/135‑1 (MQ=255)
                                                        ggtgATGAAACCACCTACGTAGTCGGTAACGCCACTGTCTTTCGGTGCGACGAAGTCGGCCAGCGACCAGTTGGGCTTGCCGTCCGGCTTGGTGGTCTGCTGGCGCAGGTGATGCAGGGTGGCCAGTGCCTGGcc                                                                                                        <  2:77878/135‑1 (MQ=255)
                                                               aaaCCACCTACGTAGTCGGTAACGCCACTGTCTTTCGGTGCGACGAAGTCGGCCAGCGACCAGTTGGGCTTGCCGTCCGGCTTGATGGTCTGCTGGCGCAGGTGATGCAGGGTGGCCAGTGCCTGGCCGTCTTCGCCGTAGACTTCGa                                                                                    <  4:15415/148‑1 (MQ=255)
                                                                aaCCACCTACGTAGTCGGTAACGCCACTGTCTTTCGGTGCGACGAAGTCGGCCAGCGACCAGTTGGGCTTGCCGTCCGGCTTGATGGTCTGCTGGCGCAGGTGATGCAGGGTGGCCAGTGCCTGGCCGTCTTCGCCGTAGACtt                                                                                       >  2:286777/1‑144 (MQ=255)
                                                                aaCCACCTACGTAGTCGGTAACGCCACTGTCTTTCGGTGCGACGAAGTCGGCCAGCGACCAGTTGGGCTTGCCGTCCGGCTTGATGGTCTGCTGGCGCAGGTGATGCAGGGTGGCCAGTGCCTGGCCGTCTTCGCCGTAGACtt                                                                                       <  1:286777/144‑1 (MQ=255)
                                                                      cTACGTAGTCGGTAACGCCACTGTCTTTCGGTGCGACGAAGTCGGCCAGCGACCAGTTGGGCTTGCCGTCCGGCTTGATGGTCTGCTGGCGCAGGTGATGCAGGGTGGCCAGTGCCTGGCCGTCTTCGCCGTAGACTTCGATgtc                                                                                <  6:192528/145‑1 (MQ=255)
                                                                      cTACGTAGTCGGTAACGCCACTGTCTTTCGGTGCGACGAAGTCGGCCAGCGACCAGTTGGGCTTGCCGTCCGGCTTGATGGTCTGCTGGCGCAGGTGATGCAGGGTGGCCAGTGCCTGGCCGTCTTCGCCGTAGACTTCGATgtc                                                                                >  5:192528/1‑145 (MQ=255)
                                                                      cTACGTAGTCGGTAACGCCACTGTCTTTCGGTGCGACGAAGTCGGCCAGCGACCAGTTGGGCTTGCCGTCCGGCTAGATGGTCTGCTGGCGCAGGTGATGCAGGGTGGCCAGTGCCTGGcc                                                                                                        >  1:140582/1‑121 (MQ=255)
                                                                      cTACGTAGTCGGTAACGCCACTGTCTTTCGGTGCGACGAAGTCGGCCAGCGACCAGTTGGGCTTGCCGTCCGGCTAGATGGTCTGCTGGCGCAGGTGATGCAGGGTGGCCAGTGCCTGGcc                                                                                                        <  2:140582/121‑1 (MQ=255)
                                                                           tAGTCGGTAACGCCACTGTCTTTCGGTGCGACGAAGTCGGCCAGCGACCAGTTGGGCTTGCCGTCCGGCTTGATGGTCTGCTGGCGCAGGTGATGCAGGGTGGCCAGTGCCTGGCCGTCTTCGCCGTAGACTTCGATGTCGTCATCGTc                                                                       >  5:350630/1‑149 (MQ=255)
                                                                             gTCGGTAACGCCACTGTCTTTCGGTGCGACGAAGTCGGCCAGCGACCAGTTGGGCTTGCCGTCCGGCTTGATGGTCTGCTGGCGCAg                                                                                                                                   <  3:94443/87‑1 (MQ=255)
                                                                              tCGGTAACGCCACTGTCTTTCGGTGCGACGAAGTCGGCCAGCGACCAGTTGGGCTTGCCGTCCGGCTTGATGGTCTGCTGGCGCAGGTGATGCAGGGTGGCCAGTGCCTGGCCGTCTTCGCCGTAGACTTCGATGTCGTCATCGTCCAc                                                                    <  1:10750/149‑1 (MQ=255)
                                                                              tCGGTAACGCCACTGTCTTTCGGTGCGACGAAGTCGGCCAGCGACCAGTTGGGCTTGCCGTCCGGCTTGATGGTCTGCTGGCGCAGGTGATGCAGGGTGGCCAGTGCCTGGCCGTCTTCGCCGTAGACTTCGATGTCGTCATCGTCCAc                                                                    >  5:196694/1‑149 (MQ=255)
                                                                              tCGGTAACGCCACTGTCTTTCGGTGCGACGAAGTCGGCCAGCGACCAGTTGGGCTTGCCGTCCGGCTTGATGGTCTGCTGGCGCAGGTGATGCAGGGTGGCCAGTGCCTGGCCGTCTTCGCCGTAGACTTCGATGTCGTCATCGTCCAc                                                                    >  6:163233/1‑149 (MQ=255)
                                                                                    aCGCCACTGTCTTTCGGTGCGACGAAGTCGGCCAGCGACCAGTTGGGCTTGCCGTCCGGCTTGATGGTCTGCTGGCGCAGGTGATGCAGGGTGGCCAGTGCCTGGCCGTCTTCGCCGTAGACTTCGATGTCGTCATCGTCCACCTGGt                                                               <  7:194805/148‑1 (MQ=255)
                                                                                      gCCACTGTCTTTCGGTGCGACGAAGTCGGCCAGCGACCAGTTGGGCTTGCCGTCCGGCTTGATGGTCTGCTGGCGCAGGTGATGCAGGGTGGCCAGTGCCTGGCCGTCTTCGCCGTAGACTTCGATGTCGTCATCGTCCACCTGGTTgg                                                            >  4:318211/1‑149 (MQ=255)
                                                                                        cACTGTCTTTCGGTGCGACGAAGTCGGCCAGCGACCAGTTGGGCTTGCCGTCCGGCTTGATGGTCTGCTGGCGCAGGTGATGCAGGGTGGCCAGTGCCTGGCCGTCTTCGCCGTAGACTTCGATGTCGTCATCGTCCACCTGGTTggcc                                                          <  6:196694/149‑1 (MQ=255)
                                                                                            gTCTTTCGGTGCGACGAAGTCGGCCAGCGACCAGTTGGGCTTGCCGTCCGGCTTGATGGTCTGCTGGCGCAGGTGATGCAGGGTGGCCAGTGCCTGGCCGTCTTCGCCGTAGACTTCGATGTCGTCATCGTCCACCTGGTTggccggc                                                       >  7:134273/1‑148 (MQ=255)
                                                                                                ttCGGTGCGACGAAGTCGGCCAGCGACCAGTTGGGCTTGCCGTCCGGCTTGATGGTCTGCTGGCGCAGGTGATGCAGGGTGGCCAGTGCCTGGCCGTCTTCGCCGTAGACTTCGATTTCGTCATCGTCCACCTGGTTGGCCGGCCAg                                                    >  2:187687/1‑147 (MQ=255)
                                                                                                   ggTGCGACGAAGTCGGCCAGCGACCAGTTGGGCTTGCCGTCCGGCTTGATGGTCTGCTGGCGCAGGTGATGCAGGGTGGCCAGTGCCTGGCCGTCTTCGCCGTAGACTTCGATGTCGTCATCGTCCACCTGGTTGGCCGGCCAGAAGcc                                               <  6:350630/149‑1 (MQ=255)
                                                                                                       cgacgaAGTCGGCCAGCGACCAGTTGGGCTTGCCGTCCGGCTTGATGGTCTGCTGGCGCAGGTGATGCAGGGTGGCCAGTGCCTGGCCGTCTTCGCCGTAGACTTCGATGTCGTCATCGTc                                                                       <  4:34102/121‑1 (MQ=255)
                                                                                                       cgacgaAGTCGGCCAGCGACCAGTTGGGCTTGCCGTCCGGCTTGATGGTCTGCTGGCGCAGGTGATGCAGGGTGGCCAGTGCCTGGCCGTCTTCGCCGTAGACTTCGATGTCGTCATCGTc                                                                       >  3:34102/1‑121 (MQ=255)
                                                                                                       cgacgaAGTCGGCCAGCGACCAGTTGGGCTTGCCGTCCGGCTTGATGGTCTGCTGGCGCAGGTGATGCAGGGTGGCCAGTGCCTGGCCGTCTTCGCCGTAGACTTCGATGTCGTCATCGTCCACCTGGTTggccggc                                                       >  5:333202/1‑137 (MQ=255)
                                                                                                       cgacgaAGTCGGCCAGCGACCAGTTGGGCTTGCCGTCCGGCTTGATGGTCTGCTGGCGCAGGTGATGCAGGGTGGCCAGTGCCTGGCCGTCTTCGCCGTAGACTTCGATGTCGTCATCGTCCACCTGGTTggccggc                                                       <  6:333202/137‑1 (MQ=255)
                                                                                                        gacgaAGTCGGCCAGCGACCAGTTGGGCTTGCCGTCCGGCTTGATGGTCTGCTGGCGCAGGTGATGCAGGGTGGCCAGTGCCTGGCCGTCTTCGCCGTAGACTTCGATGTCGTCATCGTCCACCTGGTTGGCCGGCCAGAAGCCGAATa                                          >  3:288894/1‑149 (MQ=255)
                                                                                                        gacgaAGTCGGCCAGCGACCAGTTGGGCTTGCCGTCCGGCTTGATGGTCTGCTGGCGCAGGTGATGCAGGGTGGCCAGTGCCTGGCCGTCTTCGCCGTAGACTTCGATGTCGTCATCGTCCACCTGGTTGGCCGGCCAGAAGCCGAATa                                          >  6:321362/1‑149 (MQ=255)
                                                                                                           gaAGTCGGCCAGCGACCAGTTGGGCTTGCCGTCCGGCTTGATGGTCTGCTGGCGCAGGTGATGCAGGGTGGCCAGTGCCTGGCCGTCTTCGCCGTAGACTTCGATGTCGTCATCGTCCACCTGGTTGGCCGGCCAGAAGCCGAATACcg                                       >  6:184744/1‑149 (MQ=255)
                                                                                                                 ggCCAGCGACCAGTTGGGCTTGCCGTCCGGCTTGATGGTctg                                                                                                                                            >  8:269652/1‑42 (MQ=255)
                                                                                                                  gCCAGCGACCAGTTGGGCTTGCCGTCCGGCTTGATGGTCTGCTGGCGCAGGTGATGCAGGGTGGCCAGTGCCTGGCCGTCTTCGCCGTAGACTTCGATGTCGTCATCGTCCACCTGGTTGGCCGGCCAGAAGCCGAATACCGCGCGAgc                                >  2:232184/1‑149 (MQ=255)
                                                                                                                  gCCAGCGACCAGTTGGGCTTGCCGTCCGGCTTGATGGTCTGCTGGCGCAGGTGATGCAGGGTGGCCAGTGCCTGGCCGTCTTCGCCGTAGACTTCGATGTCGTCATCGTCCACCTGGTTGGCCGGCCAGAAGCCGAATACCGCGCGAgc                                >  2:27232/1‑149 (MQ=255)
                                                                                                                     aGCGACCAGTTGGGCTTGCCGTCCGGCTTGATGGTCTGCTGGCGCAGGTGATGCAGGGTGGCCAGTGCCTGGCCGTCTTCGCCGTAGACTTCGATGTCGTCATCGTCCACCTGGTTGGCCGGCCAGAAGCCGAATACCGCGCGAGCGCt                             <  6:347616/149‑1 (MQ=255)
                                                                                                                       cGACCAGTTGGGCTTGCCGTCCGGCTTGATGGTCTGCTGGCGCAGGTGATGCAGGGTGGCCAGTGCCTGGCCGTCTTCGCCGTAGACTTCGATGTCGTCATCGTCCACCTGGTTGGCCGGCCAGAAGCCGAATACCGCGCGAGCGCTGa                           >  7:316143/1‑149 (MQ=255)
                                                                                                                        gACCAGTTGGGCTTGCCGTCCGGCTTGATGGTCTGCTGGCGCAGGTGATGCAGGGTGGCCAGTGCCTGGCCGTCTTCGCCGTAGACTTCGATGTCGTCATCGTCCACCTGGTTGGCCGGCCAGAAGCCGAATACCGCGCGAGCGCTGAt                          <  8:316143/149‑1 (MQ=255)
                                                                                                                         aCCAGTTGGGCTTGCCGTCCGGCTTGATGGTCTGCTGGCGCAGGTGATGCAGGGTGGCCAGTGCCTGGCCGTCTTCGCCGTAGACtt                                                                                       <  2:207206/87‑1 (MQ=255)
                                                                                                                         aCCAGTTGGGCTTGCCGTCCGGCTTGATGGTCTGCTGGCGCAGGTGATGCAGGGTGGCCAGTGCCTGGCCGTCTTCGCCGTAGACtt                                                                                       >  1:207206/1‑87 (MQ=255)
                                                                                                                            aGTTGGGCTTGCCGTCCGGCTTGATGGTCTGCTGGCGCAGGTGATGCAGGGTGGCCAGTGCCTGGCCGTCTTCGCCGTAGACTTCGATGTCGTCATCGTCCACCTGGTTGGCCGGCCAGAAGCCGAATACCGCGCGAGCGCTGATCAGc                      <  1:232184/149‑1 (MQ=255)
                                                                                                                                gggCTTGCCGTCCGGCTTGATGGTCTGCTGGCGCAGGTGATGCAGGGTGGCCAGTGCCTGGCCGTCTTCGCCGTAGACTTCGATGTCGTCATCGTCCACCTGGTTGGCCGGCCAGAAGCCGAATACCGCGCGAGCGCTGATCAGCTtct                  >  8:47874/1‑149 (MQ=255)
                                                                                                                                  gCTTGCCGTCCGGCTTGATGGTCTGCTGGCGCAGGTGATGCAGGGTGGCCAGTGCCTGGCCGTCTTCGCCGTAGACTTCGATGTCGTCATCGTCCACCTGGTTGGCCGGCCAGAAGCCGAATACCGCGCGAGCGCTGATCAGCTTCTCa                <  3:318211/149‑1 (MQ=255)
                                                                                                                                   cTTGCCGTCCGGCTTGATGGTCTGCTGGCGCAGGTGATGCATGGTGGCCAGTGCCTGGCCGTCTTCGCCGTAGACTTCGATGTCGTCATCGTCCACCTGGTTGGCCGGCCAGAAGCTGATCAGCGCGCGCGCGGTATTCGGCTtctggc               <  7:243564/149‑4 (MQ=38)
                                                                                                                                   cTTGCCGTCCGGCTTGATGGTCTGCTGGCGCAGGTGATGCAGGGTGGCCAGTGCCTGGCCGTCTTCGCCGTAGACTTCGATGTCGTCATCGTCCACCTGGTTGGCCGGCCAGAAGCCGAATACCGCGCGAGCGCTGATCAGCTTCTCAt               >  3:239438/1‑149 (MQ=255)
                                                                                                                                   cTTGCCGTCCGGCTTGATGGTCTGCTGGCGCAGGTGATGCAGGGTGGCCAGTGCCTGGCCGTCTTCGCCGTAGACTTCGATGTCGTCATCGTCCACCTGGTTGGCCGGCCAGAAGCCGAATACCGCGCGAGCGCTGATCAGCTTCTCAt               <  5:321362/149‑1 (MQ=255)
                                                                                                                                       ccGTCCGGCTTGATGGTCTGCTGGCGCAGGTGATGCAGGGTGGCCAGTGCCTGGCCGTCTTCGCCGTAGACTTCGATGTCGTCATCGTCCACCTGGTTGGCCGGCCAGAAGCCGAATACCGCGCGAGCGCTGATCAGCTTCTCATCGa            <  8:134273/148‑1 (MQ=255)
                                                                                                                                                  gATGGTCTGCTGGCGCAGGTGATGCAGGGTGGCCAGTGCCTGGCCGTCTTCGCCGTAGACTTCGATGTCGTCATCGTCCACCTGGTTGGCCGGCCAGAAGCCGAATACCGCGCGAGCGCTGATCAGCTTCTCATCGATCAGCTTGTc  >  2:180536/1‑147 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TCGTCGCCCTTGTCCTGGTAAGCCTTGGCCACTTCCTCGGCACCGATGCCGGCGGTGGTGATGAAACCACCTACGTAGTCGGTAACGCCACTGTCTTTCGGTGCGACGAAGTCGGCCAGCGACCAGTTGGGCTTGCCGTCCGGCTTGATGGTCTGCTGGCGCAGGTGATGCAGGGTGGCCAGTGCCTGGCCGTCTTCGCCGTAGACTTCGATGTCGTCATCGTCCACCTGGTTGGCCGGCCAGAAGCCGAATACCGCGCGAGCGCTGATCAGCTTCTCATCGATCAGCTTGTC  >  NC_002947/2711847‑2712139

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