Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,104,167 A→T 7.0% D232V (GAC→GTC)  PP_2714 → sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,104,1670AT 7.0% 186.3 / 2.7 66D232V (GAC→GTC) PP_2714sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base A (33/28);  new base T (2/3);  total (35/31)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGGCTGGGGGCTGCGGCCACTGCAGAACATGGCCAATGTCATCCGCAGGCGTCATGCAGATTCGCTGGAGCCGCTGCAACTGGTACCGCTT‑‑CCCACTGAACTGGAGC‑‑CCATGCAGGCGGCCATCAATCGCTTGCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGG  >  NC_002947/3104020‑3104314
                                                                                                                                                       |                                                                                                                                                   
tCGGCTGGGGGCTGCGGCCACTGCAGAACATGGCCAATGTCATCCGCAGGCGTCATGCAGATTCGCTGGAGCCGCTGCAACTGGTACCGCTT‑‑CCCACTGAACTGGAGC‑‑CCATGCAGGCGGCCATCAAACGCTTGCTGGCGCAGATcgag                                                                                                                                                    >  3:127173/1‑148 (MQ=255)
tCGGCTGGGGGCTGCGGCCACTGCAGAACATGGCCAATGTCATCCGCAGGCGTCATGCAGATTCGCTGGAGCCGCTGCAACTGGTACCGCTT‑‑CCCACTGAACTGGAGA‑‑CCATGCAGGCGGCCATCAATCGCTTGCTGGCGCAGATcgac                                                                                                                                                    >  4:21826/1‑149 (MQ=255)
  ggCTGGGGGCTGCGGCCACTGCAGAACATGGCCAATGTCATCCGCAGGCGTCATGCAGATTCGCTGGAGCCGCTGCAACTGGTACCGCTT‑‑CCCACTGAACTGGAGC‑‑CCATGCAGGCGGCCATCAATCGCTTGCTGGCGCAGATcgacg                                                                                                                                                   <  2:266814/148‑1 (MQ=255)
    cTGGGGGCTGCGGCCACTGCAGAACATGGCCAATGTCATCCGCAGGCGTCATGCAGATTCGCTGGAGCCGCTGCAACTGGTACCGCTT‑‑CCCACTGAACTGGAGC‑‑CCATGCAGGCGGCCATCAATCGCTTGCTGGCGCAGATCGACGATc                                                                                                                                                >  8:175695/1‑149 (MQ=255)
     tGGGGGCTGCGGCCACTGCAGAACATGGCCAATGTCATCCGCAGGCGTCATGCAGATTCGCTGGAGCCGCTGCAACTGGTACCGCTT‑‑CCCACTGAACTGGAGC‑‑CCATGCAGGCGGCCATCAATCGCTTGCTGGCGCAGATCGACGATct                                                                                                                                               >  4:114242/1‑149 (MQ=255)
       ggggCTGCGGCCACTGCAGAACATGGCGAATCGCATCCGCAGGCGTCATGCAGATTCGCTGGAGCCGCTGCAACTGGTACCGCTT‑‑CCCACTGAACTGGAGC‑‑CCATGCAGGCGGCCATCAATCGCTTGCTGGCGCAGATCGACGATctg                                                                                                                                              <  6:282482/148‑1 (MQ=255)
       ggggCTGCGGCCACTGCAGAACATGGCCAATGTCATCCGCAGGCGTCATGCAGATTCGCTGGAGCCGCTGCAACTGGTACCGCTT‑‑CCCACTGAACTGGAGC‑‑CCATGCAGGCGGCCATCAATCGCTTGCTGGCGCAGATCGACGATctg                                                                                                                                              <  6:110813/148‑1 (MQ=255)
         ggCTGCGGCCACTGCAGAACATGGCCAATGTCATCCGCAGGCGTCATGCAGATTCGCTGGAGCCGCTGCAACTGGTACCGCTTCCCCCTCTGAACTGGAGC‑‑CCATGCAGGCGGCCATCAATCGCTTGCTGGCGCAGATCGTCGATctgc                                                                                                                                             <  1:298892/149‑1 (MQ=255)
                  cACTGCAGAACATGGCCAATGTCATCCGCAGGCGTCATGCAGATTCGCTGGAGCCGCTGCAACTGGTACCGCTT‑‑CCCACTGAACTGGAGC‑‑CGATGCAGGCGGCCATCAATCGCTTGCTGGCGCAGATCGACGATCTGCTGCGGCGTGaa                                                                                                                                  <  3:21826/149‑1 (MQ=255)
                           aCATGGCCAATGTCATCCGCAGGCGTCATGCAGATTCGCTGGAGCCGCTGCAACTGGTACCGCTT‑‑CCCACTGAACTGGCGC‑‑CCATGCAGGCGGCCATCAATCGCTTGCTGGCGCAGATCGTCGATctgcgc                                                                                                                                           <  8:84373/131‑3 (MQ=255)
                           aCATGGCCAATGTCATCCGCAGGCGTCATGCAGATTCGCTGGAGCCGCTGCAACTGGTACCGCTT‑‑CCCACTGAACTGGAGC‑‑CCATGCAGGCGGCCATCAATCGCTTGCTGGCGCAGATCGTCGATctgcgc                                                                                                                                           >  7:84373/1‑129 (MQ=255)
                                 ccAATGTCATCCGCAGGCGTCATGCAGATTCGCTGGAGCCGCTGCAACTGGTACCGCTT‑‑CCCACTGAACTGGAGC‑‑CCATGCAGGCGGCTATCAATCGCTTGCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGcc                                                                                                                   >  4:209798/1‑149 (MQ=255)
                                       tCATCCGCAGGCGTCATGCAGATTCGCTGGAGCCGCTGCAACTGGTACCGCTT‑‑CCCACTGAACTGGAGC‑‑CCATGCAGGCGGCCATCAATCGCTTGCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGAtgct                                                                                                             >  7:70697/1‑149 (MQ=255)
                                             gCAGGCGTCATGCAGATTCGCTGGAGCCGCTGCAACTGGTACCGCTT‑‑CCCACTGAACTGGAGCAGGCATGCAGGCGGCCATCAATCGCTTGCTGGCGCAGAACCACGATCTGATGCGGCTTGAACATCGCTTCATTGCCCATCCTGcc                                                                                                          >  1:303016/1‑148 (MQ=255)
                                              cAGGCGTCATGCAGATTCGCTGGAGCCGCTGCAACTGGTACCGCTT‑‑CCCACTGAACTGGAGC‑‑CCATGCAGGCGGCCATCAATCGCTTGCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACg                                                                                                      <  4:97454/149‑1 (MQ=255)
                                                 gCGTCATGCAGATTCGCTGGAGCCGCTGCAACTGGTACCGCTT‑‑CCCACTGAACTGGAGC‑‑CCATGCAGGCGGCCATCAATCGCTTGCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCGGCCCACGaaa                                                                                                   >  8:260546/1‑149 (MQ=255)
                                                    tCATGCAGATTCGCTGGAGCCGCTGCAACTGGTACCGCTT‑‑CCCACTGAACTGGAGC‑‑CCATGCAGGCGGCCATCAATCGCTTGCTGGCGCAGATCGTCGATctgcgc                                                                                                                                           <  8:74542/106‑3 (MQ=255)
                                                    tCATGCAGATTCGCTGGAGCCGCTGCAACTGGTACCGCTT‑‑CCCACTGAACTGGAGC‑‑CCATGCAGGCGGCCATCAATCGCTTGCTGGCGCAGATCGTCGATctgcgc                                                                                                                                           >  7:74542/1‑104 (MQ=255)
                                                      aTGCAGATTCGCTGGAGCCGCTGCAACTGGTACCGCTT‑‑CCCACTGAACTGGAGC‑‑CCATGCAGGCGGCCATCAATCGCTTGCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCg                                                                                               >  6:269203/1‑148 (MQ=255)
                                                            aTTCGCTGGAGCCGCTGCAACTGGTACCGCTT‑‑CCCACTGAACTGGAGC‑‑CCATGCAGGCGGCCATCAATCGCTTGCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCct                                                                                        >  5:310788/1‑149 (MQ=255)
                                                             ttCGCTGGAGCCGCTGCAACTGGTACCGCTT‑‑CCCACTGAACTGGAGC‑‑CCATGCAGGCGGCCATCAATCGCTTGCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGcc                                                                                                                   <  6:292117/121‑1 (MQ=255)
                                                             ttCGCTGGAGCCGCTGCAACTGGTACCGCTT‑‑CCCACTGAACTGGAGC‑‑CCATGCAGGCGGCCATCAATCGCTTGCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGcc                                                                                                                   >  5:292117/1‑121 (MQ=255)
                                                                       ccGCTGCAACTGGTACCGCTT‑‑CCCACTGACCTGGCGC‑‑CCATGCAGGCGGCCATCAATCGCTTGCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCGAGGAAATGCGTACGCCTCTGGCCGTAtt                                                                             <  5:309681/149‑1 (MQ=255)
                                                                         gCTGCAACTGGTACCGCTT‑‑CCCACTGAACTGGAGC‑‑CCATGCAGGCGGCCATCAATCGCTTGCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTgc                                                                           <  3:114242/149‑1 (MQ=255)
                                                                           tGCAACTGGTACCGCTT‑‑CCCACTGAACTGGAGC‑‑CCATGCAGGCGGCCATCAATCGCTTGCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTgcgc                                                                         >  8:214154/1‑149 (MQ=255)
                                                                               aCTGGTACCGCTT‑‑CCCACTGAACTGGAGC‑‑CCATGCAGGCGGCCATCAATCGCTTGCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTg                                                                      >  2:196718/1‑148 (MQ=255)
                                                                               aCTGGTACCGCTT‑‑CCCACTGAACTGGAGC‑‑CCATGCAGGCGGCCATCAATCGCTTGCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGATTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGc                                                                     <  2:303016/149‑1 (MQ=255)
                                                                                   gTACCGCTT‑‑CCCAGTGAACTGGAGC‑‑CGATGCAGGCGGCCATCAATCGCTTGCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACTCCTCTGGCCGTATTGCGCCTGCAcgc                                                                 <  4:127173/149‑1 (MQ=255)
                                                                                   gTACCGCTT‑‑CCCACTGAACTGGAGC‑‑CCATGCAGGCGGCCATCAATCGCTTGCTGGCGCAGATCGACGATCTGCTgcggc                                                                                                                                       >  6:277907/1‑79 (MQ=255)
                                                                                   gTACCGCTT‑‑CCCACTGAACTGGAGC‑‑CCATGCAGGCGGCCATCAATCGCTTGCTGGCGCAGATCGACGATCTGCTgcggc                                                                                                                                       <  5:277907/79‑1 (MQ=255)
                                                                                   gTACCGCTT‑‑CCCACTGAACTGGAGC‑‑CCATGCAGGCGGCCATCAATCGCTTGCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCAcgc                                                                 >  5:274620/1‑149 (MQ=255)
                                                                                   gTACCGCTT‑‑CCCACTGAACTGGAGC‑‑CCATGCAGGCGGCCATCAATCGCTTGCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCAcgc                                                                 <  2:188618/149‑1 (MQ=255)
                                                                                    tACCGCTT‑‑CCCACTGAACTGGAGC‑‑CCATGCAGGCGGCCATCAATCGCTTGCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCAcgcg                                                                <  7:175695/149‑1 (MQ=255)
                                                                                      ccGCTT‑‑CCCACTGAACTGGAGC‑‑CCATGCAGGCGGCCATCAATCGCTTGCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTg                                                                                     <  5:205453/126‑1 (MQ=255)
                                                                                      ccGCTT‑‑CCCACTGAACTGGAGC‑‑CCATGCAGGCGGCCATCAATCGCTTGCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTg                                                                                     >  6:205453/1‑126 (MQ=255)
                                                                                       cGCTT‑‑CCCACTGAACTGGAGC‑‑CCATGCAGGCGGCCATCAATCGCTTGCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCAcgcg                                                                >  7:15454/1‑146 (MQ=255)
                                                                                               ccACTGAACTGGAGC‑‑CCATGCAGGCGGCCATCAATCGCTTGCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAAcgcg                                                       >  7:173713/1‑149 (MQ=255)
                                                                                               ccACTGAACTGGAGC‑‑CCATGCAGGCGGCCATCAATCGCTTGCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCAGTATTGCGCCTGCACGCGCAGAAcgcg                                                       >  8:307029/1‑149 (MQ=255)
                                                                                                  cTGAACTGGAGC‑‑CCATGCAGGCGGCCATCAATCGCTTGCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTg                                                    >  4:351417/1‑149 (MQ=255)
                                                                                                   tGAACTGGAGC‑‑CCATGCAGGCGGCCATCAATCGCTTGCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTgg                                                   >  7:210407/1‑149 (MQ=255)
                                                                                                     aaCTGGAGC‑‑CCATGCAGGCGGCCATCAATCGCTTGCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGgcc                                                 >  8:172320/1‑149 (MQ=255)
                                                                                                                 cATGCAGGCGGCCATCAATCGCTTGCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCa                                                                                                                        >  4:31490/1‑68 (MQ=255)
                                                                                                                 cATGCAGGCGGCCATCAATCGCTTGCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCa                                                                                                                        <  3:31490/68‑1 (MQ=255)
                                                                                                                 cATGCAGGCGGCCATCAATCGCTTGCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCGGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAat                                        >  3:15683/1‑148 (MQ=255)
                                                                                                                 cATGCAGGCGGCCATCAATCGCTTGCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAat                                        >  1:354169/1‑148 (MQ=255)
                                                                                                                  aTGCAGGCGGCCATCAATCGCTTGCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTg                                                                                                                     >  5:68896/1‑70 (MQ=255)
                                                                                                                  aTGCAGGCGGCCATCAATCGCTTGCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTg                                                                                                                     <  6:68896/70‑1 (MQ=255)
                                                                                                                    gCAGGCGGCCATCAATCGCTTGCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCAcg                                                                  <  1:354721/119‑1 (MQ=255)
                                                                                                                    gCAGGCGGCCATCAATCGCTTGCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCAcg                                                                  >  2:354721/1‑119 (MQ=255)
                                                                                                                     cAGGCGGCCATCAATCGCTTGCTGGCGCAGATCGACGTTCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAat                                        <  6:363751/144‑1 (MQ=255)
                                                                                                                     cAGGCGGCCATCAATCGCTTGCTGGCGCAGATCGACGTTCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAat                                        >  5:363751/1‑144 (MQ=255)
                                                                                                                      aGGCGGCCATCAATCGCTTGCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGcc                                  >  5:346654/1‑149 (MQ=255)
                                                                                                                        gcggcCATCAATCGCTTGCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCAGGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGcc                                  <  8:173713/147‑1 (MQ=255)
                                                                                                                        gcggcCATCAATCGCTTGCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGa                                <  6:274620/149‑1 (MQ=255)
                                                                                                                         cggcCATCAATCGCTTGCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAg                               <  5:269203/149‑1 (MQ=255)
                                                                                                                              aTCAATCGCTTGCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAgcgc                            >  7:286587/1‑147 (MQ=255)
                                                                                                                                  aTCGCTTGCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGaaa                      <  3:209798/149‑1 (MQ=255)
                                                                                                                                  aTCGCTTGCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGa                        >  3:66887/1‑147 (MQ=255)
                                                                                                                                   tCGCTTGCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCGGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAg                     >  3:166024/1‑149 (MQ=255)
                                                                                                                                         gCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCt                 <  6:310788/147‑1 (MQ=255)
                                                                                                                                         gCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCt                 <  8:15454/147‑1 (MQ=255)
                                                                                                                                         gCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCt                 <  5:51805/147‑1 (MQ=255)
                                                                                                                                         gCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTgg               >  2:173608/1‑149 (MQ=255)
                                                                                                                                          cTGTCGCAGATCGACGATCTGCTGCGGCGTGAACATAGCTTCATTGCCGATTCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTgg               <  7:260546/148‑1 (MQ=255)
                                                                                                                                           tGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGc             >  6:63507/1‑149 (MQ=255)
                                                                                                                                              cgcAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCttt          <  3:351417/149‑1 (MQ=255)
                                                                                                                                               gcAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTc         >  7:352935/1‑149 (MQ=255)
                                                                                                                                               gcAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTc         <  8:70697/149‑1 (MQ=255)
                                                                                                                                                   aTCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCGGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTtgt     <  4:166024/149‑1 (MQ=255)
                                                                                                                                                      gacgaTCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTggg  >  5:314759/1‑149 (MQ=255)
                                                                                                                                                       |                                                                                                                                                   
TCGGCTGGGGGCTGCGGCCACTGCAGAACATGGCCAATGTCATCCGCAGGCGTCATGCAGATTCGCTGGAGCCGCTGCAACTGGTACCGCTT‑‑CCCACTGAACTGGAGC‑‑CCATGCAGGCGGCCATCAATCGCTTGCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGG  >  NC_002947/3104020‑3104314

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: