Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,012,339 T→C 7.3% L561L (CTA→CTG PP_1793 ← group 2 family glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,012,3390TC 7.3% 239.5 / 6.9 82L561L (CTA→CTGPP_1793group 2 family glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base T (38/38);  new base C (3/3);  total (41/41)
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.53e-01

GCTATCCCGGTGCCAGACAAGAGAGAGGGTAGAATCAGCCCCTGTGCCGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAGTGTTTTTTCCGCCCCGAAAGAAATCCCGCCTTGAT  >  NC_002947/2012195‑2012484
                                                                                                                                                 |                                                                                                                                                   
gCTATCCCGGTGCCAGACAAGAGCGAGGGTAGAATCAGCCCCTGTGCCGGCATGAATGTGTAGCGTATTACCTTCCTCCATCTGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAgc                                                                                                                                                  <  7:193549/148‑1 (MQ=255)
  tATCCCGGTGCCAGACAAGAGAGAGGGTAGAATCAGCCCCTGTGCCGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGACGCGCTACAGc‑‑gcg                                                                                                                                             >  5:280720/1‑149 (MQ=255)
  tATCCCGGTGCCAGACAAGAGAGAGGGTAGAATCAGCCCCTGTGCCGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGACGCGCTACAGc‑‑gcg                                                                                                                                             >  5:349012/1‑149 (MQ=255)
   aTCCCGGTGCCAGACAAGAGAGAGGGTAGAATCAGCCCCTGTGCCGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGt                                                                                                                                            >  4:335838/1‑149 (MQ=255)
    tCCCGGTGCCAGACAAGAGAGAGGGTAGAATCAGCCCCTGTGCCGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGACGCGCTACAGc‑‑gcgcc                                                                                                                                           <  6:280720/149‑3 (MQ=255)
    tCCCGGTGCCAGACAAGAGAGAGGGTAGAATCAGCCCCTGTGCCGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGACGCGCTACAGc‑‑gcgcc                                                                                                                                           <  6:349012/149‑3 (MQ=255)
      ccGGTGCCAGACAAGAGAGAGGGTAGAATCAGCCCCTGTGCCGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTccc                                                                                                                                         >  5:365166/1‑149 (MQ=255)
          tGCCAGACAAGAGAGAGGGTAGAATCAGCCCCTGTGCCGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGc                                                                                                                                     <  3:79149/149‑1 (MQ=255)
             cAGACAAGAGAGAGGGTAGAATCAGCCAATGTGCCGTAATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGAGAAGCGGGACAGTCGGGCGCGCTGTAGA‑‑GCGTACCGGGCAGc                                                                                                                                  >  6:12317/1‑149 (MQ=255)
                   agagagagGGTAGAATCAGCCCCTGTGCCGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGAt                                                                                                                            <  7:203524/149‑1 (MQ=255)
                    gagagagGGTAGAATCAGCCCCTGTGCCGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATc                                                                                                                           >  8:145144/1‑149 (MQ=255)
                    gagagagGGTAGAATCAGCCCCTGTGCCGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATc                                                                                                                           <  4:247932/149‑1 (MQ=255)
                      gagagGGTAGAATCAGCCCCTGTGCCGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCg                                                                                                                          <  6:342125/148‑1 (MQ=255)
                         agGGTAGAATCAGCCCCTGTGCCGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTc                                                                                                                      >  6:386926/1‑149 (MQ=255)
                           ggTAGAATCAGCCCCTGTGCCGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTccc                                                                                                                    >  6:366274/1‑149 (MQ=255)
                                 aTCAGCCCCTGTGCCGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGACGCGCTACAGC‑‑GCGccc                                                                                                                                          <  8:282389/122‑1 (MQ=255)
                                 aTCAGCCCCTGTGCCGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGACGCGCTACAGC‑‑GCGccc                                                                                                                                          >  7:282389/1‑122 (MQ=255)
                                  tCAGCCCATGTGCCGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGACAGTCGCGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGc                                                                                                              >  2:225326/1‑148 (MQ=255)
                                  tCAGCCACTGTTCCGTAATGTATGTGTAGCGTTTTACCTGCCTCCATGGGGGATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAATGGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTGCCGGGCAGCTTTGATCGGGTCCCGCAGGCa                                                                                                             <  8:23946/149‑1 (MQ=255)
                                            tgCCGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCAc                                                                                                     <  3:190484/147‑1 (MQ=255)
                                             gCCGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCggg                                                                                                                        <  4:360244/127‑1 (MQ=255)
                                             gCCGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCggg                                                                                                                        >  3:360244/1‑127 (MQ=255)
                                                     gAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGc                                                                                                                                  >  6:365173/1‑109 (MQ=255)
                                                     gAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGc                                                                                                                                  <  5:365173/109‑1 (MQ=255)
                                                      aaTGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGcc                                                                                          >  5:234506/1‑148 (MQ=255)
                                                        tgtgTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCg                                                                                         >  1:240639/1‑147 (MQ=255)
                                                            tAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAgcgc                                                                                   >  6:46403/1‑149 (MQ=255)
                                                               cGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTg                                                                                <  6:217640/149‑1 (MQ=255)
                                                                gTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACtata                                                                                                 >  7:217928/1‑131 (MQ=255)
                                                                gTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACtata                                                                                                 <  8:217928/131‑1 (MQ=255)
                                                                   ttACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGa                                                                              >  6:191171/1‑147 (MQ=255)
                                                                      ccTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGTGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCgg                                                                                                                         >  3:45340/1‑101 (MQ=255)
                                                                      ccTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGTGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCgg                                                                                                                         <  4:45340/101‑1 (MQ=255)
                                                                        tGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAgcggc                                                                       >  8:30485/1‑149 (MQ=255)
                                                                           cTCCATCGGCCATTGCGAGGCTTTTTGTGCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAgcggcggc                                                                    <  6:217868/149‑1 (MQ=255)
                                                                              cATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGgcgc                                                                  >  6:12944/1‑148 (MQ=255)
                                                                               aTCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCaa                                                                <  1:151680/149‑1 (MQ=255)
                                                                                 cGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTg                                                                                            >  5:65674/1‑119 (MQ=255)
                                                                                 cGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTg                                                                                            <  6:65674/119‑1 (MQ=255)
                                                                                  ggCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGa                                                              <  7:145144/148‑1 (MQ=255)
                                                                                   gCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGCTa                                                            <  2:240639/149‑1 (MQ=255)
                                                                                     cATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAgcggcggc                                                                    >  3:314059/1‑139 (MQ=255)
                                                                                     cATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAgcggcggc                                                                    <  8:143270/139‑1 (MQ=255)
                                                                                     cATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAgcggcggc                                                                    <  4:314059/139‑1 (MQ=255)
                                                                                     cATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAgcggcggc                                                                    >  7:143270/1‑139 (MQ=255)
                                                                                         gCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCgtt                                                       <  5:366274/148‑1 (MQ=255)
                                                                                         gCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACGATAGCGTGCCGAAGCGcca                                                                                 >  2:360401/1‑121 (MQ=255)
                                                                                           gAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCTGGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCgttg                                                      >  1:294645/1‑149 (MQ=255)
                                                                                           gAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCgttgtt                                                    >  8:311079/1‑149 (MQ=255)
                                                                                            aGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGttgttg                                                   >  5:114484/1‑149 (MQ=255)
                                                                                             ggCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTtgttgt                                                  <  1:275218/149‑1 (MQ=255)
                                                                                               cTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTtgtgt                                                <  5:386926/149‑1 (MQ=255)
                                                                                               cTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTtgtgt                                                <  6:234506/149‑1 (MQ=255)
                                                                                                  tttGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTtgtgtg                                               <  7:30485/147‑1 (MQ=255)
                                                                                                     gTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAgc                                                                          >  6:353626/1‑117 (MQ=255)
                                                                                                     gTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAgc                                                                          <  5:353626/117‑1 (MQ=255)
                                                                                                      ttCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGAcgc                                         >  6:241579/1‑149 (MQ=255)
                                                                                                         aTCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGCGGCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGAcgcg                                        >  1:73332/1‑148 (MQ=255)
                                                                                                         aTCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGAcgcg                                        >  3:143970/1‑147 (MQ=255)
                                                                                                          tCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAg                                     <  5:241579/149‑1 (MQ=255)
                                                                                                            gCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAgtg                                   <  5:191171/149‑1 (MQ=255)
                                                                                                             ctctGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAgtgt                                  >  3:178436/1‑149 (MQ=255)
                                                                                                                tgttgtGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAGTGtt                                 <  3:335838/147‑1 (MQ=255)
                                                                                                                  ttgtGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAGTGtttttt                             >  7:293455/1‑149 (MQ=255)
                                                                                                                    gtgCGTAGCTGGCCAGTCGGGCGGGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCGCGCAGGCAAAACGCACTATATCGTGCCGAAGCGCCTTGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAGTGTTTTTTc                            <  1:225326/148‑1 (MQ=255)
                                                                                                                     tgCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAGTGTTTTTTc                            <  3:339318/147‑1 (MQ=255)
                                                                                                                     tgCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAGTGTTTTTTc                            >  4:339318/1‑147 (MQ=255)
                                                                                                                       cGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAGTGTTTTTTCCGcc                        <  4:143970/149‑1 (MQ=255)
                                                                                                                        gAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAGTGTTTTTTCCGccc                       >  5:361136/1‑149 (MQ=255)
                                                                                                                          aGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGCCGCGCAGTGTTTTTTCCGcccc                      >  6:87844/1‑148 (MQ=255)
                                                                                                                              ggCCAGTCGGTCGCGCTGTAGC‑‑GCGTCCCGGGCATCTTTGATCTGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGAAGCGCAGTGTTTTTTCCGCCCCgaaag                 <  5:12317/149‑1 (MQ=255)
                                                                                                                              ggCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAGTGTTTTTTCCGCCCCga                    >  5:158413/1‑146 (MQ=255)
                                                                                                                              ggCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAGTGTTTTTTCCGCCCCga                    >  2:127997/1‑146 (MQ=255)
                                                                                                                                   gTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAGTGTTTTTTCCGCCCCGAAAGAAATc            <  6:365166/149‑1 (MQ=255)
                                                                                                                                   gTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAGTGTTTTTTCCGCCCCGAAAGAAATc            <  6:158413/149‑1 (MQ=255)
                                                                                                                                    tCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAGTGTTTTTTCCGCCCCGAAAGAAATcc           >  4:113098/1‑149 (MQ=255)
                                                                                                                                      gggCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAGTGTTTTTTCCGCCCCGAAAGAAATCCCg         >  6:239732/1‑149 (MQ=255)
                                                                                                                                      gggCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAGTGTTTTTTCCGCCCCGAAAGAAATCCCg         <  6:114484/149‑1 (MQ=255)
                                                                                                                                        gcgcgcTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAGTGTTTTTTc                            >  7:325193/1‑128 (MQ=255)
                                                                                                                                        gcgcgcTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAGTGTTTTTTc                            <  8:325193/128‑1 (MQ=255)
                                                                                                                                           gcgcTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAGTGTTTTTTCCGCCCCGAAAGAAATCCCGCCtt     <  1:127997/149‑1 (MQ=255)
                                                                                                                                              cTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAGTGTTTTTTCCGCCCCGAAAGAAATCCCGCCTTGAt  <  5:46403/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GCTATCCCGGTGCCAGACAAGAGAGAGGGTAGAATCAGCCCCTGTGCCGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAGTGTTTTTTCCGCCCCGAAAGAAATCCCGCCTTGAT  >  NC_002947/2012195‑2012484

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: