Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,975,987 C→T 7.1% Q307* (CAG→TAG)  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,9870CT 7.1% 205.3 / 4.7 70Q307* (CAG→TAG) pheAbifunctional chorismate mutase/prephenate dehydratase
Reads supporting (aligned to +/- strand):  ref base C (36/29);  new base T (2/3);  total (38/32)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.29e-01

GAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGGCGC  >  NC_002947/1975839‑1976134
                                                                                                                                                    |                                                                                                                                                    
gaagaTCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATc                                                                                                                                                      <  7:156601/149‑1 (MQ=255)
gaagaTCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATc                                                                                                                                                      <  7:152285/149‑1 (MQ=255)
    aTCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGaa                                                                                                                                                  >  2:214772/1‑149 (MQ=255)
       gAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACgg                                                                                                                                               >  5:193699/1‑149 (MQ=255)
       gAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACg                                                                                                                                                <  8:21709/148‑1 (MQ=255)
          gACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTGATAGAACGGCac                                                                                                                                            <  5:332215/149‑2 (MQ=255)
          gACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTGATAGAACGGCac                                                                                                                                            <  3:232689/149‑2 (MQ=255)
           aCCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATc                                                                                                                                           <  6:285361/149‑1 (MQ=255)
             cGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTGATAGAACGGCacc                                                                                                                                           >  4:72328/1‑145 (MQ=255)
             cGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTGATAGAACGGCacc                                                                                                                                           <  3:72328/147‑3 (MQ=255)
             cGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGa                                                                                                                                         <  7:185942/149‑1 (MQ=255)
              gCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGAc                                                                                                                                        <  2:128680/149‑1 (MQ=255)
                 cGGCCAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGAGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCAAGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTg                                                                                                                                     >  8:82622/1‑149 (MQ=255)
                       aCTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGc                                                                                                                              <  6:198886/149‑1 (MQ=255)
                              gcgcTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCgaga                                                                                                                       >  2:152619/1‑149 (MQ=255)
                                 cTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGAccc                                                                                                                    >  1:232054/1‑149 (MQ=255)
                                  ttCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGAcccg                                                                                                                   >  3:258717/1‑149 (MQ=255)
                                    ccTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGAcccgc                                                                                                                  <  6:193699/148‑1 (MQ=255)
                                      tgatgaTCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGAcccgcccg                                                                                                               >  4:14703/1‑149 (MQ=255)
                                       gatgatCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGt                                                                                                              >  8:230028/1‑149 (MQ=255)
                                       gatgatCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGt                                                                                                              >  1:97321/1‑149 (MQ=255)
                                                  aCCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGc                                                                                                   <  1:193258/149‑1 (MQ=255)
                                                    cAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGCCCCGCATCGAGACCCGCCCGTCGCGCAGCGGc                                                                                                   >  5:341392/1‑147 (MQ=255)
                                                    cAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGc                                                                                                   >  1:124508/1‑147 (MQ=255)
                                                    cAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGc                                                                                                   >  4:52127/1‑147 (MQ=255)
                                                    cAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGc                                                                                                   >  1:154700/1‑147 (MQ=255)
                                                      ggAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGt                                                                                               <  3:14703/149‑1 (MQ=255)
                                                      ggAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGt                                                                                               <  3:52127/149‑1 (MQ=255)
                                                        aaGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTgg                                                                                             <  7:230028/149‑1 (MQ=255)
                                                          gTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGAc                                                                                           <  2:232054/149‑1 (MQ=255)
                                                           tACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGTGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGAc                                                                                           <  2:290292/148‑1 (MQ=255)
                                                             ccgccgACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCt                                                                                         >  2:264235/1‑148 (MQ=255)
                                                                    ccGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGttc                                                                                 <  1:264235/149‑1 (MQ=255)
                                                                       gCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGttctt                                                                               >  7:50384/1‑148 (MQ=255)
                                                                          acgacAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATc                                                                           <  2:124508/149‑1 (MQ=255)
                                                                              cAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGAc                                                                                           >  7:13814/1‑129 (MQ=255)
                                                                              cAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGAc                                                                                           <  8:13814/129‑1 (MQ=255)
                                                                              cAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACt                                                                       >  5:173198/1‑149 (MQ=255)
                                                                               aaGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCaa                                                                                                 <  5:164987/122‑1 (MQ=255)
                                                                               aaGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCaa                                                                                                 >  6:164987/1‑122 (MQ=255)
                                                                                    cTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGa                                                                         >  2:318385/1‑141 (MQ=255)
                                                                                    cTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGa                                                                         <  1:318385/141‑1 (MQ=255)
                                                                                               tttCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACcgc                                                      >  7:2441/1‑149 (MQ=255)
                                                                                                   gATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGAcc                                                  >  1:13403/1‑149 (MQ=255)
                                                                                                      gAGCAACAAGAAAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGTCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCAACGCGACCCGc                                               >  5:74270/1‑149 (MQ=255)
                                                                                                        gCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCt                                              <  8:50384/148‑1 (MQ=255)
                                                                                                          aacaaGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACgt                                                                                     >  7:154869/1‑107 (MQ=255)
                                                                                                          aacaaGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACgt                                                                                     <  8:154869/107‑1 (MQ=255)
                                                                                                           acaaGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGAt                                           >  6:273423/1‑148 (MQ=255)
                                                                                                           acaaGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCCGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGAt                                           >  3:116222/1‑148 (MQ=255)
                                                                                                             aaGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGAGTCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGAta                                          >  1:99212/1‑147 (MQ=255)
                                                                                                               gCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAAAGGCATCGAACAGACCCGCATCGAGACCCGCCCGTCGAGCACCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGTCCACCACCGCGACCCGCTGATCAAgg                                      >  4:76097/1‑149 (MQ=255)
                                                                                                                  aGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACCTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCgg                                   >  7:156715/1‑149 (MQ=255)
                                                                                                                         ttGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGc                                               <  3:290781/130‑1 (MQ=255)
                                                                                                                         ttGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGc                                               >  4:290781/1‑130 (MQ=255)
                                                                                                                            cATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAg                           >  6:182511/1‑147 (MQ=255)
                                                                                                                              tGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGa                       <  1:152619/149‑1 (MQ=255)
                                                                                                                                aGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGAt                      <  7:82622/148‑1 (MQ=255)
                                                                                                                                  ttgagggTGCCGTTCTATCACAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCa                    <  3:172848/143‑1 (MQ=255)
                                                                                                                                    gCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGc                  >  4:93602/1‑148 (MQ=255)
                                                                                                                                       gatGCCGTTCTGATAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCag               >  8:185292/3‑148 (MQ=255)
                                                                                                                                        gTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCagga             <  2:154700/149‑1 (MQ=255)
                                                                                                                                        gTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCag               >  4:25376/1‑147 (MQ=255)
                                                                                                                                          gCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCagg              >  2:55163/1‑146 (MQ=255)
                                                                                                                                          gCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCagg              <  1:55163/146‑1 (MQ=255)
                                                                                                                                           ccGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGc          >  7:202219/1‑149 (MQ=255)
                                                                                                                                                cTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCtgtg      <  2:99212/148‑1 (MQ=255)
                                                                                                                                                cTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCtgtg      <  4:258717/148‑1 (MQ=255)
                                                                                                                                                  aTCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGgcg   >  8:773/1‑149 (MQ=255)
                                                                                                                                                  aTCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGgcg   >  4:43561/1‑149 (MQ=255)
                                                                                                                                                   tCAGAACGTCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGgcgc  <  6:74270/149‑1 (MQ=255)
                                                                                                                                                   tCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGgcg   <  8:2441/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGGCGC  >  NC_002947/1975839‑1976134

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: