Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,552,399:1 +G 9.0% 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 9.0% 248.3 / 17.7 67coding (56/978 nt)PP_3137group 2 family glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base . (34/27);  new base G (3/3);  total (37/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.95e-01

GTAGATGTGCCTGACAATTCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACATGAGTATCCCGTTCACCGCCGCACGCGCC  >  NC_002947/3552252‑3552547
                                                                                                                                                    |                                                                                                                                                    
gTAGATGTGCCTGACAATTCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGa                                                                                                                                                       <  1:149997/148‑1 (MQ=255)
gTAGATGTGCCTGACAATTCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGa                                                                                                                                                       <  8:2762/148‑1 (MQ=255)
      gtgCCTGACAATTCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGc                                                                                                                                               >  5:270484/1‑149 (MQ=255)
      gtgCCTGACAATTCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGc                                                                                                                                               >  2:145683/1‑149 (MQ=255)
       tgCCTGACAATTCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCTCGAGGCGCTggc                                                                                                                                              <  1:153644/149‑3 (MQ=255)
       tgCCTAACAATTCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGAGCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGcc                                                                                                                                              >  8:78398/1‑149 (MQ=255)
          cTGACAATTCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTg                                                                                                                                           >  7:129649/1‑149 (MQ=255)
           tGACAATTCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGcc                                                                                                                                              >  3:210042/1‑145 (MQ=255)
           tGACAATTCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGcc                                                                                                                                              <  4:210042/145‑1 (MQ=255)
               aaTTCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGCCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCg                                                                                                                                      >  3:8228/1‑149 (MQ=255)
                aTTCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGc                                                                                                                                     <  2:302336/149‑1 (MQ=255)
                 ttCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCt                                                                                                                                    <  5:371354/149‑1 (MQ=255)
                    aCTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGgc                                                                                                                                 >  2:68716/1‑149 (MQ=255)
                    aCTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGgc                                                                                                                                 <  6:42893/149‑1 (MQ=255)
                     cTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGgcg                                                                                                                                >  7:276765/1‑149 (MQ=255)
                          aGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACgg                                                                                                                           <  5:263404/149‑1 (MQ=255)
                                ggTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGa                                                                                                                     >  3:177156/1‑149 (MQ=255)
                                      ggCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGt                                                                                                               >  2:73881/1‑149 (MQ=255)
                                         catcatCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCAt                                                                                                            <  8:78403/149‑1 (MQ=255)
                                         catcatCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCAt                                                                                                            <  8:276765/149‑1 (MQ=255)
                                          atcatcAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATg                                                                                                           >  2:101860/1‑149 (MQ=255)
                                                      cGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCg                                                                                               >  6:135912/1‑149 (MQ=255)
                                                       ggCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCgg                                                                                              <  4:177156/149‑1 (MQ=255)
                                                        gCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCggg                                                                                             <  8:129649/149‑1 (MQ=255)
                                                             cGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCAt                                                                                                            >  2:43959/1‑129 (MQ=255)
                                                                 cgccgTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGa                                                                                    <  1:73881/149‑1 (MQ=255)
                                                                  gccgTGGCATCAGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGAt                                                                                   <  5:169415/149‑1 (MQ=255)
                                                                    cgTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCTGGGGGGTGGTAATCGGCCGCAATCAAGGCCAGCGCCACGAGCGCTGACTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATgg                                                                                 >  1:360787/1‑149 (MQ=255)
                                                                      tGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGtt                                                                               <  3:363565/149‑1 (MQ=255)
                                                                        gCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCg                                                                             >  5:289180/1‑149 (MQ=255)
                                                                             ggTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCa                                                                        >  5:161359/1‑149 (MQ=255)
                                                                              gTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCGGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAg                                                                       >  2:50896/1‑149 (MQ=255)
                                                                              gTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAg                                                                       >  6:300131/1‑149 (MQ=255)
                                                                                   gTGAGGCCTATGGCAACGGTGGTCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGgc                                                                  >  5:336728/1‑149 (MQ=255)
                                                                                   gTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGTCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGgc                                                                  >  2:146342/1‑149 (MQ=255)
                                                                                        gCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGgccgcc                                                              >  1:111189/1‑148 (MQ=255)
                                                                                             tGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCtt                                                        <  3:88677/149‑1 (MQ=255)
                                                                                               gCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTgg                                                      >  3:84938/1‑149 (MQ=255)
                                                                                               gCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTgg                                                      >  2:203724/1‑149 (MQ=255)
                                                                                               gCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTg                                                       <  6:270484/148‑1 (MQ=255)
                                                                                               gCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTATCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTg                                                       <  6:190922/148‑1 (MQ=255)
                                                                                                cAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGc                                                     <  1:68716/149‑1 (MQ=255)
                                                                                                cAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGc                                                     <  6:289180/149‑1 (MQ=255)
                                                                                                    ggTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTgg                                                 <  4:84938/149‑1 (MQ=255)
                                                                                                                 gtgggAAGCGGCCGCACTCCAGCCCAGCTCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTCCGTCGACTCGCGTTCCACCGATGGTTCGCTGCAGGTGGCCGCCAACCTTGGTGTGGCAGTACTGGCGCt                                    <  6:112237/149‑1 (MQ=255)
                                                                                                                    gtAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTgg                                  <  5:300131/148‑1 (MQ=255)
                                                                                                                    gtAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGa                                 <  3:227477/149‑1 (MQ=255)
                                                                                                                       aTCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACAt                              <  7:78398/149‑1 (MQ=255)
                                                                                                                        tCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGTCGTGGCAGTACTGGCGCTGGACAt                              >  3:40292/1‑148 (MQ=255)
                                                                                                                          ggCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTTGCAGTACTGGCGCTGGACATGAg                           <  1:145683/149‑1 (MQ=255)
                                                                                                                          ggCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACATGa                            <  2:111189/148‑1 (MQ=255)
                                                                                                                          ggCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACATGa                            <  1:101860/148‑1 (MQ=255)
                                                                                                                               cAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACATGAGTATcc                      <  1:203724/149‑1 (MQ=255)
                                                                                                                                 aTGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACATGAGTATCCCg                    >  6:382879/1‑149 (MQ=255)
                                                                                                                                  tGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACATGAGTATCCCGt                   >  4:194699/1‑149 (MQ=255)
                                                                                                                                   acaaggcAGCGCTCGAGGCGCTGGCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACATGAGTATCCCGtt                  >  6:42568/7‑149 (MQ=255)
                                                                                                                                     aaggcAGCGCTCGAGGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACATGa                            >  6:282818/5‑137 (MQ=255)
                                                                                                                                     aaggcAGCGCTCGAGGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACATGa                            <  5:282818/133‑1 (MQ=255)
                                                                                                                                      aggcAGCGCTCGAGGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACATGAGTATCCCGTTCAc               >  1:184690/4‑149 (MQ=255)
                                                                                                                                        ccAGCGCCTCGAGCGCTTCCTTGATTCGCTGGCGCACGGGGCCGACTAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCAGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACATGAGTATCCCGTTCAccg             >  8:163263/1‑149 (MQ=255)
                                                                                                                                        ccAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTAGCTGCAGTTGGCCGCCAACCTTGGCGTTGCAGTACTGGCGCTGGACATGAGTATCCCGTTCAccg             >  8:349500/1‑149 (MQ=255)
                                                                                                                                         cAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTTGCGTGGCAGTACTGGCGCTGGACATGa                            >  5:83223/1‑133 (MQ=255)
                                                                                                                                         cAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACATGa                            <  6:83223/133‑1 (MQ=255)
                                                                                                                                         cAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACATGAGTATCCCGTTCAccgc            >  7:334961/1‑149 (MQ=255)
                                                                                                                                         cAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACATGAGTATCCCGTTCAccgc            >  3:174803/1‑149 (MQ=255)
                                                                                                                                           agcgcTCGAGGCGCTGGCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACATGAGTATCCCGTTCAccgcc           <  5:42568/145‑1 (MQ=255)
                                                                                                                                             gcCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACATGAGTATCCCGTTCACCGCCGCa        >  5:55185/1‑149 (MQ=255)
                                                                                                                                               cTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACATGAGTATCCCGTTCACCGCCGCAcg      >  3:289500/1‑149 (MQ=255)
                                                                                                                                                tCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACATGAGTATCCCGTTCACCGCCGCAcgc     >  6:128605/1‑149 (MQ=255)
                                                                                                                                                   aGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACATGAGTATCCCGTTCACCGCCGCACGCGcc  >  1:223789/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GTAGATGTGCCTGACAATTCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACATGAGTATCCCGTTCACCGCCGCACGCGCC  >  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: