Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,552,394 Δ1 bp 6.5% coding (51/978 nt) PP_3137 → group 2 family glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,552,3940C. 6.5% 247.6 / 10.3 77coding (51/978 nt)PP_3137group 2 family glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base C (36/36);  new base . (3/2);  total (39/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.10e-01

GGTAGATGTGCCTGACAATTCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACATGAGTATCCCGTTCACCGC  >  NC_002947/3552251‑3552537
                                                                                                                                                  |                                                                                                                                                 
ggTAGATGTGCCTGACAATTCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTc                                                                                                                                                <  4:268057/147‑1 (MQ=255)
ggTAGATGTGCCTGACAATTCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTc                                                                                                                                                >  3:268057/1‑147 (MQ=255)
 gTAGATGTGCCTGACAATTCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGa                                                                                                                                              <  5:416660/148‑1 (MQ=255)
 gTAGATGTGCCTGACAATTCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAg                                                                                                                                            >  8:346099/1‑149 (MQ=255)
 gTAGATGTGCCTGACAATTCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAg                                                                                                                                            >  2:131325/1‑149 (MQ=255)
      tgtgCCTGACAATTCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTg                                                                                                                                      >  5:418942/1‑149 (MQ=255)
        tgCCTGACAATTCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGcc                                                                                                                                    <  6:418942/149‑1 (MQ=255)
           cTGACAATTCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTg                                                                                                                                 >  1:205048/1‑149 (MQ=255)
             gACAATTCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGAGCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCt                                                                                                                                                 <  5:273562/133‑1 (MQ=255)
             gACAATTCACTGGCAGTGGGGGTGCCGGACATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCt                                                                                                                                                 >  6:273562/1‑133 (MQ=255)
                aaTTCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCg                                                                                                                            >  6:125425/1‑149 (MQ=255)
                 tttCACTGGCAGTGGGGGTGCCTGCCATCATCAAGGCCCGTCATCCGGCCGCCGTTGCATCGGTCGAGTTAGTCCTATGGCAACTTTGATCGTGGTGGTTGTAATC‑‑GGACGCAATGAAGGCCAGCGCCTCGATCGCTGCCTTGTTTCGc                                                                                                                           <  2:208279/148‑1 (MQ=255)
                  ttCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCt                                                                                                                          >  3:164507/1‑149 (MQ=255)
                  ttCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCt                                                                                                                          >  2:171136/1‑149 (MQ=255)
                  ttCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCACt                                                                                                                          >  5:320899/1‑149 (MQ=255)
                  ttCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCAGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTc                                                                                                                                                <  5:34721/129‑1 (MQ=255)
                  ttCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGACGAGTGAGGCCTATGGCAACGGTGATCAGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTc                                                                                                                                                >  6:34721/1‑129 (MQ=255)
                    cACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTg                                                                                                                         <  5:336993/148‑1 (MQ=255)
                    cACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATACGGCCTCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTg                                                                                                                         <  6:332275/148‑1 (MQ=255)
                     aCTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGgc                                                                                                                       >  7:126481/1‑149 (MQ=255)
                             tGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGGCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACgggg                                                                                                               <  8:265310/149‑1 (MQ=255)
                              gggggTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCAGCGCTCGAGGCGCTggcctt                                                                                                                                 <  7:229359/129‑6 (MQ=255)
                              gggggTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCAGCGCTCGAGGCGCTggcctt                                                                                                                                 >  8:229359/1‑124 (MQ=255)
                                gggTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCg                                                                                                            <  2:205048/149‑1 (MQ=255)
                                 ggTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGa                                                                                                           >  5:216038/1‑149 (MQ=255)
                                     ccGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACaaa                                                                                                       <  6:32285/149‑1 (MQ=255)
                                       ggCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGatccag                                                                                                      >  2:131770/1‑143 (MQ=255)
                                          catcatCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCAt                                                                                                  >  7:19845/1‑149 (MQ=255)
                                          catcatCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCAt                                                                                                  <  7:256275/149‑1 (MQ=255)
                                            tcatcaAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGt                                                                                                <  7:346099/149‑1 (MQ=255)
                                            tcatcaAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGt                                                                                                <  4:164507/149‑1 (MQ=255)
                                              atcaAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTAc                                                                                              >  6:299584/1‑149 (MQ=255)
                                              atcaAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTAc                                                                                              >  3:15058/1‑149 (MQ=255)
                                                    gCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGa                                                                                                           <  5:149882/130‑1 (MQ=255)
                                                    gCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGa                                                                                                           <  5:156302/130‑1 (MQ=255)
                                                    gCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGa                                                                                                           >  6:156302/1‑130 (MQ=255)
                                                    gCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGa                                                                                                           >  6:149882/1‑130 (MQ=255)
                                                      ccGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGTGGCCGACAAAGTCATGTACGTCGACTc                                                                                      <  6:19973/149‑1 (MQ=255)
                                                      ccGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTc                                                                                      >  8:335825/1‑149 (MQ=255)
                                                       cGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCg                                                                                     >  7:267672/1‑149 (MQ=255)
                                                       cGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCg                                                                                     >  8:89900/1‑149 (MQ=255)
                                                           aTCCGGCCGCCGTTTCATGGGTCGTTTGAGGCCTTTTGCAACGGTTTTTGGGGTGTAAGTACATCTAAATCATAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCgg                                                                                    <  1:131325/149‑1 (MQ=255)
                                                             ccGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTc                                                                                >  6:73503/1‑148 (MQ=255)
                                                               ggCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCAc                                                                             <  5:299584/149‑1 (MQ=255)
                                                               ggCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCAc                                                                             <  5:125425/149‑1 (MQ=255)
                                                                  cgccgTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGa                                                                          <  3:336892/149‑1 (MQ=255)
                                                                   gccgTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGAt                                                                         <  8:126481/149‑1 (MQ=255)
                                                                    ccgTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATg                                                                        <  5:73503/149‑1 (MQ=255)
                                                                                tCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGt                                                            >  1:254918/1‑149 (MQ=255)
                                                                                     tGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGgcc                                                       <  6:216038/149‑1 (MQ=255)
                                                                                     tGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCACTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGATCGCTGCAGTTGgc                                                        <  6:320899/148‑1 (MQ=255)
                                                                                      gAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGgccg                                                      >  1:125067/1‑149 (MQ=255)
                                                                                           cTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAAc                                                 <  2:125067/149‑1 (MQ=255)
                                                                                           cTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCAAGGGGCCGAAAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAAc                                                 >  6:86856/1‑149 (MQ=255)
                                                                                            tATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAAcc                                                >  8:15106/1‑149 (MQ=255)
                                                                                               ggCAACGGTAATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACATTg                                             >  5:36012/1‑149 (MQ=255)
                                                                                                gCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTgg                                            <  1:142914/149‑1 (MQ=255)
                                                                                                     ggTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACTTTGGCGTgg                                       <  2:336540/149‑1 (MQ=255)
                                                                                                     ggTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTgg                                       <  8:267672/149‑1 (MQ=255)
                                                                                                       tGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGGCGCTGCCTTGt                                                                                                                                >  1:91174/1‑59 (MQ=255)
                                                                                                                  gtggtAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGAc                                                                                        <  3:271754/87‑1 (MQ=255)
                                                                                                                  gtggtAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGAc                                                                                        >  4:271754/1‑87 (MQ=255)
                                                                                                                    ggtAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTg                         >  7:17182/1‑148 (MQ=255)
                                                                                                                     gtAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTgg                                            <  7:18823/128‑1 (MQ=255)
                                                                                                                     gtAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTgg                                            >  8:18823/1‑128 (MQ=255)
                                                                                                                       aaTC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACa                     >  8:329897/1‑149 (MQ=255)
                                                                                                                         aTC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCGGTACTGGCGCTGGACAt                    <  2:254918/149‑1 (MQ=255)
                                                                                                                         aTC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACAt                    <  7:329897/149‑1 (MQ=255)
                                                                                                                         aTC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACAt                    <  4:15058/149‑1 (MQ=255)
                                                                                                                              ggCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACATGa                  <  7:15106/148‑1 (MQ=255)
                                                                                                                              ggCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACATGAg                 >  5:13263/1‑149 (MQ=255)
                                                                                                                                         aaggcAGCGCTCGAGGCGCTGGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACATGAGTATCCCGtt        >  8:225278/5‑148 (MQ=255)
                                                                                                                                          ggCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTATGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACATGa                  <  1:398541/136‑1 (MQ=255)
                                                                                                                                          ggCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTATGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACATGa                  >  2:398541/1‑136 (MQ=255)
                                                                                                                                             cAGCGCCTCGAGCACTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGTCTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACATGAGTATCCCGTTCAccgc  >  4:175390/1‑149 (MQ=255)
                                                                                                                                             cAGCGCTCGAGGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCAc                                                                             <  2:196939/74‑1 (MQ=255)
                                                                                                                                             cAGCGCTCGAGGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCAc                                                                             >  1:196939/1‑74 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GGTAGATGTGCCTGACAATTCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACATGAGTATCCCGTTCACCGC  >  NC_002947/3552251‑3552537

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: