Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,975,985 A→G 7.3% Y306C (TAT→TGT) ‡ 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,9850AG 7.3% 194.7 / 5.1 68Y306C (TAT→TGT) ‡pheAbifunctional chorismate mutase/prephenate dehydratase
Reads supporting (aligned to +/- strand):  ref base A (34/29);  new base G (2/3);  total (36/32)
Fisher's exact test for biased strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

GAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGGCG  >  NC_002947/1975839‑1976133
                                                                                                                                                  |                                                                                                                                                     
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                                                                                                                                           <  3:232689/149‑2 (MQ=255)
          gACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTGATAGAACGGCac                                                                                                                                           <  5:332215/149‑2 (MQ=255)
           aCCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATc                                                                                                                                          <  6:285361/149‑1 (MQ=255)
             cGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTGATAGAACGGCacc                                                                                                                                          <  3:72328/147‑3 (MQ=255)
             cGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTGATAGAACGGCacc                                                                                                                                          >  4:72328/1‑145 (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                                                                                                             >  1:97321/1‑149 (MQ=255)
                                       gatgatCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGt                                                                                                             >  8:230028/1‑149 (MQ=255)
                                                  aCCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGc                                                                                                  <  1:193258/149‑1 (MQ=255)
                                                    cAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGCCCCGCATCGAGACCCGCCCGTCGCGCAGCGGc                                                                                                  >  5:341392/1‑147 (MQ=255)
                                                    cAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGc                                                                                                  >  4:52127/1‑147 (MQ=255)
                                                    cAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGc                                                                                                  >  1:124508/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                                                                                          <  8:13814/129‑1 (MQ=255)
                                                                              cAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGAc                                                                                          >  7:13814/1‑129 (MQ=255)
                                                                              cAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACt                                                                      >  5:173198/1‑149 (MQ=255)
                                                                               aaGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCaa                                                                                                >  6:164987/1‑122 (MQ=255)
                                                                               aaGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCaa                                                                                                <  5:164987/122‑1 (MQ=255)
                                                                                    cTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGa                                                                        <  1:318385/141‑1 (MQ=255)
                                                                                    cTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGa                                                                        >  2:318385/1‑141 (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             <  1:55163/146‑1 (MQ=255)
                                                                                                                                          gCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCagg             >  2:55163/1‑146 (MQ=255)
                                                                                                                                           ccGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGc         >  7:202219/1‑149 (MQ=255)
                                                                                                                                                cTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCtgtg     <  2:99212/148‑1 (MQ=255)
                                                                                                                                                cTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCtgtg     <  4:258717/148‑1 (MQ=255)
                                                                                                                                                  aTCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGgcg  >  4:43561/1‑149 (MQ=255)
                                                                                                                                                  aTCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGgcg  >  8:773/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGGCG  >  NC_002947/1975839‑1976133

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: