Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,975,986 T→A 7.7% Y306* (TAT→TAA) ‡ 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,9860TA 7.7% 173.0 / 3.1 61Y306* (TAT→TAA) ‡pheAbifunctional chorismate mutase/prephenate dehydratase
Reads supporting (aligned to +/- strand):  ref base T (31/25);  new base A (3/2);  total (34/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

AGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATC‑‑ATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGGCGC  >  NC_002947/1975838‑1976134
                                                                                                                                                      |                                                                                                                                                    
agaagaTCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATC‑‑ATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTAt                                                                                                                                                      <  6:352912/149‑1 (MQ=255)
   agaTCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATC‑‑ATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCa                                                                                                                                                    <  3:367154/148‑1 (MQ=255)
      tCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATC‑‑ATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAAc                                                                                                                                                >  2:42093/1‑149 (MQ=255)
      tCGAAGACCGCCCGGACAACTCCACGCGCTTCCAGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATC‑‑ATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAAc                                                                                                                                                >  1:50390/1‑149 (MQ=255)
        gAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACGTCGATC‑‑ATCGTTTCGATGAGCAACAAGCCTGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACgg                                                                                                                                              <  6:276067/149‑1 (MQ=255)
           gACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATC‑‑ATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTGATAGAACGGCac                                                                                                                                           <  2:11643/149‑2 (MQ=255)
             ccGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATC‑‑ATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTGATAGAACGGCacca                                                                                                                                         <  1:347214/149‑4 (MQ=255)
                cccGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATC‑‑ATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGAcc                                                                                                                                      <  6:58817/149‑1 (MQ=255)
                    gACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATC‑‑ATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGAc                                                                                                                                  <  5:324788/149‑1 (MQ=255)
                         cTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATC‑‑ATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCa                                                                                                                             >  8:79809/1‑149 (MQ=255)
                         cTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATC‑‑ATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCa                                                                                                                             >  3:225982/1‑149 (MQ=255)
                         cTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATC‑‑ATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCa                                                                                                                             >  5:319332/1‑149 (MQ=255)
                         cTCCACGCGCTTCATGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACATCGATC‑‑ATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGc                                                                                                                              >  2:357752/1‑148 (MQ=255)
                                 gcTTCCTGATGATCGGTAACCAGGAAGTACCTCCGACCGGCGACGACAACAAATCGATC‑‑ATACTTTCGATGAGCAACAAGCCACGTTCGATGCATGAGTTGCTGGTGCACTTATAACAGAACGGCATCGACCTGCCCCGCATCGAGAcc                                                                                                                     >  8:239669/1‑149 (MQ=255)
                                 gcTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATC‑‑ATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGAcc                                                                                                                     >  2:251072/1‑149 (MQ=255)
                                   ttCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATCTGATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGTTGCCGTTCTATCAGAACGGCATCGACATGCCCCGCATCGAGAcc                                                                                                                     >  1:102360/1‑149 (MQ=255)
                                    tCCTGATGATTGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATC‑‑ATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGAcccgc                                                                                                                  >  4:282026/1‑149 (MQ=255)
                                       tgatgaTCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATC‑‑ATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGAcccgcccg                                                                                                               >  3:191440/1‑149 (MQ=255)
                                             tCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATC‑‑ATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCgcgc                                                                                                         >  5:246940/1‑149 (MQ=255)
                                               ggTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATC‑‑ATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAg                                                                                                       <  8:178716/149‑1 (MQ=255)
                                                    ccAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATC‑‑ATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGc                                                                                                   >  1:343880/1‑148 (MQ=255)
                                                     cAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATC‑‑ATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTc                                                                                                             >  7:165588/1‑137 (MQ=255)
                                                     cAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATC‑‑ATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTc                                                                                                             <  8:165588/137‑1 (MQ=255)
                                                               cgccgACCGGCGACGACAAGACCTCGATC‑‑ATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTAc                                                                                       <  1:42093/149‑1 (MQ=255)
                                                                gccgACCGGCGACGACAAGACCTCGATC‑‑ATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACg                                                                                      >  5:203509/1‑149 (MQ=255)
                                                                      cGGCGACGACAAGACCTCGATC‑‑ATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGttct                                                                                >  3:74669/1‑149 (MQ=255)
                                                                           acgacAAGACCTCGATC‑‑ATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATc                                                                           <  2:226735/149‑1 (MQ=255)
                                                                               cAAGACCTCGATC‑‑ATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACt                                                                       >  1:354033/1‑149 (MQ=255)
                                                                                aaGACCTCGATC‑‑ATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGAc                                                                        <  7:41287/147‑1 (MQ=255)
                                                                                aaGACCTCGATC‑‑ATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGAc                                                                        >  8:41287/1‑147 (MQ=255)
                                                                                aaGACCTCGATC‑‑ATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACtt                                                                      >  3:321442/1‑149 (MQ=255)
                                                                                   aCCTCGATC‑‑ATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTtcgt                                                                   <  2:102360/149‑1 (MQ=255)
                                                                                    ccTCGATC‑‑ATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTtcgtc                                                                  <  3:282026/149‑1 (MQ=255)
                                                                                    ccTCGATC‑‑ATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTtcgt                                                                   <  7:239669/148‑1 (MQ=255)
                                                                                    ccTCGATC‑‑ATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTtcgt                                                                   >  8:170700/1‑148 (MQ=255)
                                                                                    ccTCGATC‑‑ATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTtcgt                                                                   <  7:170700/148‑1 (MQ=255)
                                                                                     cTCGATC‑‑ATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTtcgtcg                                                                 >  7:156741/1‑149 (MQ=255)
                                                                                        gATC‑‑ATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGcc                                                              >  8:43543/1‑149 (MQ=255)
                                                                                        gATC‑‑ATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGcc                                                              <  6:319332/149‑1 (MQ=255)
                                                                                                cGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAg                                                                                                       >  8:346058/1‑102 (MQ=255)
                                                                                                cGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAg                                                                                                       <  7:346058/102‑1 (MQ=255)
                                                                                                        tGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCg                                                <  6:203509/149‑1 (MQ=255)
                                                                                                             aacaaGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGAt                                           >  4:244581/1‑149 (MQ=255)
                                                                                                              acaaGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATc                                          <  1:251072/149‑1 (MQ=255)
                                                                                                              acaaGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATc                                          <  4:74669/149‑1 (MQ=255)
                                                                                                                      ggTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGt                                  <  7:79809/149‑1 (MQ=255)
                                                                                                                            ttGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGTCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGa                            <  1:357752/149‑1 (MQ=255)
                                                                                                                            ttGCATGAGTTGCTGGGGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGa                            <  6:246940/149‑1 (MQ=255)
                                                                                                                              gCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAg                           >  2:15280/1‑148 (MQ=255)
                                                                                                                               cATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGc                          <  4:225982/148‑1 (MQ=255)
                                                                                                                                  gAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGa                       >  6:269329/1‑148 (MQ=255)
                                                                                                                                    gTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCa                    <  7:43543/149‑1 (MQ=255)
                                                                                                                                       gCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGc                  <  4:191440/148‑1 (MQ=255)
                                                                                                                                         tGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCag               >  4:267342/1‑149 (MQ=255)
                                                                                                                                          ggTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCag               >  2:247903/1‑148 (MQ=255)
                                                                                                                                           atGCCGTTCTGATAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCag               >  7:250443/2‑147 (MQ=255)
                                                                                                                                           atGCCGTTCTGATAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCag               >  8:130664/2‑147 (MQ=255)
                                                                                                                                             gCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCCTCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCaggagg           >  5:52723/1‑149 (MQ=255)
                                                                                                                                             gCCGGTCTATCAGAACGGCAGCGCCCGGACCCGCCTCGAGACGCGCCCGGCGCGCAGCGGCAAGTGGACCTGCGGGTTGTTAATCGCCGTCGTGGGGCCCCCCCCCGCCCCGCGGAGCACGGCGGGGGTGGCGCCGATCAGCCaggcgg           >  6:359162/1‑149 (MQ=255)
                                                                                                                                              ccGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGc                                               >  1:371102/1‑112 (MQ=255)
                                                                                                                                              ccGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGc                                               <  2:371102/112‑1 (MQ=255)
                                                                                                                                                     aTCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGgcg   >  6:51094/1‑149 (MQ=255)
                                                                                                                                                      tCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGgcgc  >  2:184286/1‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
AGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATC‑‑ATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGGCGC  >  NC_002947/1975838‑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: