Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,099,295 G→T 6.1% R91S (CGT→AGT)  kdsD ← D‑arabinose 5‑phosphate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,099,2950GT 6.1% 208.9 / 3.8 66R91S (CGT→AGT) kdsDD‑arabinose 5‑phosphate isomerase
Reads supporting (aligned to +/- strand):  ref base G (33/29);  new base T (2/2);  total (35/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 = 1.19e-01

GACCTCGGCGGCCTGGGCCAGCGGCGAGTCCGGGTTGCCGGTCAGGCTGATCAGCTTGATGCCAAGGCGCTTGATGAACGGCAGCAGGGTGACGATTTCGGCGGTGCTGCCGGAGTTCGACAGCGCAAGGATGACATCATCACGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGGCGATCTTGTTGCCGATGTGCCCGGACTTGCCCATGCCGACCACGACGACCCGGCCCTTGCTG  >  NC_002947/1099152‑1099440
                                                                                                                                               |                                                                                                                                                 
gACCTCGGCGGCCTGGGCCAGCGGCGAGTCCGGGTTGCCGGTCAGGCTGATCAGCTTGATGCCAAGGCGCTTGATGAACGGCAGCAGGGTGACGATTTCGGCGGTGCTGCCGGAGTTCGACAGCGCAAGGATGACATCATCACGGGTGa                                                                                                                                              >  4:217071/1‑149 (MQ=255)
  ccTCGGCGGCCTGGGCCAGCGGCGAGTCCGGGTTGCCGGTCAGGCTGATCAGCTTGATGCCAAGGCGCTTGATGAACGGCAGCAGGGTGACGATTTCGGCGGTGCTGCCGGAGTTCGACAGCGCAAGGATGACATCATCACGGGTGAt                                                                                                                                             >  8:316755/1‑148 (MQ=255)
  ccTCGGCGGCCTGGGCCAGCGGCGAGTCCGGGTTGCCGGTCAGGCTGATCAGCTTGATGCCAAGGCGCTTGATGAACGGCAGCAGGGTGACGATTTCGGCGGTGCTGCCGGAGTTCGACAGCGCAAGGATGACATCATCACGGGTGAt                                                                                                                                             >  3:16042/1‑148 (MQ=255)
          gcCTGGGCCAGCGGCGAGTCCGGGTTGCCGGTCAGGCTGATCAGCTTGATGCCAAGGCGCTTGATGAACGGCAGCAGGGTGACGATTTCGGCGGTGCTGCCGGAGTTCGACAGCGCAAGGATGACATCATCACGGGTGATCATGccc                                                                                                                                      <  8:354320/147‑1 (MQ=255)
          gcCTGGGCCAGCGGCGAGTCCGGGTTGCCGGTCAGGCTGATCAGCTTGATGCCAAGGCGCTTGATGAACGGCAGCAGGGTGACGATTTCGGCGGTGCTGCCGGAGTTCGACAGCGCAAGGATGACATCATCACGGGTGATCATGccc                                                                                                                                      >  7:354320/1‑147 (MQ=255)
          gcCTGGGCCAGCGGCGAGTCCGGGTTGCCGGTCAGGCTGATCAGCTTGATGCCAAGGCGCTTGATGAACGGCAGCAGGGTGACGATTTCGGCGGTGCTGCCGGAGTTCGACAGCGCAAGGATGACATCATCACGGGTGATCATGccc                                                                                                                                      <  6:70955/147‑1 (MQ=255)
          gcCTGGGCCAGCGGCGAGTCCGGGTTGCCGGTCAGGCTGATCAGCTTGATGCCAAGGCGCTTGATGAACGGCAGCAGGGTGACGATTTCGGCGGTGCTGCCGGAGTTCGACAGCGCAAGGATGACATCATCACGGGTGATCATGccc                                                                                                                                      >  5:70955/1‑147 (MQ=255)
                            tCCGGGTTGCCGGTCAGGCTGATCAGCTTGATGCCAAGGCGCTTGATGAACGGCAGCAGGGTGACGATTTCGGCGGTGCTGCCGGAGTTCGACAGCGCAAGGATGACATCATCACGGGTGATCATGCCCATGTCGCCATGGCTGGATTc                                                                                                                  >  6:266578/1‑149 (MQ=255)
                              cGGGTTGCCGGTCAGGCTGATCAGCTTGATGCCAAGGCGCTTGATGAACGGCAGCAGGGTGACGATTTCGGCGGTGCTGCCGGAGTTCGACAGCGCAAGGATGACATCATCACTGGTGATCATGCCCATGTCGCCATGGCTGGCTTCgg                                                                                                                >  4:52415/1‑149 (MQ=255)
                              cGGGTTGCCGGTCAGGCTGATCAGCTTGATGCCAAGGCGCTTGATGAACGGCAGCAGGGTGACGATTTCGGCGGTGCTGCCGGAGTTCGACAGCGCAAGGATGACATCATCACTGGTGATCATGCCCATGTCGCCATGGCTGGCTTCgg                                                                                                                >  4:57367/1‑149 (MQ=255)
                                   tGCCGGTCAGGCTGATCAGCTTGATGCCAAGGCGCTTGATGAACGGCAGCAGGGTGACGATTTCGGCGGTGCTGCCGGAGTTCGACAGCGCAAGGATGACATCATCACGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGa                                                                                                           >  5:133561/1‑149 (MQ=255)
                                      cGGTCAGGCTGATCAGCTTGATGCCAAGGCGCTTGATGAACGGCAGCAGGGTGACGATTTCGGCGGTGCTGCCGGAGTTCGACAGCGCAAGGATGACATCATCACGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGc                                                                                                        <  5:266578/149‑1 (MQ=255)
                                       ggTCAGGCTGATCAGCTTGATGCCAAGGCGCTTGATGAACGGCAGCAGGGTGACGATTTCGGCGGTGCTGCCGGAGTTCGACAGCGCAAGGATGACATCATCACggg                                                                                                                                                 <  6:286967/107‑1 (MQ=255)
                                       ggTCAGGCTGATCAGCTTGATGCCAAGGCGCTTGATGAACGGCAGCAGGGTGACGATTTCGGCGGTGCTGCCGGAGTTCGACAGCGCAAGGATGACATCATCACggg                                                                                                                                                 >  5:286967/1‑107 (MQ=255)
                                        gTCAGGCTGATCAGCTTGATGCCAAGGCGCTTGATGAACGGCAGCAGGGTGACGATTTCGGCGGTGCTGCCGGAGTTCGACAGCGCAAGGATGACATCATCACGGGTGATCATGCCCATGTCGCCATCGCTGGCATCGCCCGGATGCa                                                                                                       >  6:349818/1‑148 (MQ=255)
                                             gCTGATCAGCTTGATGCCAAGGCGCTTGATGAACGGCAGCAGGGTGACGATTTCGGCGGTGCTGCCGGAGTTCGACAGCGCAAGGATGACATCATCACGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGaaaa                                                                                                 <  3:80598/149‑1 (MQ=255)
                                             gCTGATCAGCTTGATGCCAAGGCGCTTGATGAACGGCAGCAGGGTGACGATTTCGGCGGTGCTGCCGGAGTTCGACAGCGCAAGGATGACATCATCACGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGaaa                                                                                                  >  4:80598/1‑148 (MQ=255)
                                              cTGATCAGCTTGATGCCAAGGCGCTTGATGAACGGCAGCAGGGTGACGATTTCGGCGGTGCTGCCGGAGTTCGACAGCGCAAGGATGACATCATCACGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGaaaaa                                                                                                >  5:128425/1‑149 (MQ=255)
                                              cTGATCAGCTTGATGCCAAGGCGCTTGATGAACGGCAGCAGGGTGACGATTTCGGCGGTGCTGCCGGAGTTCGACAGCGCAAGGATGACATCATCACGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGaaaaa                                                                                                <  6:133561/149‑1 (MQ=255)
                                               tGATCAGCTTGATGCCAAGGCGCTTGATGAACGGCAGCAGGGTGACGATTTCGGCGGTGCTGCCGGAGTTCGACAGCGCAAGGATGACATCATCACGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAAc                                                                                               >  5:349310/1‑149 (MQ=255)
                                                  tCAGCTTGATGCCAAGGCGCTTGATGAACGGCAGCAGGGTGACGATTTCGGCGGTGCTGCCGGAGTTCGACAGCGCAAGGATGACATCATCACGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGaa                                                                                                   <  7:82611/142‑1 (MQ=255)
                                                  tCAGCTTGATGCCAAGGCGCTTGATGAACGGCAGCAGGGTGACGATTTCGGCGGTGCTGCCGGAGTTCGACAGCGCAAGGATGACATCATCACGGGTGATCATGCCCATGTCGCCATGCCTGGCTTCGGCCGGATGCACGaa                                                                                                   >  8:82611/1‑142 (MQ=255)
                                                  tCAGCGTGATGCCAAGGCGCTTGATGAACGGCAGCAGGGTGACGATTTCGGCGGTGCTGCCGGAGTTCGACAGCGCAAGGATGACATCATCACGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGc                                                                                             <  6:339698/148‑1 (MQ=255)
                                                         gATGCCAAGGCGCTTGATGAACGGCAGCAGGGTGACGATTTCGGCGGTGCTGCCGGAGTTCGACAGCGCAAGGATGACATCATCACGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGgtgt                                                                                       <  5:187714/147‑1 (MQ=255)
                                                             ccAAGGCGCTTGATGAACGGCAGCAGGGTGACGATTTCGGCGGTGCTGCCGGAGTTCGACAGCGCAAGGATGACATCATCACGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGt                                                                                 >  3:211945/1‑149 (MQ=255)
                                                               aaGGCGCTTGATGAACGGCAGCAGGGTGACGATTTCGGCGGTGCTGCCGGAGTTCGACAGCGCAAGGATGACATCATCACGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGc                                                                               >  1:218353/1‑149 (MQ=255)
                                                                    gcTTGATGAACGGCAGCAGGGTGACGATTTCGGCGGTGCTGCCGGAGTTCGACAGCGCAAGGATGACATCATCACGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGc                                                                           >  6:354811/1‑148 (MQ=255)
                                                                     cTTGATGAACGGCAGCAGGGTGACGATTTCGGCGGTGCTGCCGGAGTTCGACAGCGCAAGGATGACATCATCACGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCa                                                                         <  6:349310/149‑1 (MQ=255)
                                                                        gatgaACGGCAGCAGGGTGACGATTTCGGCGGTGCTGCCGGAGTTCGACAGCGCAAGGATGACATCATCACGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAgg                                                                       <  5:349818/148‑1 (MQ=255)
                                                                         atgaACGGCAGCAGGGTGACGATTTCGGCGGTGCTGCCGGAGTTCGACAGCGCAAGGATGACATCATCACGGGTGAACATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGt                                                                     >  7:195860/1‑149 (MQ=255)
                                                                                 cagcagGGTGACGATTTCGGCGGTGCTGCCGGAGTTCGACAGCGCAAGGATGACATCATCACTGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGGCGa                                                             <  3:57367/149‑1 (MQ=255)
                                                                                 cagcagGGTGACGATTTCGGCGGTGCTGCCGGAGTTCGACAGCGCAAGGATGACATCATCACTGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGGCGa                                                             <  3:52415/149‑1 (MQ=255)
                                                                                  agcagGGTGACGATTTCGGCGGTGCTGCCGGAGTTCGACAGCGCAAGGATGACATCATCACGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGGCGAt                                                            <  6:128425/149‑1 (MQ=255)
                                                                                   gcagGGTGACGATTTCGGCGGTGCTGCCGGAGTTCGACAGCGCAAGGATGACATCATCACGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGGCGATc                                                           <  7:316755/149‑1 (MQ=255)
                                                                                    cagGGTGACGATTTCGGCGGTGCTGCCGGAGTTCGACAGCGCAAGGATGACATCATCACGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGGCGATCt                                                          >  3:271774/1‑149 (MQ=255)
                                                                                       ggTGACGATTTCGGCGGTGCTGCCGGAGTTCGACAGCGCAAGGATGACATCATCACGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGGCGATCttgt                                                       <  1:227743/149‑1 (MQ=255)
                                                                                         tGACGATTTCGGCGGTGCTGCCGGAGTTCGACAGCGCAAGGATGACATCATCACGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGGCGATCttgttg                                                     >  1:360490/1‑149 (MQ=255)
                                                                                         tGACGATTTCGGCGGTGCTGCCGGAGTTCGACAGCGCAAGGATGACATCATCACGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGGCGATCttgttg                                                     <  3:217071/149‑1 (MQ=255)
                                                                                           aCGATTTCGGCGGTGCTGCCGGAGTTCGACAGCGCAAGGATGACATCATCACGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGGCGATCTTGTTGcc                                                   >  3:138867/1‑149 (MQ=255)
                                                                                            cGATTTCGGCGGTGCTGCCGGAGTTCGACAGCGCAAGGATGACATCATCACGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGGCGATCTTGTTGcc                                                   >  4:349234/1‑148 (MQ=255)
                                                                                             gATTTCGGCGGTGCTGCCGGAGTTCGACAGCGCAAGGATGACATCATCACGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGGCGATCTTGTTGCCGa                                                 >  1:335793/1‑149 (MQ=255)
                                                                                             gATTTCGGCGGTGCTGCCGGAGTTCGACAGCGCAAGGATGACATCATCACGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGGCGATCTTGTTGCCGa                                                 >  6:167074/1‑149 (MQ=255)
                                                                                             gATTTCGGCGGTGCTGCCGGAGTTCGACAGCGCAAGGATGACATCATCACGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGGCGATCTTGTTGCCGa                                                 >  1:149129/1‑149 (MQ=255)
                                                                                                 tCGGCGGTGCTGCCGGAGTTCGACAGCGCAAGGATGACATCATCACGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGGCGATCTTGTTGCCGa                                                 <  6:286217/145‑1 (MQ=255)
                                                                                                 tCGGCGGTGCTGCCGGAGTTCGACAGCGCAAGGATGACATCATCACGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGGCGATCTTGTTGCCGa                                                 >  5:286217/1‑145 (MQ=255)
                                                                                                 tCGGCGGTGCTGCCGGAGTTCGACAGCGCAAGGATGACATCATCACGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGGCGATCTTGTTGCCGAtgt                                              <  5:354811/148‑1 (MQ=255)
                                                                                                   ggcggTGCTGCCGGAGTTCGACAGCGCAAGGATGACATCATCACGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGGCGATCTTGTTGCCGATGTGcc                                           <  8:144840/149‑1 (MQ=255)
                                                                                                   ggcggTGCTGCCGGAGTTCGACAGCGCAAGGATGACATCATCACGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGGCGATCTTGTTGCCGATGTGcc                                           >  4:229326/1‑149 (MQ=255)
                                                                                                       gTGCTGCCGGAGTTCGACAGCGCAAGGATGACATCATCACGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGGCGATCTTGTTGCCGATGTGCCCGGa                                       <  4:211945/149‑1 (MQ=255)
                                                                                                             ccGGAGTTCGACAGCGCAAGGATGACATCATCACGGGTGATCATGCCCATGTCGCCATGGTTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGGCGATCttgtt                                                      <  1:301294/128‑1 (MQ=255)
                                                                                                             ccGGAGTTCGACAGCGCAAGGATGACATCATCACGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGGCGATCttgtt                                                      >  2:301294/1‑128 (MQ=255)
                                                                                                             ccGGAGTTCGACAGCGCAAGGATGACATCATCACGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGGCGATCTTGTTGCCGATTTGCCCGGACTTGcc                                 <  8:195860/149‑1 (MQ=255)
                                                                                                             ccGGAGTTCGACAGCGCAAGGATGACATCATCACGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGGCGATCTTGTTGCCGATGTGCCCGGACTTGc                                  <  2:360490/148‑1 (MQ=255)
                                                                                                                   ttCGACAGCGCAAGGATGACATCATCACGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGGCGATCTTGTTGCCGATGTGCCCGGACTTGCCCATGcc                           >  4:96404/1‑149 (MQ=255)
                                                                                                                      gACAGCGCAAGGATGACATCATCACGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGGCGATCTTGTTGCCGATGTGCCCGGACTTGCCCATg                             <  7:14032/144‑1 (MQ=255)
                                                                                                                      gACAGCGCAAGGATGACATCATCACGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGGCGATCTTGTTGCCGATGTGCCCGGACTTGCCCATg                             >  8:14032/1‑144 (MQ=255)
                                                                                                                      gACAGCGCAAGGATGACATCATCACGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGGCGATCTTGTTGCCGATGTGCCCGGACTTGCCCATGCCGAc                        >  4:294444/1‑149 (MQ=255)
                                                                                                                       aCAGCGCAAGGATGACATCATCACGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGGCGATCTTGTTGCCGATGTGCCCGGACTTGCCCATGCCGAcc                       <  2:335793/149‑1 (MQ=255)
                                                                                                                           cgcAAGGATGACATCATCACGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGGCGATCTTGTTGCCGATGTGCCCGGACTTGCCCATGCCGACCacga                   <  3:349234/149‑1 (MQ=255)
                                                                                                                                  aTGACATCATCACGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGGCGATCTTGTTGCCGATGTGCCCGGACTTGCCCATGCCGACCACGACGAccc             <  3:229326/148‑1 (MQ=255)
                                                                                                                                   tGACATCATCACGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGGCGATCTTGTTGCCGATGTGCCCGGACTTGCCCATGCCGACCACGACGACCCgg           <  4:16042/149‑1 (MQ=255)
                                                                                                                                   tGACATCATCACGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGGCGATCTTGTTGCCGATGTGCCCGGACTTGCCCATGCCGACCACGACGACCCgg           <  2:149129/149‑1 (MQ=255)
                                                                                                                                     aCATCATCACGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGGCGATCTTGTTGCCGATGTGCCCg                                         >  2:177721/1‑117 (MQ=255)
                                                                                                                                     aCATCATCACGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGGCGATCTTGTTGCCGATGTGCCCg                                         <  1:177721/117‑1 (MQ=255)
                                                                                                                                      catcatCACGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGGCGATCTTGTTGCCGATGTGCCCGGACTTGCCCATGCCGACCACGACGACCCGGccc        >  3:77034/1‑149 (MQ=255)
                                                                                                                                      catcatCACGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGGCGATCTTGTTGCCGATGTGCCCGGACTTGCCCATGCCGACCACGACGACCCGGccc        >  6:115117/1‑149 (MQ=255)
                                                                                                                                           tcaCGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGGCGATCTTGTTGCCGATGTGCCCGGACTTGCCCATGCCGACCACGACGACCCGGCCCTTGCt   >  5:330495/1‑149 (MQ=255)
                                                                                                                                             aCGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGGCGATCTTGTTGCCGATGTGCCCGGACTTGCCCATGCCGACCACGACGACCCGGCCCTTGCTg  <  3:294444/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
GACCTCGGCGGCCTGGGCCAGCGGCGAGTCCGGGTTGCCGGTCAGGCTGATCAGCTTGATGCCAAGGCGCTTGATGAACGGCAGCAGGGTGACGATTTCGGCGGTGCTGCCGGAGTTCGACAGCGCAAGGATGACATCATCACGGGTGATCATGCCCATGTCGCCATGGCTGGCTTCGGCCGGATGCACGAAAAACGCCGGTGTGCCGGTGCTGGCCAGGGTCGCGGCGATCTTGTTGCCGATGTGCCCGGACTTGCCCATGCCGACCACGACGACCCGGCCCTTGCTG  >  NC_002947/1099152‑1099440

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: ‑