Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,975,828 C→G 5.1% R254G (CGT→GGT)  pheA → bifunctional chorismate mutase/prephenate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,975,8280CG 5.1% 258.7 / 2.8 78R254G (CGT→GGT) pheAbifunctional chorismate mutase/prephenate dehydratase
Reads supporting (aligned to +/- strand):  ref base C (36/38);  new base G (2/2);  total (38/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.91e-01

CCGAAAGTGGCTGGACGCTCACTACCCGAACGTGGAGCGCGTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTG  >  NC_002947/1975686‑1975974
                                                                                                                                              |                                                                                                                                                  
ccGAAAGTGGCTGGACGCTCACTACCCGAACGTGGAGCGCGTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTgg                                                                                                                                              >  3:117279/1‑149 (MQ=255)
 cGAAAGTGGCTGGACGCTCACTACCCGAACGTGGAGCGCGTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGc                                                                                                                                             <  8:17775/149‑1 (MQ=255)
  gAAAGTGGCTGGACGCTCACTACCCGAACGTGGAGCGCGTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGcc                                                                                                                                            >  5:332977/1‑149 (MQ=255)
    aaGTGGCTGGACGCTCACTACCCGAACGTGGAGCGCGTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCga                                                                                                                                          <  7:213188/149‑1 (MQ=255)
          cTGGACGCTCACTACCCGAACGTGGAGCGCGTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGAt                                                                                                                                    <  6:306959/149‑1 (MQ=255)
             gACGCTCACTACCCGAACGTGGAGCGCGTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATc                                                                                                                                   <  8:241844/147‑1 (MQ=255)
             gACGCTCACTACCCGAACGTGGAGCGCGTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCg                                                                                                                                  <  4:117279/148‑1 (MQ=255)
              aCGCTCACTACCCGAACGTGGAGCGCGTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCg                                                                                                                                  >  1:336799/1‑147 (MQ=255)
              aCGCTCACTACCCGAACGTGGAGCGCGTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCg                                                                                                                                  >  1:342180/1‑147 (MQ=255)
              aCGCTCACTACCCGAACGTGGAGCGCGTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGaa                                                                                                                                <  5:74327/149‑1 (MQ=255)
                gCTCACTACCCGAACGTGGAGCGCGTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGa                                                                                                                              >  3:292610/1‑149 (MQ=255)
                  tCACTACCCGAACGTGGAGCGCGTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGAcc                                                                                                                            >  4:254158/1‑149 (MQ=255)
                    aCTACCCGAACGTGGAGCGCGTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGc                                                                                                                          >  2:223304/1‑149 (MQ=255)
                      tACCCGAACGTGGAGCGCGTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGaa                                                                                                                                >  8:346603/1‑141 (MQ=255)
                      tACCCGAACGTGGAGCGCGTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGaa                                                                                                                                <  7:346603/141‑1 (MQ=255)
                        cccGAACGTGGAGCGCGTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCg                                                                                                                       <  3:254158/148‑1 (MQ=255)
                         ccGAACGTGGAGCGCGTGTCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGa                                                                                                                     <  3:9678/149‑1 (MQ=255)
                         ccGAACGTGGAGCGCGTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGa                                                                                                                     >  8:336027/1‑149 (MQ=255)
                           gAACGTGGAGCGCGTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCgag                                                                                                                                         <  4:74213/127‑1 (MQ=255)
                           gAACGTGGAGCGCGTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCgag                                                                                                                                         >  3:74213/1‑127 (MQ=255)
                                  gAGCGCGTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCa                                                                                                             <  4:292610/148‑1 (MQ=255)
                                  gAGCGCGTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCAc                                                                                                            >  7:188871/1‑149 (MQ=255)
                                  gAGCGCGTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGCAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCAc                                                                                                            <  7:304435/149‑1 (MQ=255)
                                    gcgcGTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCAcgc                                                                                                          <  6:212580/149‑1 (MQ=255)
                                    gcgcGTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCAcgc                                                                                                          >  3:332410/1‑149 (MQ=255)
                                        gTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGaa                                                                                                                                >  4:190249/1‑123 (MQ=255)
                                        gTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGaa                                                                                                                                <  3:190249/123‑1 (MQ=255)
                                          ggCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTcc                                                                                                    <  3:57955/149‑1 (MQ=255)
                                            ccGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCtg                                                                                                  <  1:206205/149‑1 (MQ=255)
                                              gTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCtgat                                                                                                >  6:147298/1‑149 (MQ=255)
                                               tttCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCtgat                                                                                                >  2:35511/1‑148 (MQ=255)
                                                  cGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTgatgat                                                                                             <  4:332410/148‑1 (MQ=255)
                                                       aaTGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTAcc                                                                                                    <  8:340371/136‑1 (MQ=255)
                                                       aaTGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTAcc                                                                                                    >  7:340371/1‑136 (MQ=255)
                                                        aTGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAAc                                                                                      >  1:26240/1‑149 (MQ=255)
                                                          gCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCg                                                                                           <  8:233512/142‑1 (MQ=255)
                                                          gCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCg                                                                                           >  7:233512/1‑142 (MQ=255)
                                                            cGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAgg                                                                                  <  3:201749/149‑1 (MQ=255)
                                                             gAAGCGGCCAAGCGGGTCAAGCGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTCTACCGACTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCCGACAACTCCACGCGCTTCCTGATGAACGGAAACCAgg                                                                                  >  4:106677/1‑148 (MQ=255)
                                                                    ccAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGaa                                                                                >  2:55643/1‑143 (MQ=255)
                                                                    ccAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGaa                                                                                <  1:55643/143‑1 (MQ=255)
                                                                       aGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGaa                                                                                                                                >  6:329198/1‑92 (MQ=255)
                                                                       aGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGaa                                                                                                                                <  5:329198/92‑1 (MQ=255)
                                                                       aGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTAccgccg                                                                       >  3:9711/1‑149 (MQ=255)
                                                                          gggTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAg                                                                                   >  3:2122/1‑134 (MQ=255)
                                                                          gggTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAg                                                                                   <  4:2122/134‑1 (MQ=255)
                                                                                 gggTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGcga                                                               <  2:14310/147‑1 (MQ=255)
                                                                                 gggTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGcga                                                               >  1:14310/1‑147 (MQ=255)
                                                                                   gTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCgacgac                                                           <  1:223304/149‑1 (MQ=255)
                                                                                      aGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACaa                                                         <  7:336027/148‑1 (MQ=255)
                                                                                       gTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGa                                                       >  8:166852/1‑149 (MQ=255)
                                                                                         ggAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGAcc                                                     >  6:310224/1‑149 (MQ=255)
                                                                                          gAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCt                                                    <  6:332977/149‑1 (MQ=255)
                                                                                           aaCTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTc                                                   >  3:222223/1‑149 (MQ=255)
                                                                                           aaCTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTc                                                   <  2:342180/149‑1 (MQ=255)
                                                                                           aaCTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTc                                                   <  2:336799/149‑1 (MQ=255)
                                                                                            aCTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCAcg                                                                                                           <  8:191397/92‑1 (MQ=255)
                                                                                            aCTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCAcg                                                                                                           >  7:191397/1‑92 (MQ=255)
                                                                                                ggcggcGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGa                                                                      <  3:347900/125‑1 (MQ=255)
                                                                                                ggcggcGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGa                                                                      >  4:347900/1‑125 (MQ=255)
                                                                                                       aTCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTc                                       >  3:24908/1‑149 (MQ=255)
                                                                                                              gCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGc                                <  1:35511/149‑1 (MQ=255)
                                                                                                              gCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGAAACGATGatc                                <  8:188871/149‑3 (MQ=255)
                                                                                                                gATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGcaa                              >  7:342383/1‑149 (MQ=255)
                                                                                                                     ggcggcCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCAcgc                                                                                                          >  2:110751/1‑68 (MQ=255)
                                                                                                                          ccAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTAccgc                                                                         >  8:150523/1‑96 (MQ=255)
                                                                                                                          ccAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTAccgc                                                                         <  7:150523/96‑1 (MQ=255)
                                                                                                                          ccAACCTGTACGGCATGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGTTGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGt                    <  3:106677/149‑1 (MQ=255)
                                                                                                                                ttctcGGCCAGACGGGTCAGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTg              >  2:323725/5‑149 (MQ=255)
                                                                                                                                ttctcGGCCAGACGGGTCAGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTg              >  2:318410/5‑149 (MQ=255)
                                                                                                                                  ctcGGCCAGACGGGTCAGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATTATCGTTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGc             <  1:323725/146‑1 (MQ=255)
                                                                                                                                  ctcGGCCAGACGGGTCAGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGc             <  1:318410/146‑1 (MQ=255)
                                                                                                                                     ggCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCAt           >  4:218104/1‑147 (MQ=255)
                                                                                                                                       ccTGACCCGTCTGGCCGAGAAGCTCGAACACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGt       <  4:230217/149‑1 (MQ=255)
                                                                                                                                        cTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGtt      <  8:342383/149‑1 (MQ=255)
                                                                                                                                           aCCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTg     <  2:26240/147‑1 (MQ=255)
                                                                                                                                           aCCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCt   >  4:94846/1‑149 (MQ=255)
                                                                                                                                           aCCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCt   >  1:56428/1‑149 (MQ=255)
                                                                                                                                            cccGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTg  <  5:310224/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
CCGAAAGTGGCTGGACGCTCACTACCCGAACGTGGAGCGCGTGGCCGTTTCGAGCAATGCCGAAGCGGCCAAGCGGGTCAAGGGTGAGTGGAACTCGGCGGCGATCGCCGGCGATATGGCGGCCAACCTGTACGGCCTGACCCGTCTGGCCGAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTG  >  NC_002947/1975686‑1975974

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