Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,471,553 G→C 5.5% P425R (CCG→CGG)  alg8 ← glycosyltransferase Alg8

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,471,5530GC 5.5% 227.0 / 2.2 71P425R (CCG→CGG) alg8glycosyltransferase Alg8
Reads supporting (aligned to +/- strand):  ref base G (29/38);  new base C (2/2);  total (31/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

TGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCAGCGGGTAGGCCGGGCCTACGTTGTGGCCGGAGCACAGCAGCATGATGGTGAGGATCAGGCGGGTAATGCCGATCCACAGCAGGTACACCAGCAAGAAGGCCATGCCGAACTTGAGGCTGGCGATTACCGCCACGGTCAGGCCCAGCAGGCTGGTCCACA  >  NC_002947/1471407‑1471700
                                                                                                                                                  |                                                                                                                                                   
tgtTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCAGCGGGTAGGCCggg                                                                                                                                                   <  6:344825/149‑1 (MQ=255)
    gAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCAGCGGGTAGGCCGGGCCTa                                                                                                                                               >  4:205690/1‑149 (MQ=255)
      aCCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCAGCGGGTAGGCCGGGCCTACg                                                                                                                                             >  6:47201/1‑149 (MQ=255)
        cATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCAGCGGGTAGGCCGGGCCTACGtt                                                                                                                                           >  6:197374/1‑149 (MQ=255)
         tttgttgAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCAGCGGGTAGGCCGGGCCTACGTtg                                                                                                                                          >  2:56645/2‑149 (MQ=255)
             ttgAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCAGCGGGTAGGCCGGGCCTACGTTGTGGc                                                                                                                                      >  1:100251/1‑149 (MQ=255)
                aaaGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATTAGCGCGCCGATGATCTGGTTGTAATAGAGAATCAGCGGGTAGGCCGGGCCTACGTTGTGGCCg                                                                                                                                    <  8:214405/148‑1 (MQ=255)
                aaaGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCAGCGGGTAGGCCGGGCCTACGTTGTGGCCg                                                                                                                                    >  4:48856/1‑148 (MQ=255)
                 aaGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCAGCGGGTAGGCCGGGCCTACGTTGTGGCCGGa                                                                                                                                  <  2:203704/149‑1 (MQ=255)
                  aGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCAGCGGGTAGGCCGGGCCTACGTTGTGGCCGGAg                                                                                                                                 <  7:145830/149‑1 (MQ=255)
                        cGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCAGCGGGTAGGCCGGGCCTACGTTGTGGCCGGAGCAcagc                                                                                                                           <  6:29517/149‑1 (MQ=255)
                        cGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCAGCGGGTAGGCCGGGCCTACGTTGTGGCCGGAGCAcagc                                                                                                                           <  6:24245/149‑1 (MQ=255)
                              cACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCAGCGGGTAGGCCGGGCCTACGTTGTGGCCGGAGCACAGCAGCat                                                                                                                      >  7:338619/1‑148 (MQ=255)
                               aCGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCAGCGGGTAGGCCGGGCCTACGTTGTGGCCGGAGCACAGCAGCatga                                                                                                                    <  5:197374/149‑1 (MQ=255)
                               aCGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCAGCGGGTAGGCCGGGCCTACGTTGTGGCCGGAGCACAGCAGCatga                                                                                                                    >  2:316500/1‑149 (MQ=255)
                               aCGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGAGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCAGCGGGTAGGCCGGGCCTACGTTGTGGCCGGAGCACAGCAGCatga                                                                                                                    >  3:122817/1‑149 (MQ=255)
                                     gAGGGCAGTGGGCTGAAGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCAGCGGGTAGGCCGGGCCTACGTTGTGGCCGGAGCACAGCAGCATGATGGTg                                                                                                               <  8:320956/148‑1 (MQ=255)
                                       gggCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCAGCGGGTAGGCCGGGCCTACGTTGTGGCCGGAGCACAGCAGCATGATGGTGAg                                                                                                             <  5:47201/148‑1 (MQ=255)
                                                  tGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCAGCGGGTAGGCCGGGCCTACGTTGTGGCCGGAGCACAGCAGCATGATGGTGAGGATCAGGCggg                                                                                                  >  6:162/1‑148 (MQ=255)
                                                  tGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCAGCGGGTAGGCCGGGCCTACGTTGTGGCCGGAGCACAGCAGCATGATGGTGAGGATCAGGCggg                                                                                                  >  5:350805/1‑148 (MQ=255)
                                                  tGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCAGCGGGTAGGCCGGGCCTACGTTGTGGCCGGAGCACAGCAGCATGATGGTGAGGATCAGGCggg                                                                                                  >  4:3800/1‑148 (MQ=255)
                                                    aCGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCAGCGGGTAGGCCGGGCCTACGTTGTGGCCGGAGCACAGCAGCATGATGGTGAGGATCAGGCGGGTaa                                                                                               >  7:71097/1‑149 (MQ=255)
                                                      gctgccACGAGTGCTTGTCGATGCGGAAGAACACGTAGATCTTCATCAGCGGTCCGTTGATCTGGTTGTAATAGAGAATCAGCGGGTAGGCCGGGCCTACGTTGTGGCCGGAGCACAGCAGCAGGATGGTGAGGATCAGGCGGGTAATg                                                                                             <  4:87368/145‑1 (MQ=255)
                                                      gcgTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCGGCGCGCCGATGATCTGGTTGTAATAGAGAATCAGCGGGTAGGCCGGGCCTACGTTGTGGCCGGAGCACAGCAGCATGATGGTGa                                                                                                              <  2:40641/132‑1 (MQ=255)
                                                           cACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCAGCGGGTAGGCCGGGCCTACGTTGTGGCCGGAGCACAGCAGCATGATGGTGAGGATCAGGCGGGTAATGCCGAt                                                                                        <  5:162/149‑1 (MQ=255)
                                                           cACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCAGCGGGTAGGCCGGGCCTACGTTGTGGCCGGAGCACAGCAGCATGATGGTGAGGATCAGGCGGGTAATGCCGAt                                                                                        <  4:122817/149‑1 (MQ=255)
                                                            aCGACTGCTTGTCGAGTAGGAATAACACGTAGATCTTCATCAGCGCTACGATGATCTGGTTGTAATAGAGAATCAGCGGGTAGTCCGGGCCTACGTTGTGGCCGGAGCACAGCAGCATGATGGTGAGGATCAGGCGGGTAATGCCGaac                                                                                       <  1:80430/149‑3 (MQ=255)
                                                            aCGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCAGCGGGTAGGCCGGGCCTACGTTGTGGCCGGAGCACAGCAGCATGATGGTGAGGATCAGGCgg                                                                                                   <  7:307141/137‑1 (MQ=255)
                                                            aCGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCAGCGGGTAGGCCGGGCCTACGTTGTGGCCGGAGCACAGCAGCATGATGGTGAGGATCAGGCgg                                                                                                   >  8:307141/1‑137 (MQ=255)
                                                            aCGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCAGCGGGTAGGCCGGGCCTACGTTGTGGCCGGAGCACAGCAGCATGATGGTGAGGATCAGGCGGGTAATGCCGATc                                                                                       >  4:333317/1‑149 (MQ=255)
                                                            aCGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCAGCGGGTAGGCCGGGCCTACGTTGTGGCCGGAGCACAGCAGCATGATGGTGAGGATCAGGCGGGTAATGCCGATc                                                                                       >  6:318399/1‑149 (MQ=255)
                                                                cTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCAGCGGGTAGGCCGGGCCTACGTTGTGGCCGGAGCACAGCAGCATGATGGTGAGGATCAGGCGGGTAATGCCGATCcaca                                                                                   <  5:303130/149‑1 (MQ=255)
                                                                cTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCAGCGGGTAGGCCGGGCCTACGTTGTGGCCGGAGCACAGCAGCATGATGGTGAGGATCAGGCGGGTAATGCCGATCca                                                                                     >  8:26105/1‑147 (MQ=255)
                                                                 tGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCAGCGGGTAGGCCGGGCCTACGTTGTGGCCGGAGCACAGCAGCATGATGGTGAGGATCAGGCGGGTAATGCCGATCCAcag                                                                                  >  5:335164/1‑149 (MQ=255)
                                                                   cTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCAGCGGGTAGGCCGGGCCTACGTTGTGGCCGGAGCACAGCAGCATGATGGTGa                                                                                                              <  6:72253/119‑1 (MQ=255)
                                                                   cTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCAGCGGGTAGGCCGGGCCTACGTTGTGGCCGGAGCACAGCAGCATGATGGTGa                                                                                                              >  5:72253/1‑119 (MQ=255)
                                                                   cTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCAGCGGGTAGGCCGGGCCTACGTTGTGGCCGGAGCACAGCAGCATGATGGTGAGGATCAGGCGGGTAATGCCGATCCAcagca                                                                                <  3:333317/149‑1 (MQ=255)
                                                                   cTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCAGCGGGTAGGCCGGGCCTACGTTGTGGCCGGAGCACAGCAGCATGATGGTGAGGATCAGGCGGGTAATGCCGATCCAcag                                                                                  <  1:316500/147‑1 (MQ=255)
                                                                      gTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCAGCGGGTAGGCCGGGCCTACGTTGTGGCCGGAGCACAGCAGCATGATGGTGAGGATCAGGCGGGTAATGCCGATCCAcagcag                                                                               <  6:350805/147‑1 (MQ=255)
                                                                       tCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCAGCGGGTAGGCCGGGCCTACGTTGTGGCCGGAGCACAGCAGCATGATGGTGAGGATCAGGCGGGTAATGCCGa                                                                                         <  2:274600/136‑1 (MQ=255)
                                                                       tCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCAGCGGGTAGGCCGGGCCTACGTTGTGGCCGGAGCACAGCAGCATGATGGTGAGGATCAGGCGGGTAATGCCGa                                                                                         >  1:274600/1‑136 (MQ=255)
                                                                        cGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCAGCGTGTAGGCCGGGCCTACGTTGTGGCCGGAGCACAGCAGCATGATGGTGAGGATCAGTCGGGTAATGCCGATCCACAGCAGGTa                                                                            <  5:280403/148‑1 (MQ=255)
                                                                        cGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCAGCGGGTAGGCCGGGCCTACGTTGTGGCCGGAGCACAGCAGCATGATGGTGAGGATCAGGCGGGTAATGCCGATCCACAGCAGGTa                                                                            >  3:162937/1‑148 (MQ=255)
                                                                         gAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCAGCGGGTAGGCCGGGCCTACGTTGTGGCCGGAGCACAGCAGCATGATGGTGAGGATCAGGCGGGTAATGCCGATCCACAGCAGGTaca                                                                          >  1:152178/1‑149 (MQ=255)
                                                                          aGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCAGCGGGTAGGCCGGGCCTACGTTGTGGCCGGAGCACAGCAGCATGATGGTGAGGATCAGGCGGGTAATGCCGATCCACAGCAGGTacac                                                                         <  3:48856/149‑1 (MQ=255)
                                                                                     cacGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCAGCGGGTAGGCCGGGCCTACGTTGTGGCCGGAGCACAGCAGCATGATGGTGAGGATCAGGCGGGTAATGCCGATCCACAGCAGGTACACCAGCAAGAAgg                                                              <  8:338619/149‑1 (MQ=255)
                                                                                      ctgaaggtCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCAGCGGGTAGGCCGGGCCTACGTTGTGGCCGGAGCACAGCAGCATGATGGTGAGGATCAGGCGGGTAATGCCGATCCACAGCAGGTACACCAGCAAGAAgg                                                              <  8:70193/141‑1 (MQ=255)
                                                                                                 catcaGCGCGCCGATGATCTGGTTGTAATAGAGAATCAGCGGGTAGGCCGGGCCTACGTTGTGGCCGGAGCACAGCAGCATGATGGTGAGGATCAGGCGGGTAATGCCGATCCACAGCAGGTACACCAGCAAGAAGGCCATGCCGAAc                                                   <  3:205690/148‑1 (MQ=255)
                                                                                                     aGCGCGCCGATGATCTGGTTGTAATAGAGAATCAGCGGGTAGGCCGGGCCTACGTTGTGGCCGGAGCACAGCAGCATGATGGTGAGGATCAGGCGGGTAATGCCGATCCACAGCAGGTACACCAGCAAGAAGGCCATGCCGAACTTGa                                               <  8:71097/148‑1 (MQ=255)
                                                                                                       cgcgcCGATGATCTGGTTGTAATAGAGAATCAGCGGGTAGGCCGGGCCTACGTTGTGGCCGGAGCACAGCAGCATGATGGTGAGGATCAGGCGGGTAATGCCGATCCACAGCAGGTACACCAGCAAGAAGGCCATGCCGAACTTGAGGc                                            <  6:335164/149‑1 (MQ=255)
                                                                                                           ccGATGATCTGGTTGTAATAGAGAATCAGCGGGTAGGCCGGGCCTACGTTGTGGCCGGAGCACAGCAGCATGATGGTGAGGATCAGGCGGGTAATGCCGATCCACAGCAGGTACACCAGCAAGAAGGCCATGCCGAACTTGAGGCTGGc                                        >  6:321683/1‑149 (MQ=255)
                                                                                                            cGATGATCTGGTTGTAATAGAGAATCAGCGGGTAGGCCGGGCCTACGTTGTGGCCGGAGCACAGCAGCATGATGGTGAGGATCAGGCGGGTAATGCCGATCCACAGCAGGTACACCAGCAAGAAGGCCATGCCGAACTTGAGGCTGGc                                        >  6:36933/1‑148 (MQ=255)
                                                                                                            cGATGATCTGGTTGTAATAGAGAATCAGCGGGTAGGCCGGGCCTACGTTGTGGCCGGAGCACAGCAGCATGATGGTGAGGATCAGGCGGGTAATGCCGATCCACAGCAGGCACACCAGCAAGAAGGCCATGCCGAACTTGAGGCTGGc                                        >  5:113268/1‑148 (MQ=255)
                                                                                                             gatgatCTGGTTGTAATAGAGAATCAGCGGGTAGGCCGGGCCTACGTTGTGGCCGGAGCACAGCAGCATGATGGTGAGGATCAGGCGGGTAATGCCGATCCACAGCAGGCACACCAGCAAGAAGGCCATGCCGAACTTGAGGCTGGCGa                                      <  6:113268/149‑1 (MQ=255)
                                                                                                                   cTGGTTGTAATAGAGAATCAGCGGGTAGGCCGGGCCTACGTTGTGGCCGGAGCACAGCAGCATGATGGTGAGGATCAGGCGGGTAATGCCGATCCACAGCAGGTACACCAGCAAGAAGGCCATGCCGAACTTGAGGCTGGCGATTACCg                                >  7:158544/1‑149 (MQ=255)
                                                                                                                         gTAATAGAGAATCAGCGGGTAGGCCGGGCCTACGTTGTGGCCGGAGCACAGCAGCATGATGGTGAGGATCAGGCGGGTAATGCCGATCCACAGCAGGTACACCAGCAAGAAGGCCATGCCGAACTTGAGGCTGGCGATTACCGCCACgg                          >  2:46506/1‑149 (MQ=255)
                                                                                                                          tAATAGAGAATCAGCGGGTAGGCCGGGCCTACGTTGTGGCCGGAGCACAGCAGCATGATGGTGAGGATCAGGCGGGTAATGCCGATCCACAGCAGGTACACCAGCAAGAAGGCCATGCCGAACTTGAGGCTGGCGATTACCGCCACGGt                         <  5:318399/149‑1 (MQ=255)
                                                                                                                               gagaATCAGCGGGTAGGCCGGGCCTACGTTGTGGCCGGAGCACAGCAGCATGATGGTGAGGATCAGGCGGGTAATGCCGATCCACAGCAGGTACACCAGCAAGAAGGCCATGCCGAACTTGAGGCTGGCGATTACCGCCACGGTCAGGc                    <  4:162937/149‑1 (MQ=255)
                                                                                                                                agaATCAGCGGGTAGGCCGGGCCTACGTTGTGGCCGGAGCACAGCAGCATGATGGTGAGGATCAGGCGGGTAATGCCGATCCACAGCAGGTACAcc                                                                        <  3:287604/96‑1 (MQ=255)
                                                                                                                                agaATCAGCGGGTAGGCCGGGCCTACGTTGTGGCCGGAGCACAGCAGCATGATGGTGAGGATCAGGCGGGTAATGCCGATCCACAGCAGGTACAcc                                                                        >  4:287604/1‑96 (MQ=255)
                                                                                                                                  aaTCAGCGGGTAGGGCGGGCCTACGTTGTGGCCGGAGCACAGCAGCATGATTGTGAGGATCAGGCGTGTAATGCCGATCCACAGCAGGTACACCAGCAAGAAGGCCATGCCGAACTTGAGGCTGGCGATTACCGCCACGGTCAGGccc                  <  3:247961/148‑1 (MQ=255)
                                                                                                                                  aaTCAGCGGGTAGGCCGGGCCTACGTTGTGGCCGGAGCACAGCAGCATGATGGTGAGGATCAGGCGGGTAATGCCGATCCACAGCAGGTACACCAGCAAGAAGGCCATGCCGAACTTGAGGCTGGCGATTACCGCCACGGTCAGGCCca                 <  1:46506/149‑1 (MQ=255)
                                                                                                                                     cacaacgTAGGCCCGGCCTACGTTGTGGCCGGAGCACAGCAGCATGATGGTGAGGATCAGGCGGGTAATGCCGATCCACAGCAGGTACACCAGCAAGAAGGCCATGCCGAACTTGAGGCTGGCGATTACCGCCACGGTCAGGCCca                 >  5:69063/7‑146 (MQ=255)
                                                                                                                                     cacaacgTAGGCCCGGCCTACGTTGTGGCCGGAGCACAGCAGCATGATGGTGAGGATCAGGCGGGTAATGCCGATCCACAGCAGGTACACCAGCAAGAAGGCCATGCCGAACTTGAGGCTGGCGATTACCGCCACGGTCAGGCCca                 <  6:69063/140‑1 (MQ=255)
                                                                                                                                      acaacgTAGGCCCGGCCTACGTTGTGGCCGGAGCACAGCAGCATGATGGTGAGGATCAGGCGGGTAATGCCGATCCACAGCAGGTACACCAGCAAGAAGGCCATGCCGaa                                                    >  8:141618/6‑110 (MQ=255)
                                                                                                                                      acaacgTAGGCCCGGCCTACGTTGTGGCCGGAGCACAGCAGCATGATGGTGAGGATCAGGCGGGTAATGCCGATCCACAGCAGGTACACCAGCAAGAAGGCCATGCCGaa                                                    <  7:141618/105‑1 (MQ=255)
                                                                                                                                        cGGGTAGGCCGGGCCTACGTTGTGGCCGGAGCACAGCAGCATGATGGTGAGGATCAGGCGGGTAATGCCGATCCACAGCAGGTACACCAGCAAGAAGGCCATGCCGAACTTGAGGCTGGCGATTACCGCCACGGTCAGGCCCAGCAGGc           <  1:56645/149‑1 (MQ=255)
                                                                                                                                           gTAGGCCGGGCCTACGTTGTGGCCGGAGCACAGCAGCATGATGGTGAGGATCAGGCGGGTAATGCCGATCCACAGCAGGTACACCAGCAAGAAGGCCATGCCGAACTTGAGGCTGGCGATTACCGCCACGGTCAGGCCCAGCAGGCTgg        <  2:350877/149‑1 (MQ=255)
                                                                                                                                            tAGGCCGGGCCTACGTTGTGGCCGGAGCACAGCAGCATGATGGTGAGGATCAGGCGGGTAATGCCGATCCACAGCAGGTACACCAGCAAGAAGGCCATGCCGAACTTGAGGCTGGCGATTACCGCCACGGTCAGGCCCAGCAGGCTGGt       <  3:3800/149‑1 (MQ=255)
                                                                                                                                               gCCGGGCCTACGTTGTGGCCGGAGCACAGCAGCATGATGGTGAGGATCAGGCGGGTAATGCCGATCCACAGCAGGTACACCAGCAAGAAGGCCATGCCGAACTTGAGGCTGGCGATTACCGCCACg                           <  6:236745/126‑1 (MQ=255)
                                                                                                                                               gCCGGGCCTACGTTGTGGCCGGAGCACAGCAGCATGATGGTGAGGATCAGGCGGGTAATGCCGATCCACAGCAGGTACACCAGCAAGAAGGCCATGCCGAACTTGAGGCTGGCGATTACCGCCACg                           >  5:236745/1‑126 (MQ=255)
                                                                                                                                                 cGGGCCTACGTTGTGGCCGGAGCACAGCAGCATGATGGTGAGGATCAGGCGGGTAATGCCGATCCACAGCAGGTACACCAGCAAGAAGGCCATGCCGAACTTGAGGCTGGCGATTACCGCCACGGTCAGGCCCAGCAGGCTGGTCcaca  <  5:36933/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCAGCGGGTAGGCCGGGCCTACGTTGTGGCCGGAGCACAGCAGCATGATGGTGAGGATCAGGCGGGTAATGCCGATCCACAGCAGGTACACCAGCAAGAAGGCCATGCCGAACTTGAGGCTGGCGATTACCGCCACGGTCAGGCCCAGCAGGCTGGTCCACA  >  NC_002947/1471407‑1471700

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: