Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,684,806 C→T 6.9% R696H (CGC→CAC)  PP_4988 ← signal transduction histidine kinase CheA/CheY‑like receiver

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,684,8060CT 6.9% 215.0 / 3.5 72R696H (CGC→CAC) PP_4988signal transduction histidine kinase CheA/CheY‑like receiver
Reads supporting (aligned to +/- strand):  ref base C (33/34);  new base T (2/3);  total (35/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCGGGCTGTTGCTGCCAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCATTTGCAGTGGCAGT  >  NC_002947/5684661‑5684952
                                                                                                                                                 |                                                                                                                                                  
gCCGGGCTGTTGCTGCCAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGgcggcg                                                                                                                                                   <  1:283254/147‑1 (MQ=255)
gCCGGGCTGTTGCTGCCAGCGCGACAACAACGTCTCTGAATCGAGCAGGATGTCCCTGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTTCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGgcggcg                                                                                                                                                   <  7:217088/147‑1 (MQ=255)
  cGGGCTGTTGCTGCCAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTcc                                                                                                                                               >  5:21604/1‑149 (MQ=255)
   gggCTGTTGCTGCCAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCa                                                                                                                                              >  4:192782/1‑149 (MQ=255)
    ggCTGTTGCTGCCAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCAc                                                                                                                                             >  4:14210/1‑149 (MQ=255)
        gTTGCTGCCAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGt                                                                                                                                                 >  3:268032/1‑141 (MQ=255)
        gTTGCTGCCAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGt                                                                                                                                                 <  4:268032/141‑1 (MQ=255)
         ttGCTGCCAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCt                                                                                                                                        >  6:54603/1‑149 (MQ=255)
         ttGCTGCCAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCt                                                                                                                                        >  7:125440/1‑149 (MQ=255)
            ctgcCAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGTTGGCCAGGCCGTGGACCGCCGCCTgg                                                                                                                                      <  6:192532/148‑1 (MQ=255)
            ctgcCAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGTTGGCCAGGCCGTGGACCGCCGCCTgg                                                                                                                                      <  6:197041/148‑1 (MQ=255)
            ctgcCAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGc                                                                                                                                     <  5:13982/149‑1 (MQ=255)
             tgcCAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTg                                                                                                                                       >  2:251559/1‑146 (MQ=255)
             tgcCAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTg                                                                                                                                       <  1:251559/146‑1 (MQ=255)
              gcCAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCa                                                                                                                                              >  6:216057/1‑138 (MQ=255)
              gcCAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCa                                                                                                                                              <  5:216057/138‑1 (MQ=255)
               ccAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCaa                                                                                                                                  >  2:107718/1‑149 (MQ=255)
               ccAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCaa                                                                                                                                  >  7:93380/1‑149 (MQ=255)
                cAGCGCGACAACAACGTCTCTGCATCGAGCAGGGTGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAAc                                                                                                                                 <  7:135865/149‑1 (MQ=255)
                cAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAAc                                                                                                                                 <  3:192782/149‑1 (MQ=255)
                cAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAAc                                                                                                                                 <  2:55868/149‑1 (MQ=255)
                    gcgACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGTTGGCCAGGCCGTGGACCGCCGCCTGGc                                                                                                                                     <  8:259438/141‑1 (MQ=255)
                    gcgACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGTTGGCCAGGCCGTGGACCGCCGCCTGGc                                                                                                                                     >  7:259438/1‑141 (MQ=255)
                    gcgACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGgcggcg                                                                                                                                                   <  2:258531/127‑1 (MQ=255)
                    gcgACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGgcggcg                                                                                                                                                   >  1:258531/1‑127 (MQ=255)
                    gcgACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGc                                                                                                                             <  7:253503/149‑1 (MQ=255)
                      gACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTtctc                                                                                                                           >  7:90937/1‑149 (MQ=255)
                         aacaacGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGAt                                                                                                                        <  2:229724/149‑1 (MQ=255)
                           caacGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGAt                                                                                                                        >  3:251746/1‑147 (MQ=255)
                                      cATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCcgcg                                                                                                              <  5:54603/146‑1 (MQ=255)
                                         cGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCg                                                                                                        >  7:202238/1‑149 (MQ=255)
                                            gCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCgggg                                                                                                     <  6:21604/149‑1 (MQ=255)
                                                gATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATgg                                                                                                 <  1:107718/149‑1 (MQ=255)
                                                ctggcccatgccatgGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATgg                                                                                                 >  4:179670/6‑149 (MQ=255)
                                                  tGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATggcg                                                                                               >  5:41608/1‑149 (MQ=255)
                                                   gcccatgccatgGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATggcgg                                                                                              <  3:179670/147‑1 (MQ=255)
                                                   gTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATggcg                                                                                               <  8:125440/148‑1 (MQ=255)
                                                    tCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATggcggc                                                                                             >  3:272017/1‑149 (MQ=255)
                                                      catggcatggGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATggcggc                                                                                             >  2:38719/1‑147 (MQ=255)
                                                      catggcatggGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGcc                                                                                            <  8:182425/148‑1 (MQ=255)
                                                          gcatggGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAgggg                                                                                       >  7:121276/1‑149 (MQ=255)
                                                          gcatggGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAggg                                                                                        <  8:93380/148‑1 (MQ=255)
                                                          gcatggGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAggg                                                                                        <  3:14210/148‑1 (MQ=255)
                                                               ggCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGc                                                                                  >  2:303967/1‑149 (MQ=255)
                                                                gCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAg                                                                                          >  2:180915/1‑140 (MQ=255)
                                                                gCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAg                                                                                          <  1:180915/140‑1 (MQ=255)
                                                                     gAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTTCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGa                                                                            <  8:90937/149‑1 (MQ=255)
                                                                      aaGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAg                                                                           <  3:248845/149‑1 (MQ=255)
                                                                            gCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTgcgc                                                                     <  1:38719/149‑1 (MQ=255)
                                                                              gACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGcc                                                                    >  3:286032/1‑148 (MQ=255)
                                                                                         gATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATgg                                                        <  8:121276/149‑1 (MQ=255)
                                                                                          cggccGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATgg                                                        >  8:17813/5‑148 (MQ=255)
                                                                                           tCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGaa                                                      <  6:311623/149‑1 (MQ=255)
                                                                                               gCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGaacaac                                                  <  3:141042/149‑1 (MQ=255)
                                                                                                              cATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGAcc                                              <  8:101704/138‑1 (MQ=255)
                                                                                                              cATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGAcc                                              >  7:101704/1‑138 (MQ=255)
                                                                                                                tGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAg                                 <  4:290513/149‑1 (MQ=255)
                                                                                                                    gCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGa                                                           >  4:270577/1‑119 (MQ=255)
                                                                                                                    gCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGa                                                           <  3:270577/119‑1 (MQ=255)
                                                                                                                    gCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATg                                                         <  1:58275/121‑1 (MQ=255)
                                                                                                                    gCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATg                                                         >  2:58275/1‑121 (MQ=255)
                                                                                                                        gTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACATGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGccattcca                         >  1:219243/1‑149 (MQ=255)
                                                                                                                           ttCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCAttt                      >  2:148857/1‑149 (MQ=255)
                                                                                                                            tCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCAttt                      <  4:251746/148‑1 (MQ=255)
                                                                                                                            tCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTc                     >  5:198723/1‑149 (MQ=255)
                                                                                                                              aTGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCgg                   >  5:163607/1‑149 (MQ=255)
                                                                                                                              aTGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCg                    <  6:41608/148‑1 (MQ=255)
                                                                                                                                 cAGGTTGGCCAGGCCGTGGACCGCCGCCTGGCCAGCCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCa                >  1:158173/1‑149 (MQ=38)
                                                                                                                                    gCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCAttt             >  8:111818/1‑149 (MQ=255)
                                                                                                                                     ttGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCATTTg            >  5:203338/2‑149 (MQ=255)
                                                                                                                                     cTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCATTTg            >  3:107393/1‑149 (MQ=255)
                                                                                                                                     cTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCATTTg            <  8:202238/149‑1 (MQ=255)
                                                                                                                                        gCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCg                                                                                     <  2:203695/73‑1 (MQ=255)
                                                                                                                                        gCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCg                                                                                     >  1:203695/1‑73 (MQ=255)
                                                                                                                                         ccAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCATTTGCAGt        <  2:219243/149‑1 (MQ=255)
                                                                                                                                          cAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCAttt             <  3:135961/143‑1 (MQ=255)
                                                                                                                                             gcggcgGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCATTTGCAGTgg      >  5:222401/1‑147 (MQ=255)
                                                                                                                                               ggcggTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCATTTGCAGTGGCAGt  <  4:272017/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GCCGGGCTGTTGCTGCCAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCATTTGCAGTGGCAGT  >  NC_002947/5684661‑5684952

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: