Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,684,794 C→T 9.3% S700N (AGC→AAC)  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,7940CT 9.3% 212.6 / 8.0 75S700N (AGC→AAC) PP_4988signal transduction histidine kinase CheA/CheY‑like receiver
Reads supporting (aligned to +/- strand):  ref base C (32/36);  new base T (4/3);  total (36/39)
Fisher's exact test for biased strand distribution p-value = 7.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGCGTCACGCTGGCCGGGCTGTTGCTGCCAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCATTTG  >  NC_002947/5684648‑5684942
                                                                                                                                                  |                                                                                                                                                    
gagCGTCACGCTGGCCGGGCTGTTGCTGCCAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTg                                                                                                                                                    <  8:192521/149‑1 (MQ=255)
 agCGTCACGCTGGCCGGGCTGTTGCTGCCAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTgg                                                                                                                                                   <  4:243724/149‑1 (MQ=255)
 agCGTCACGCTGGCCGGGCTGTTGCTGCCAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTgg                                                                                                                                                   <  5:147778/149‑1 (MQ=255)
     tCACACTGGCCGGGCTGTTGCTGCCAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAg                                                                                                                                               >  1:229724/1‑149 (MQ=255)
       aCGCTGGCCGGGCTGTTGCTGCCAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGTTGGCCAggc                                                                                                                                             >  5:192532/1‑149 (MQ=255)
       aCGCTGGCCGGGCTGTTGCTGCCAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGTTGGCCAggc                                                                                                                                             >  5:197041/1‑149 (MQ=255)
             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:197041/148‑1 (MQ=255)
                         ctgcCAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGTTGGCCAGGCCGTGGACCGCCGCCTgg                                                                                                                            <  6:192532/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                                                                                                                                    <  5:216057/138‑1 (MQ=255)
                           gcCAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCa                                                                                                                                    >  6:216057/1‑138 (MQ=255)
                            ccAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCaa                                                                                                                        >  7:93380/1‑149 (MQ=255)
                            ccAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCaa                                                                                                                        >  2:107718/1‑149 (MQ=255)
                             cAGCGCGACAACAACGTCTCTGCATCGAGCAGGGTGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAAc                                                                                                                       <  7:135865/149‑1 (MQ=255)
                             cAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAAc                                                                                                                       <  2:55868/149‑1 (MQ=255)
                             cAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAAc                                                                                                                       <  3:192782/149‑1 (MQ=255)
                                 gcgACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGTTGGCCAGGCCGTGGACCGCCGCCTGGc                                                                                                                           >  7:259438/1‑141 (MQ=255)
                                 gcgACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGTTGGCCAGGCCGTGGACCGCCGCCTGGc                                                                                                                           <  8:259438/141‑1 (MQ=255)
                                 gcgACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGgcggcg                                                                                                                                         >  1:258531/1‑127 (MQ=255)
                                 gcgACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGgcggcg                                                                                                                                         <  2:258531/127‑1 (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                                                                              <  3:14210/148‑1 (MQ=255)
                                                                       gcatggGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAggg                                                                              <  8:93380/148‑1 (MQ=255)
                                                                            ggCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGc                                                                        >  2:303967/1‑149 (MQ=255)
                                                                             gCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAg                                                                                <  1:180915/140‑1 (MQ=255)
                                                                             gCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAg                                                                                >  2:180915/1‑140 (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                                    >  7:101704/1‑138 (MQ=255)
                                                                                                                           cATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGAcc                                    <  8:101704/138‑1 (MQ=255)
                                                                                                                             tGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAg                       <  4:290513/149‑1 (MQ=255)
                                                                                                                                 gCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGa                                                 >  4:270577/1‑119 (MQ=255)
                                                                                                                                 gCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGa                                                 <  3:270577/119‑1 (MQ=255)
                                                                                                                                 gCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATg                                               >  2:58275/1‑121 (MQ=255)
                                                                                                                                 gCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATg                                               <  1:58275/121‑1 (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)
                                                                                                                                                  cTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCATTTg  >  3:107393/1‑149 (MQ=255)
                                                                                                                                                  cTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCATTTg  <  8:202238/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GAGCGTCACGCTGGCCGGGCTGTTGCTGCCAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGGCCTGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCATTTG  >  NC_002947/5684648‑5684942

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