Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,552,394 Δ1 bp 7.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. 7.5% 213.1 / 10.5 67coding (51/978 nt)PP_3137group 2 family glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base C (31/31);  new base . (3/2);  total (34/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.15e-01

CTGGTAGATGTGCCTGACAATTCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACATGAGTATCCCGTTCACCGCCGCA  >  NC_002947/3552249‑3552541
                                                                                                                                                 |                                                                                                                                                    
cTGGTAGATGTGCCTGACAATTCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTc                                                                                                                                                   <  6:143518/149‑1 (MQ=255)
cTGGTAGATGTGCCTAACAATTCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTc                                                                                                                                                   <  1:50861/149‑1 (MQ=255)
   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                                                                                                                                        <  4:210042/145‑1 (MQ=255)
              tGACAATTCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGcc                                                                                                                                        >  3:210042/1‑145 (MQ=255)
                  aaTTCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGCCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCg                                                                                                                                >  3:8228/1‑149 (MQ=255)
                   aTTCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGc                                                                                                                               <  2:302336/149‑1 (MQ=255)
                    ttCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCt                                                                                                                              <  5:371354/149‑1 (MQ=255)
                       aCTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGgc                                                                                                                           <  6:42893/149‑1 (MQ=255)
                       aCTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGgc                                                                                                                           >  2:68716/1‑149 (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                                                >  2:203724/1‑149 (MQ=255)
                                                                                                  gCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTgg                                                >  3:84938/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      >  3:174803/1‑149 (MQ=255)
                                                                                                                                            cAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACATGAGTATCCCGTTCAccgc      >  7:334961/1‑149 (MQ=255)
                                                                                                                                                gcCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACATGAGTATCCCGTTCACCGCCGCa  >  5:55185/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CTGGTAGATGTGCCTGACAATTCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGAGCGCTGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACATGAGTATCCCGTTCACCGCCGCA  >  NC_002947/3552249‑3552541

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: