Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,556,415 C→G 6.1% P278A (CCC→GCC)  PP_3140 → group 2 family glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,556,4150CG 6.1% 212.8 / 3.0 65P278A (CCC→GCC) PP_3140group 2 family glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base C (31/30);  new base G (2/2);  total (33/32)
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.68e-01

GGCCTGGACCTTTGCCAGCAGGGCAAACCGCCTGCTTTCTGCCCCGAGGCACTGGTCACTAGCCAGTTCCCACTTAGCCAGCAGGGCCTCAACAGCCAACGTACCCGCTGGGAGCACGGCCATCTGGGCTTGATGCTTGCCGACGCACCCAAGCGGGCGATGGCTGCCATCACCCAGCGTAACGCTAGCCTGGCGGCCTTGACCCTGGACTTGCTGGTGCCGCCCCTGGCGTTGCTGGTACTGACCTTGCTCGGGCTTAACCTGGTGACCTGGCTGGCGTACCTGCTGTTCGGCC  >  NC_002947/3556268‑3556562
                                                                                                                                                   |                                                                                                                                                   
ggCCTGGACCTTTGCCAGCAGGGCAAACCGCCTGCTTTCTGCCCCGAGGCACTGGTCACTAGCCAGTTCCCACTTAGCCAGCAGGGCCTCAACAGCCAACGTACCCGCTGGGAGCACGGCCATCTGGGCTTGATGCTTGCCGACGCAcc                                                                                                                                                    >  2:139866/1‑149 (MQ=255)
       aCCTTTGCCAGCAGGGCAAACCGCCTGCTTTCTGCCCCGAGGCACTGGTCACTAGCCAGTTCCCACTTAGCCAGCAGGGCCTCAACAGCCAACGTACCCGCTGGGAGCACGGCCATCTGGGCTTGATGCTTGCCGACGCACCCAAGCgg                                                                                                                                             <  6:374790/149‑1 (MQ=255)
        ccTTTGCCAGCAGGGCAAACCGCCTGCTTTCTGCCCCGAGGCACTGGTCACTAGCCAGTTCCCACTTAGCCAGCAGGGCCTCAACAGCCAACGTACCCGCTGGGAGCACGGCCATCTGGGCTTGATGCTTGCCGACGCACCCAAGCggg                                                                                                                                            <  4:344444/149‑1 (MQ=255)
         cattGCCAGCAGGGCAAACCGCCTGCTTTCTGCCCCGAGGCACTGGTCACTAGCCAGTTCCCACTTAGCCAGCAGGGCCTCAACAGCCAACGTACCCGCTGGGAGCACGGCCATCTGGGCTTGATGCTTGCCGACGCACACAAGCGGGc                                                                                                                                           >  8:297879/3‑149 (MQ=255)
         cTTTGCCAGCAGGGCAAACCGCCTGCTTTCTGCCCCGAGGCACTGGTCACTAGCCAGTTCCCACTTAGCCAGCAGGGCCTCAACAGCCAACGTACCCGCTGGGAGCACGGCCATCTGGGCTTGATGCTTGCCGACGCACCCAAGCGGGc                                                                                                                                           >  4:129478/1‑149 (MQ=255)
                     ggCAAACCGCCTGCTTTCTGCCCCGAGGCACTGGTCACTAGCCAGTTCCCACTTAGCCAGCAGGGCCTCAACAGCCAACGTACCCGCTGGGAGCACGGCCATCTGGGCTTGATGCTTGCCGACGCACCCAAGCGGGCGATg                                                                                                                                       <  1:335462/141‑1 (MQ=255)
                     ggCAAACCGCCTGCTTTCTGCCCCGAGGCACTGGTCACTAGCCAGTTCCCACTTAGCCAGCAGGGCCTCAACAGCCAACGTACCCGCTGGGAGCACGGCCATCTGGGCTTGATGCTTGCCGACGCACCCAAGCGGGCGATg                                                                                                                                       >  2:335462/1‑141 (MQ=255)
                       cAAACCGCCTGCTTTCTGCCCCGAGGCACTGGTCACTAGCCAGTTCCCACTTAGCCAGCAGGGCCTCAACAGCCAACGTACCCGCTGGGAGCACGGCCATCTGGGCTTGATGCTTGCCGACGCACCCAAGCGGGCGATGGCTGCCATCa                                                                                                                             <  3:75917/149‑1 (MQ=255)
                       cAAACCGCCTGCTTTCTGCCCCGAGGCACTGGTCACTAGCCAGTTCCCACTTAGCCAGCAGGGCCTCAACAGCCAACGTACCCGCTGGGAGCACGGCCATCTGGGCTTGATGCTTGCCGACGCACCCAAGCGGGCGATGGCTGCCATCa                                                                                                                             <  1:139866/149‑1 (MQ=255)
                                tGCTTTCTGCCCCGAGGCACTGGTCACTAGCCAGTTCCCACTTAGCCAGCAGGGCCTCAACAGCCAACGTACCCGCTGGGAGCACGGCCATCTGGGCTTGATGCTTGCCGACGCACCCAAGCGGGCGATGGCTGCCATCACCCAGCGt                                                                                                                     >  4:77800/1‑148 (MQ=255)
                                   tttCTGCCCCGAGGCACTGGTCACTAGCCAGTTCCCACTTAGCCAGCAGGGCCTCAACAGCCAACGTACCCGCTGGGAGCACGGCCATCTGGGCTTGATGCTTGCCGACGCACCCAAGCGGGCGATGGCTGCCATCACCCAGCGTAACg                                                                                                                 >  2:227419/1‑149 (MQ=255)
                                         ccccGAGGCACTGGTCACTAGCCAGTTCCCACTTAGCCAGCAGGGCCTCAACAGCCAACGTACCCGCTGGGAGCACGGCCATCTGGGCTTGATGCTTGCCGACGCACCCAAGCGGGCGATGGCTGCCATCACCCAGCGTAACGCTAGcc                                                                                                           >  7:160524/1‑149 (MQ=255)
                                            cGAGGCACTGGTCACT‑‑CCAGTTCCCACTTAGCCAGCAGGGCCTCAACAGCCAACGTACCCGCTGGGAGCACGGCCATCTGGGCTTGATGCTTGCCGACGCACCCAAGCGGGCGATGGCTGCCATCACCCAGCGTAACGCTAGCCTggcg                                                                                                      <  3:319714/149‑1 (MQ=255)
                                                gCACTGGTCACTAGCCAGTTCCCACTTAGCCAGCAGGGCCTCAACAGCCAACGTACCCGCTGGGAGCACGGCCATCTGGGCTTGATGCTTGCCGACGCACCCAAGCGGGCGATGGCTGCCATCACCCAGCGTAACGCTAGCCTggcgg                                                                                                     <  8:261056/148‑1 (MQ=255)
                                                 cACTGCTCACTAGCCAGTTCCCACTTAGCCAGCAGGGCCTCAACAGCCAACGTACCCGCTGGGAGCACGGCCATCTGGGCTTGATGCTTGCCGACGCACCCAAGCGGGCGATGGCTGCCATCACCCAGCGTAACGCTAGCCTGGCGGcc                                                                                                   <  8:248369/149‑1 (MQ=255)
                                                  aCTGGTCACTAGCCAGTTCCCACTTAGCCAGCAGGGCCTCAACAGCCAACGTACCCGCTGGGAGCACGGCCATCTGGGCTTGATGCTTGCCGACGCACCCAAGCGGGCGATGGCTGCCATCACCCAGCGTAACGCTAGCCTGGCGGCCt                                                                                                  >  7:301084/1‑149 (MQ=255)
                                                      gTCACTAGCCAGTTCCCACTTAGCCAGCAGGGCCTCAACAGCCAACGTACCCGCTGGGAGCACGGCCATCTGGGCTTGATGCTTGCCGACGCACCCAAGCGGGCGATGGCTGCCATCACCCAGCGTAACGCTAGCCTGGCGGCCTTGAc                                                                                              >  5:203558/1‑149 (MQ=255)
                                                         aCTAGCCAGTTCCCACTTAGCCAGCAGGGCCTCAACAGCCAACGTACCCGCTGGGAGCACGGCCATCTGGGCTTGATGCTTGCCGACGCACCCAAGCGGGCGATGGCTGCCATCACCCAGCGTAACGCTAGCCTGGCGGCCTTGACCCt                                                                                           >  5:89803/1‑149 (MQ=255)
                                                             gCCAGTTCCCACTTAGCCAGCAGGGCCTCAACAGCCAACGTACCCGCTGGGAGCACGGCCATCTGGGCTTGATGCTTGCCGACGCACCCAAGCGGGCGATGGCTGCCATCACCCAGCGTAACGCTAGCCTGGCGGCCTTGACCCTGGAc                                                                                       >  8:282853/1‑149 (MQ=255)
                                                               cAGTTCCCACTTAGCCAGCAGGGCCTCAACAGCCAACGTACCCGCTGGGAGCACGGCCATCTGGGCTTGATGCTTGCCGACGCACCCAAGCGGGCGATGGCTGCCATCACCCAGCGTAACGCTAGCCTGGCGGCCTTGACCCTGGACtt                                                                                     >  6:9080/1‑149 (MQ=255)
                                                               cAGTTCCCACTTAGCCAGCAGGGCCTCAACAGACAACGTACCCGCTGGGAGCACGGCCATCTGGGCTTGATGCTTGCCGACGCACCCAAGCGGGCGATGGCTGCCATCACCCAGAGTAACGCTAGCCTGGCGGCCTTGATCCTGGACtt                                                                                     >  2:282214/1‑149 (MQ=255)
                                                                aGTTCCCACTTAGCCAGCAGGGCCTCAACAGCCAACGTACCCGCTGGGAGCACGGCCATCTGGGCTTGATGCTTGCCGACGCACCCAAGCGGGCGATGGCTGCCATCACCCAGCGTAACGCTAGCCTGGCGGCCTTGACCCTGGACTTg                                                                                    <  7:298906/149‑1 (MQ=255)
                                                                    cccACTTAGCCAGCAGGGCCTCAACAGCCAACGTACCCGCTGGGAGCACGGCCATCTGGGCTTGATGCTTGCCGACGCACCCAAGCGGGCGATGGCTGCCATCACCCAGCGTAACGCTAGCCTGGCGGCCTTGACCCT‑GACTTGCTGGt                                                                               >  8:245220/1‑149 (MQ=255)
                                                                     ccACTTAGCCAGCAGGGCCTCAACAGCCAACGTACCCGCTGGGAGCACGGCCATCTGGGCTTGATGCTTGCCGACGCACCCAAGCGGGCGATGGCTGCGATCACCCAGCGTAACGCTAGCCTGGCGGCCTTGACCCTGGACTTGCTgg                                                                                <  5:332815/148‑1 (MQ=255)
                                                                       aCTTAGCCAGCAGGGCCTCAACAGCCAACGTACCCGCTGGGAGCACGGCCATCTGGGCTTGATGCTTGCCGACGCACCCAAGCGGGCGATGGCTGCCATCACCCAGCGTAACGCTAGCCTGGCGGCCTTGACCCTGGACTTGCTGGTgc                                                                             <  3:208603/149‑1 (MQ=255)
                                                                       aCTTAGCCAGCAGGGCCTCAACAGCCAACGTACCCGCTGGGAGCACGGCCATCTGGGCTTGATGCTTGCCGACGCACCCAAGCGGGCGATGGCTGCCATCACCCAGCGTAACGCTAGCCTGGCGGCCTTGACCCTGGACTTGCTGGTgc                                                                             >  1:206256/1‑149 (MQ=255)
                                                                           aGCCAGCAGGGCCTCAACAGCCAACGTACCCGCTGGGAGCACGGCCATCTGGGCTTGATGCTTGCCGACGCACCCAAGCGGGCGATGGCTGCCATCACCCAGCGTAACGCTAGCCTGGCGGCCTTGACCCTGGACTTGCTGGTgccgc                                                                          <  8:301084/148‑1 (MQ=255)
                                                                               agcagGGCCTCAACAGCCAACGTACCCGCTGGGAGCACGGCCATCTGGGCTTGATGCTTGCCGACGCACCCAAGCGGGCGATGGCTGCCATCACCCAGCGTAACGCTGGCCTGGCGGCCTTGACCCTGGACTTGCTGGTGCCGCCCCTg                                                                     >  5:326479/1‑149 (MQ=255)
                                                                                      ccTCAACAGCCAACGTACCCGCTGGGAGCACGGCCATCTGGGCTTGATGCTTGCCGACGCACCCAAGCgg                                                                                                                                             >  8:231667/1‑70 (MQ=255)
                                                                                      ccTCAACAGCCAACGTACCCGCTGGGAGCACGGCCATCTGGGCTTGATGCTTGCCGACGCACCCAAGCgg                                                                                                                                             <  7:231667/70‑1 (MQ=255)
                                                                                      ccTCAACAGCCAACGTACCCGCTGGGAGCACGGCCATCTGGGCTTGATGCTTGCCGACGCACCCAAGCGGGCGATGGCTGCCATCACCCAGCGTAACGCTAGCCTGGCGGCCTTGACCCTGGACTTGCTGGTGCCGCCCCTGGCGTTGc                                                              <  5:9080/149‑1 (MQ=255)
                                                                                       cTCAACAGCCAACGTACCCGCTGGGAGAACGGCCATCTGGGCTTGATGCTTGCCGACGCACCCAAGCGGGCGATGGCTGACATCAACCAGCGTAACGATATCCTGGCGGCCTTGACCCTCGACTAGCTGGTGCCGCCCCTGGCGTTCCt                                                             >  8:85408/1‑149 (MQ=255)
                                                                                          aaCAGCCAACGTACCCGCTGGGAGCACGGCCATCTGGGCTTGTTGCTTGCCGACGCACCCAAGCGGGCGATGGCTGCCATCACCCAGCGTAACGCTAGCCTGGCGGCCTTGACCCT‑GACTTGCTGGTGCCGCCCCTGGCGTTGCTgg                                                           <  7:245220/147‑1 (MQ=255)
                                                                                          aaCAGCCAACGTACCCGCTGGGAGCACGGCCATCTGGGCTTGATGCTTGCCGACGCACCCAAGCGGGCGATGGCTGCCATCACCCAGCGTAACGCTAGCCTGGCGGCCTTGACCCTGGACTTGCTGGTGCCGCCCCTGGCGTTGCTg                                                            <  7:282853/147‑1 (MQ=255)
                                                                                           aCAGCCAACGTACCCGCTGGGAGCACGGCCATCTGGGCTTGATGCTTGCCGACGCACCCAAGCGGGCGATGGCTGCCATCACCCAGCGTAACGCTAGCCTGGCGGCCTTGACCCTGGACTTGCTGGTGCCGCCCCTGGCGTTGCTGGTa                                                         <  3:129478/149‑1 (MQ=255)
                                                                                           aCAGCCAACGTACCCGCTGGGAGCACGGCCATCTGGGCTTGATGCTTGCCGACGCACCCAAGCGGGCGATGGCTGCCATCACCCAGCGTAACGCTAGCCTGGCGGCCTTGACCCTGGACTTGCTGGTGCCGCCCCTGGCGTTGCTGGTa                                                         <  2:206256/149‑1 (MQ=255)
                                                                                               ccAACGTACCCGCTGGGAGCACGGCCATCTGGGCTTGATGCTTGCCGACGCACCCAAGCGGGCGATGGCTGCCATCACCCAGCGTAACGCTAGCCTGGCGGCCTTGACCCTGGACTTGCTGGTGCCGCCCCTGGCGTTGCTGGTACTGa                                                     <  6:203558/149‑1 (MQ=255)
                                                                                                cAACGTACCCGCTGGGAGCACGGCCATCTGGGCTTGATGCTTGCCGACGCACCCAAGCGGGCGATGGCTGCCATCACCCAGCGTAACGCTAGCCTGGCGGCCTTGACCCTGGACTTGCTGGTGCCGCCCCTGGCGTTGCTGGTAc                                                        <  6:187526/145‑1 (MQ=255)
                                                                                                cAACGTACCCGCTGGGAGCACGGCCATCTGGGCTTGATGCTTGCCGACGCACCCAAGCGGGCGATGGCTGCCATCACCCAGCGTAACGCTAGCCTGGCGGCCTTGACCCTGGACTTGCTGGTGCCGCCCCTGGCGTTGCTGGTAc                                                        >  5:187526/1‑145 (MQ=255)
                                                                                                       cccGCTGGGAGCACGGCCATCTGGGCTTGATGCTTGCCGACGCACCCAAGCGGGCGATGGCTGCCATCACCCAGCGTAACGCTAGCCTGGCGGCCTTGACCCTGGACTTGCTGGTGCCGCCCCTGGCGTTGCTGGTACTGACCTTGCTc                                             >  3:255138/1‑149 (MQ=255)
                                                                                                          gCTGGGAGCACGGCCATCTGGGCTTGATGCTTGCCGACGCACCCAAGCGGGCGATGGCTGCCATCACCCAGCGTAACGCTAGCCTGGCGGCCTTGATCCTGGACTTGCTGGTGCCGCCCCTGGCGTTGCTGGTACTGACCTTGCTCggg                                          <  1:282214/149‑1 (MQ=255)
                                                                                                              ggAGCACGGCCATCTGGGCTTGATGCTTGCCGACGCACCCAAGCGGGCGATGGCTGCCATCACCCAGCGTAACGCTAGCCTGGCGGCCTTGACCCTGGACTTGCTGGTGCCGCCCCTGGCGTTGCTGGTACTGACCTTGCTCGGGCtt                                       <  4:255138/148‑1 (MQ=255)
                                                                                                                aGCACGGCCATCTGGGCTTGATGCTTGCCGACGCACCCAAGCGGGCGATGGCTGCCATCACCCAGCGTAACGCTGGCCTGGCGGCCTTGACCCTGGACTTGCTGGTGCCGCCCCTGGCGTTGCTGGTACTGACCTTGCTCGGGCTTAAc                                    <  6:326479/149‑1 (MQ=255)
                                                                                                                 gCACGGCCATCTGGGCTTGATGCTTGCCGACGCACCCAAGCGGGCGATGGCTGCCATCACCCAGCGTAACGCTAGCCTGGCGGCCTTGACCCTGGACTTGCTGGTGCCGCCCCTGGCGTTGCTGGTACTGACCTTGCTCGGGCTTAAcc                                   <  3:229449/149‑1 (MQ=255)
                                                                                                                 gCACGGCCATCTGGGCTTGATGCTTGCCGACGCACCCAAGCGGGCGATGGCTGCCATCACCCAGCGTAACGCTAGCCTGGCGGCCTTGACCCTGGACTTGCTGGTGCCGCCCCTGGCGTTGCTGGTACTGACCTTGCTCGGGCTTAAcc                                   >  4:203392/1‑149 (MQ=255)
                                                                                                                   aCGGCCATCTGGGCTTCATGCTTGCCCACGCACCCAAGCCGGCGATGGCTGCCATCACCCAGCGTAACGCTAGCCAGGCGGCCTTGACCCTGGACTTGCTGGTCCCGCCCCTGGCCTTGCTGGTACTGACATTGCTCGGGCTTAACCTg                                 >  6:246461/1‑149 (MQ=255)
                                                                                                                     ggCCATCTGGGCTTGATGCTTGCCGACGCACCCAAGCGGGCGATGGCTGCCATCACCCAGCGTAACGCTAGCCTGGCGGCCTTGACCCTGGACTTGCTGGTGCCGCCCCTGGCGTTGc                                                              <  5:67019/118‑1 (MQ=255)
                                                                                                                     ggCCATCTGGGCTTGATGCTTGCCGACGCACCCAAGCGGGCGATGGCTGCCATCACCCAGCGTAACGCTAGCCTGGCGGCCTTGACCCTGGACTTGCTGGTGCCGCCCCTGGCGTTGc                                                              >  6:67019/1‑118 (MQ=255)
                                                                                                                     ggCCATCTGGGCTTGATGCTTGCCGACGCACCCAAGCGGGCGATGGCTGCCATCACCCAGCGTAACGCTAGCCTGGCGGCCTTGACCCTGGACTTGCTGGTGCCGCCCCTGGCGTTGCTGGTACTGACCTTGCTCGGGCTTAACCTGGt                               <  7:103529/149‑1 (MQ=255)
                                                                                                                     ggCCATCTGGGCTTGATGCTTGCCGACGCACCCAAGCGGGCGATGGCTGCCATCACCCAGCGTAACGCTAGCCTGGCGGCCTTGACCCTGGACTTGCTGGTGCCGCCCCTGGCGTTGCTGGTACTGACCTTGCTCGGGCTTAACCTGGt                               >  8:103529/1‑149 (MQ=255)
                                                                                                                        cATCTGGGCTTGATGCTTGCCGACGCACCCAAGCGGGCGATGGCTGCCATCACCCAGCGTAACGCTAGCCTGGCGGCCTTGACCCTGGACTTGCTGGTGCCGCCCCTGGCGTTGCTgg                                                           >  5:249988/1‑118 (MQ=255)
                                                                                                                        cATCTGGGCTTGATGCTTGCCGACGCACCCAAGCGGGCGATGGCTGCCATCACCCAGCGTAACGCTAGCCTGGCGGCCTTGACCCTGGACTTGCTGGTGCCGCCCCTGGCGTTGCTgg                                                           <  6:249988/118‑1 (MQ=255)
                                                                                                                                cTTGATGCTTGCCGACGCACCCAAGCGGGCGATGGCTGCCATCACCCAGCGTAACGCTAGCCTGGCGGCCTTGACCCTGGACTTGCTGGTGCCGCCCCTGGCGTTGCTGGTACTGACCTTGCTCGGGCTTAACCTGGTGACctggctgg                    >  3:346108/1‑149 (MQ=255)
                                                                                                                                 ttGATGCTTGCCGACGCACCCAAGCGGGCGATGGCTGCCATCACCCAGCGTAACGCTAGCCTGGCGGCCTTGACCCTGGACTTGCTGGTGCCGCCCCTGGCGTTGCTGGTACTGACCTTGCTCGGGCtt                                       >  5:239934/1‑129 (MQ=255)
                                                                                                                                 ttGATGCTTGCCGACGCACCCAAGCGGGCGATGGCTGCCATCACCCAGCGTAACGCTAGCCTGGCGGCCTTGACCCTGGACTTGCTGGTGCCGCCCCTGGCGTTGCTGGTACTGACCTTGCTCGGGCtt                                       <  6:239934/129‑1 (MQ=255)
                                                                                                                                 ttGATGCTTGCCGACGCACCCAAGCGGGCGATGGCTGCCATCACCCAGCGTAACGCTAGCCTGGCGGCCTTGACCCTGGACTTGCTGGTGCCGCCCCTGGCGTTGCTGGTACTGACCTTGCTCGGGCTTAACCTGGTGACCtggctggc                   >  6:373297/1‑149 (MQ=255)
                                                                                                                                 ttGATGCTTGCCGACGCACCCAAGCGGGCGATGGCTGCCATCACCCAGCGTAACGCTAGCCTGGCGGCCTTGACCCTGGACTTGCTGGTGCCGCCCCTGGCGTTGCTGGTACTGACCTTGCTCGGGCTTAACCTGGTGACCtggctggc                   >  8:116934/1‑149 (MQ=255)
                                                                                                                                   gATGCTTGCCGACGCACCCAAGCGGGCGATGGCTGCCATCACCCAGCGTAACGCTAGCCTGGCGGCCTTGACCCTGGACTTGCTGGTGCCGCCCCTGGCGTTGCTGGTACTGACCtt                                                 <  7:136260/117‑1 (MQ=255)
                                                                                                                                   gATGCTTGCCGACGCACCCAAGCGGGCGATGGCTGCCATCACCCAGCGTAACGCTAGCCTGGCGGCCTTGACCCTGGACTTGCTGGTGCCGCCCCTGGCGTTGCTGGTACTGACCtt                                                 >  8:136260/1‑117 (MQ=255)
                                                                                                                                      gCTTGCCGACGCACCCAAGCGGGCGATGGCTGCCATCACCCAGCGTAACGCTAGCCTGGCGGCCTTGACCCTGGACTTGCTGGTGCCGCCCCTGGCGTTGCTGGTACTGACCTTGCTCGGGCTTAACCTGGTGACCTGGCTGGCGTAcc              >  3:113429/1‑149 (MQ=255)
                                                                                                                                        ttGCCGACGCAGCCAAGCGGGCGATGGCTGCCATCACCCAGCGTAACGCTAGCCTGGCGGCCTTGACCCTGGACTTGCTGGTGCCGCCCCTGGCGTTGCTGGTACTGACCTTGCTCGGGCTTAACCTGGTGAc                            <  5:264003/133‑1 (MQ=255)
                                                                                                                                        ttGCCGACGCAGCCAAGCGGGCGATGGCTGCCATCACCCAGCGTAACGCTAGCCTGGCGGCCTTGACCCTGGACTTGCTGGTGCCGCCCCTGGCGTTGCTGGTACTGACCTTGCTCGGGCTTAACCTGGTGAc                            >  6:264003/1‑133 (MQ=255)
                                                                                                                                        ttGCCGACGCAGCCAAGCGGGCGATGGCTGCCATCACCCAGCGTAACGCTAGCCTGGCGGACTTGACCCTTGACTTGCTGGTGCCGCCCCTGGCGTTGCTGGTACTGACCTTGCTCGGGCTTAACCTGGTGAc                            <  5:264009/133‑1 (MQ=255)
                                                                                                                                        ttGCCGACGCAGCCAAGCGGGCGATGGCTGCCATCACCAAGCGTAACGCTAGCCTGGCGGCCTTGACCCTGGACTTGCTGGTGCCGCCCCTGGCGTTGCTGGTACTGAACTTGCTCTGGCTTAAACTGGTGAc                            >  6:264009/1‑133 (MQ=255)
                                                                                                                                           ccGACGCACCCAAGCGGGCGATGGCTGCCATCACCCAGCGTAACGCTAGCCTGGCGGCCTTGACCCTGGACTTGCTGGTGCCGCCCCTGGCGTTGCTGGTACTGACCTTGCTCGGGCTTAACCTGGTGACCTGGCTGGCGTACctgctg         <  7:116934/149‑1 (MQ=255)
                                                                                                                                                  aCCCAAGCGGGCGATGGCTGCCATCACCCAGCGTAACGCTAGCCTGGCGGCCTTGACCCTGGACTTGCTGGTGCCGCCCCTGGCGTTGCTGGTACTGACCTTGCTCGGGCTTAACCTGGTGACCTGGCTGGCGTACCTGCTGTTCGGcc  >  1:198627/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GGCCTGGACCTTTGCCAGCAGGGCAAACCGCCTGCTTTCTGCCCCGAGGCACTGGTCACTAGCCAGTTCCCACTTAGCCAGCAGGGCCTCAACAGCCAACGTACCCGCTGGGAGCACGGCCATCTGGGCTTGATGCTTGCCGACGCACCCAAGCGGGCGATGGCTGCCATCACCCAGCGTAACGCTAGCCTGGCGGCCTTGACCCTGGACTTGCTGGTGCCGCCCCTGGCGTTGCTGGTACTGACCTTGCTCGGGCTTAACCTGGTGACCTGGCTGGCGTACCTGCTGTTCGGCC  >  NC_002947/3556268‑3556562

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: ‑