Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,552,399:1 +G 8.6% coding (56/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,3991.G 8.6% 265.9 / 16.2 70coding (56/978 nt)PP_3137group 2 family glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base . (32/32);  new base G (4/2);  total (36/34)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GTAGATGTGCCTGACAATTCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACATGAGTATCCCGTTCACCGCCGCACGCGCC  >  NC_002947/3552252‑3552547
                                                                                                                                                       |                                                                                                                                                     
gTAGATGTGCCTGACAATTCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGa                                                                                                                                                        <  5:416660/148‑1 (MQ=255)
gTAGATGTGCCTGACAATTCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAg                                                                                                                                                      >  2:131325/1‑149 (MQ=255)
gTAGATGTGCCTGACAATTCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAg                                                                                                                                                      >  8:346099/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)
               aaTTCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCg                                                                                                                                      >  6:125425/1‑149 (MQ=255)
                tttCACTGGCAGTGGGGGTGCCTGCCATCATCAAGGCCCGTCATCCGGCCGCCGTTGCATCGGTCGAGTTAGTCCTATGGCAACTTTGATCGTGGTGGTTGTAATC‑‑GGACGCAATGAAGGCCAGCGCCTCGATCGCTGCCTTGTTTCGc                                                                                                                                     <  2:208279/148‑1 (MQ=255)
                 ttCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCt                                                                                                                                    >  2:171136/1‑149 (MQ=255)
                 ttCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCt                                                                                                                                    >  3:164507/1‑149 (MQ=255)
                 ttCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCACt                                                                                                                                    >  5:320899/1‑149 (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                                                                                                                                           >  8:229359/1‑124 (MQ=255)
                             gggggTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCAGCGCTCGAGGCGCTggcctt                                                                                                                                           <  7:229359/129‑6 (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                                                                                                          <  4:164507/149‑1 (MQ=255)
                                           tcatcaAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGt                                                                                                          <  7:346099/149‑1 (MQ=255)
                                             atcaAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTAc                                                                                                        >  3:15058/1‑149 (MQ=255)
                                             atcaAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTAc                                                                                                        >  6:299584/1‑149 (MQ=255)
                                                   gCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGa                                                                                                                     >  6:156302/1‑130 (MQ=255)
                                                   gCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGa                                                                                                                     >  6:149882/1‑130 (MQ=255)
                                                   gCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGa                                                                                                                     <  5:149882/130‑1 (MQ=255)
                                                   gCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGa                                                                                                                     <  5:156302/130‑1 (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                                                                                                  >  4:271754/1‑87 (MQ=255)
                                                                                                                 gtggtAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGAc                                                                                                  <  3:271754/87‑1 (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)
                                                                                                                                                  cTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACATGAGTATCCCGTTCACCGCCGCAcg      >  6:111037/1‑149 (MQ=255)
                                                                                                                                                   tCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCGGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACATGAGTATCCCGTTCACCGCCGCAcgc     >  8:385015/1‑149 (MQ=255)
                                                                                                                                                      aGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACATGAGTATCCCGTTCACCGCCGCACGCGcc  >  6:122472/1‑149 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
GTAGATGTGCCTGACAATTCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATC‑‑GGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACATGAGTATCCCGTTCACCGCCGCACGCGCC  >  NC_002947/3552252‑3552547

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: