Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,012,338 G→A 7.3% Q562* (CAG→TAG)  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,3380GA 7.3% 240.6 / 6.2 82Q562* (CAG→TAG) PP_1793group 2 family glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base G (38/38);  new base A (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 = 2.46e-01

GATGGCTATCCCGGTGCCAGACAAGAGAGAGGGTAGAATCAGCCCCTGTGCCGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAGTGTTTTTTCCGCCCCGAAAGAAATCCCGCCTTGAT  >  NC_002947/2012191‑2012484
                                                                                                                                                    |                                                                                                                                                    
gATGGCTATCCCGGTGCCAGACAAGAGAGAGGGTAGAATCAGCCCCTGTGCCGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTg                                                                                                                                                      >  5:217640/1‑148 (MQ=255)
    gCTATCCCGGTGCCAGACAAGAGCGAGGGTAGAATCAGCCCCTGTGCCGGCATGAATGTGTAGCGTATTACCTTCCTCCATCTGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAgc                                                                                                                                                  <  7:193549/148‑1 (MQ=255)
      tATCCCGGTGCCAGACAAGAGAGAGGGTAGAATCAGCCCCTGTGCCGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGACGCGCTACAGc‑‑gcg                                                                                                                                             >  5:349012/1‑149 (MQ=255)
      tATCCCGGTGCCAGACAAGAGAGAGGGTAGAATCAGCCCCTGTGCCGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGACGCGCTACAGc‑‑gcg                                                                                                                                             >  5:280720/1‑149 (MQ=255)
       aTCCCGGTGCCAGACAAGAGAGAGGGTAGAATCAGCCCCTGTGCCGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGt                                                                                                                                            >  4:335838/1‑149 (MQ=255)
        tCCCGGTGCCAGACAAGAGAGAGGGTAGAATCAGCCCCTGTGCCGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGACGCGCTACAGc‑‑gcgcc                                                                                                                                           <  6:349012/149‑3 (MQ=255)
        tCCCGGTGCCAGACAAGAGAGAGGGTAGAATCAGCCCCTGTGCCGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGACGCGCTACAGc‑‑gcgcc                                                                                                                                           <  6:280720/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                                                                                                                                          >  7:282389/1‑122 (MQ=255)
                                     aTCAGCCCCTGTGCCGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGACGCGCTACAGC‑‑GCGccc                                                                                                                                          <  8:282389/122‑1 (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                                                                                                                        >  3:360244/1‑127 (MQ=255)
                                                 gCCGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCggg                                                                                                                        <  4:360244/127‑1 (MQ=255)
                                                         gAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGc                                                                                                                                  <  5:365173/109‑1 (MQ=255)
                                                         gAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGc                                                                                                                                  >  6:365173/1‑109 (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                                                                                                 <  8:217928/131‑1 (MQ=255)
                                                                    gTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACtata                                                                                                 >  7:217928/1‑131 (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                                                                                            <  6:65674/119‑1 (MQ=255)
                                                                                     cGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTg                                                                                            >  5:65674/1‑119 (MQ=255)
                                                                                      ggCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGa                                                              <  7:145144/148‑1 (MQ=255)
                                                                                       gCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGCTa                                                            <  2:240639/149‑1 (MQ=255)
                                                                                         cATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAgcggcggc                                                                    <  4:314059/139‑1 (MQ=255)
                                                                                         cATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAgcggcggc                                                                    >  7:143270/1‑139 (MQ=255)
                                                                                         cATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAgcggcggc                                                                    >  3:314059/1‑139 (MQ=255)
                                                                                         cATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAgcggcggc                                                                    <  8:143270/139‑1 (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                                                <  6:234506/149‑1 (MQ=255)
                                                                                                   cTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTtgtgt                                                <  5:386926/149‑1 (MQ=255)
                                                                                                      tttGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTtgtgtg                                               <  7:30485/147‑1 (MQ=255)
                                                                                                         gTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAgc                                                                          <  5:353626/117‑1 (MQ=255)
                                                                                                         gTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAgc                                                                          >  6:353626/1‑117 (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                            >  4:339318/1‑147 (MQ=255)
                                                                                                                         tgCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAGTGTTTTTTc                            <  3:339318/147‑1 (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                    >  2:127997/1‑146 (MQ=255)
                                                                                                                                  ggCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAGTGTTTTTTCCGCCCCga                    >  5:158413/1‑146 (MQ=255)
                                                                                                                                       gTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAGTGTTTTTTCCGCCCCGAAAGAAATc            <  6:158413/149‑1 (MQ=255)
                                                                                                                                       gTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAGTGTTTTTTCCGCCCCGAAAGAAATc            <  6:365166/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)
                                                                                                                                                    |                                                                                                                                                    
GATGGCTATCCCGGTGCCAGACAAGAGAGAGGGTAGAATCAGCCCCTGTGCCGTCATGAATGTGTAGCGTATTACCTGCCTCCATCGGCCATTGCGAGGCTTTTTGTTCATCGCTCTGTTGTGCGAAGCGGGCCAGTCGGGCGCGCTGTAGC‑‑GCGTCCCGGGCAGCTTTGATCGGGTCCCGCAGGCAAAACGCACTATAGCGTGCCGAAGCGCCTGGATAGCGGCGGCGCAAGATAGCGTTGTTGTGTGCGACGCGCAGTGTTTTTTCCGCCCCGAAAGAAATCCCGCCTTGAT  >  NC_002947/2012191‑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: