Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,975,985 A→G 6.5% 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 6.5% 179.1 / 2.5 61Y306C (TAT→TGT) ‡pheAbifunctional chorismate mutase/prephenate dehydratase
Reads supporting (aligned to +/- strand):  ref base A (31/26);  new base G (2/2);  total (33/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.37e-01

GAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATC‑‑ATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGGCG  >  NC_002947/1975837‑1976133
                                                                                                                                                      |                                                                                                                                                    
gagaAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATC‑‑ATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTa                                                                                                                                                      <  1:267830/149‑1 (MQ=255)
 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                                                                                                                            >  5:319332/1‑149 (MQ=255)
                          cTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATC‑‑ATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCa                                                                                                                            >  3:225982/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                                                                       >  8:41287/1‑147 (MQ=255)
                                                                                 aaGACCTCGATC‑‑ATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGAc                                                                       <  7:41287/147‑1 (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:170700/148‑1 (MQ=255)
                                                                                     ccTCGATC‑‑ATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTtcgt                                                                  >  8:170700/1‑148 (MQ=255)
                                                                                     ccTCGATC‑‑ATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTtcgt                                                                  <  7:239669/148‑1 (MQ=255)
                                                                                      cTCGATC‑‑ATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTtcgtcg                                                                >  7:156741/1‑149 (MQ=255)
                                                                                         gATC‑‑ATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGcc                                                             <  6:319332/149‑1 (MQ=255)
                                                                                         gATC‑‑ATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGcc                                                             >  8:43543/1‑149 (MQ=255)
                                                                                                 cGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAg                                                                                                      <  7:346058/102‑1 (MQ=255)
                                                                                                 cGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAg                                                                                                      >  8:346058/1‑102 (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              >  8:130664/2‑147 (MQ=255)
                                                                                                                                            atGCCGTTCTGATAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCag              >  7:250443/2‑147 (MQ=255)
                                                                                                                                              gCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCCTCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCaggagg          >  5:52723/1‑149 (MQ=255)
                                                                                                                                              gCCGGTCTATCAGAACGGCAGCGCCCGGACCCGCCTCGAGACGCGCCCGGCGCGCAGCGGCAAGTGGACCTGCGGGTTGTTAATCGCCGTCGTGGGGCCCCCCCCCGCCCCGCGGAGCACGGCGGGGGTGGCGCCGATCAGCCaggcgg          >  6:359162/1‑149 (MQ=255)
                                                                                                                                               ccGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGc                                              <  2:371102/112‑1 (MQ=255)
                                                                                                                                               ccGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGc                                              >  1:371102/1‑112 (MQ=255)
                                                                                                                                                      aTCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGgcg  >  6:51094/1‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
GAGAAGATCGAAGACCGCCCGGACAACTCCACGCGCTTCCTGATGATCGGTAACCAGGAAGTACCGCCGACCGGCGACGACAAGACCTCGATC‑‑ATCGTTTCGATGAGCAACAAGCCAGGTGCGTTGCATGAGTTGCTGGTGCCGTTCTATCAGAACGGCATCGACCTGACCCGCATCGAGACCCGCCCGTCGCGCAGCGGCAAGTGGACCTACGTGTTCTTCATCGACTTCGTCGGCCACCACCGCGACCCGCTGATCAAGGCGGTGCTCGAGCAGATCAGCCAGGAGGCTGTGGCG  >  NC_002947/1975837‑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: