Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,684,806 C→T 8.5% 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 8.5% 160.3 / 3.4 59R696H (CGC→CAC) PP_4988signal transduction histidine kinase CheA/CheY‑like receiver
Reads supporting (aligned to +/- strand):  ref base C (21/32);  new base T (3/2);  total (24/35)
Fisher's exact test for biased strand distribution p-value = 6.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.62e-01

GCCGGGCTGTTGCTGCCAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCG‑GGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGG‑CCT‑GGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCATTTGCAGT  >  NC_002947/5684661‑5684946
                                                                                                                                                  |                                                                                                                                              
gCCGGGCTGTTGCTGCCAGCGCGACAACAACGTCTCTGGATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCG‑GGTCTTCATGCAGGCTGGCCAGgcggcg                                                                                                                                               <  4:82847/147‑1 (MQ=255)
   gggCTGTTGCTGCCAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGCATCACGCTTGCTGCGGCCGGCATGGTGCG‑GGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCa                                                                                                                                          <  5:267789/149‑1 (MQ=255)
        gTTGCTGCCAGCGCGACAACAACGTTTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCG‑GGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGG‑cc                                                                                                                                    >  2:300825/1‑149 (MQ=255)
            ctgcCAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCG‑GGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGG‑CCT‑gg                                                                                                                                <  5:157950/148‑1 (MQ=255)
            ctgcCAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCG‑GGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGG‑CCT‑gg                                                                                                                                <  3:209098/148‑1 (MQ=255)
            ctgcCAGCGCGACAACAACGTCTCTGCATCGAGCAGGAAGTCAATGGAATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCATGCTGAGGCCGGAATGGTGCG‑GGTCTTCATACAGGCTGGCGAGGCGGCGGTCCACGG‑CCT‑GGc                                                                                                                               >  8:198338/1‑149 (MQ=255)
               ccAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCG‑GGTCTTCATGCAGGCTGGCCAGGCCGTGGACCGCCG‑CCT‑GGCCa                                                                                                                             <  2:305521/148‑1 (MQ=255)
                cAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCG‑GGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGG‑CCT‑GGCCAAc                                                                                                                           <  2:256736/149‑1 (MQ=255)
                     cgACAACAACGTCTCTGCATCCAGCAGGATGTCCATCGCATGCGCCAGGAAGCTGCCGACCAACTGCGGATCACGCTTGCTGCGGCCGTCATGGTGCG‑GCTCATCAAGCAGGCTGGCCAGGCGGCGGTCCACGG‑CCT‑GGCCAACCTGCt                                                                                                                      <  8:304925/149‑1 (MQ=255)
                       acaacaACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCG‑GGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGG‑CCT‑GGCCAACCTGCTc                                                                                                                     >  2:117284/1‑148 (MQ=255)
                        caacaaCGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCG‑GGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGG‑CCT‑GGCCAACCTGCTCGa                                                                                                                   >  6:122908/1‑149 (MQ=255)
                         aacaacGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCG‑GGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGG‑CCT‑GGCCAACCTGCTCGAt                                                                                                                  >  2:246214/1‑149 (MQ=255)
                           caacGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCG‑GGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGG‑CCT‑GGCCAACCTGCTCGATCa                                                                                                                <  5:122908/149‑1 (MQ=255)
                            aacGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCG‑GGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGG‑CCT‑GGCCAACCTGCt                                                                                                                      <  5:351797/142‑1 (MQ=255)
                            aacGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCG‑GGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGG‑CCT‑GGCCAACCTGCt                                                                                                                      >  6:351797/1‑142 (MQ=255)
                                          gAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCG‑GGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGG‑CCT‑GGCCAACCTGCTCGATCAGTTCCGCGGCACCgg                                                                                                 >  8:59236/1‑149 (MQ=255)
                                          gAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCG‑GGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGG‑CCT‑GGCCAACCTGCTCGATCAGTTCCGCGGCACCgg                                                                                                 <  2:43439/149‑1 (MQ=255)
                                           aGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCG‑GGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGG‑CCT‑GGCCAACCTGCGCGATCAGTTCCGCGGCACCggg                                                                                                >  8:321319/1‑149 (MQ=255)
                                           aGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCG‑GGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGG‑CCT‑GGCCAACCTGCACGATCAGTTCCGCGGCACCggg                                                                                                >  4:128559/1‑149 (MQ=255)
                                               ggATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCG‑GGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGG‑CCT‑GGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATg                                                                                            >  2:36657/1‑149 (MQ=255)
                                                gATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCG‑GGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGG‑CCT‑GGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATgg                                                                                           <  7:300281/149‑1 (MQ=255)
                                                 aTGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCG‑GGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGG‑CCT‑GGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATggc                                                                                          <  2:111381/149‑1 (MQ=255)
                                                    tCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCG‑GGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGG‑CCT‑GGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATggcggc                                                                                       >  3:57749/1‑149 (MQ=255)
                                                    tCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCG‑GGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGG‑CCT‑GGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATggcgg                                                                                        <  6:377168/148‑1 (MQ=255)
                                                          gcatggGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCG‑GGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGG‑CCT‑GGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAggg                                                                                  <  1:233711/148‑1 (MQ=255)
                                                               ggCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCG‑GGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGG‑CCT‑GGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTg                                                                             <  2:133080/148‑1 (MQ=255)
                                                               ggCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCG‑GGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGG‑CCT‑GGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTg                                                                             <  6:169479/148‑1 (MQ=255)
                                                                gCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCGTGGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGG‑CCT‑GGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGc                                                                            >  1:301424/1‑149 (MQ=255)
                                                                 ccAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCG‑GGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGG‑CCTCGGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCt                                                                           <  1:300825/149‑1 (MQ=255)
                                                                        gCTGGCGAACAACTGGGGATCACGCTTGGTGCGGCCGGCATGGTTCG‑GGTCTTCCTGCATGCTGGCCATGCGGGGGTCCACGGTCCT‑GGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAg                                                                     <  1:117284/148‑1 (MQ=255)
                                                                          tGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCG‑GGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGG‑CCT‑GGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTgc                                                                 >  2:385584/1‑149 (MQ=255)
                                                                            gCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCG‑GGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGG‑CCT‑GGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTgcgc                                                               <  7:179972/149‑1 (MQ=255)
                                                                            gCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCG‑GGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGG‑CCT‑GGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTgcg                                                                <  2:352516/148‑1 (MQ=255)
                                                                              gACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCG‑GGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGG‑CCT‑GGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGcc                                                              >  6:114157/1‑148 (MQ=255)
                                                                               aCCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCG‑GGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGG‑CCT‑GGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCaa                                                            <  2:301424/149‑1 (MQ=255)
                                                                                 cAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCG‑GGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGG‑CCT‑GGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGc                                                          <  4:294686/149‑1 (MQ=255)
                                                                                       cGGATCACGCTTGCTGCGGCCGGCATGGTGCG‑GGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGG‑CCT‑GGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGa                                                     <  6:62038/148‑1 (MQ=255)
                                                                                        ggATCACGCTTGCTGCGGCCGGCATGGTGCG‑GGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGG‑CCT‑GGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATg                                                   >  6:384890/1‑149 (MQ=255)
                                                                                             aCGCTTGCTGCGGCCGGCATGGTGCG‑GGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGG‑CCT‑GGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGaac                                               <  4:57749/148‑1 (MQ=255)
                                                                                             aCGCTTGCTGCGGCCGGCATGGTGCG‑GGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGG‑CCT‑GGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGaac                                               <  6:191199/148‑1 (MQ=255)
                                                                                                    ctgcGGCCGGCATGGTGCG‑GGTCTTCATGCAGGCTTGCCAGGGGGCGGTCCAGGG‑CCT‑GGCCAACATGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGAcc                                        <  7:198338/148‑1 (MQ=255)
                                                                                                       cggccggcATGGTGCG‑GGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGT‑CCT‑GGCCAACCTGCACGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCgcgc                                    <  3:128559/149‑1 (MQ=255)
                                                                                                       cggccggcATGGTGCG‑GGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGG‑CCT‑GGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGATGCGCCAAGCCCAGATGGAACAACGACcgc                                      <  1:36657/147‑1 (MQ=255)
                                                                                                             gcATGGTGCG‑GGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGG‑CCT‑GGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTc                              >  6:47549/1‑149 (MQ=255)
                                                                                                             gcATGGTGCG‑GGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGG‑CCT‑GGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTc                              >  2:37350/1‑149 (MQ=255)
                                                                                                               aTGGTGCG‑GGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGG‑CCT‑GGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGaacaac                                            >  3:212327/1‑133 (MQ=255)
                                                                                                               aTGGTGCG‑GGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGG‑CCT‑GGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGaacaac                                            <  4:212327/133‑1 (MQ=255)
                                                                                                                tGGTGCG‑GGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGG‑CCT‑GGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAg                           <  5:47549/149‑1 (MQ=255)
                                                                                                                  gTGCG‑GGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGG‑CCT‑GGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGcc                         <  7:59236/149‑1 (MQ=255)
                                                                                                                     cG‑GGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGG‑CCT‑GGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGccatt                      >  4:186806/1‑149 (MQ=255)
                                                                                                                              cATGCAGGCTGGCCAGGCGGCGGTCCACGG‑CCT‑GGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCAttt                <  5:384890/147‑1 (MQ=255)
                                                                                                                               aTGCAGGCTGGCCAGGCGGCGGTCCACGG‑CCT‑GGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCg              <  1:385584/148‑1 (MQ=255)
                                                                                                                                tGCAGGCTGGCCAGGCGGCGGTCCACGG‑CCT‑GGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGc            >  6:118621/1‑149 (MQ=255)
                                                                                                                                    ggttGGCCAGGCCGTGGACCGCCG‑CCT‑GGCCAGCCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGAc                                         >  1:180068/4‑116 (MQ=37)
                                                                                                                                    ggttGGCCAGGCCGTGGACCGCCG‑CCT‑GGCCAGCCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGAc                                         <  2:180068/113‑1 (MQ=37)
                                                                                                                                    ggttGGCCAGGCCGTGGACCGCCG‑CCT‑GGCCAGCCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCAtt        >  7:164536/4‑149 (MQ=38)
                                                                                                                                     gttGGCCAGGCCGTGGACCGCCG‑CCT‑GGCCAGCCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCAttt       >  6:265992/3‑149 (MQ=39)
                                                                                                                                     gCTGGCCAGGCGGCGGTCCACGG‑CCT‑GGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCAttt       <  3:186806/149‑1 (MQ=255)
                                                                                                                                         gCCAGGCGGCGGTCCACGG‑CCT‑GGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCATTTGCAg   <  5:114157/149‑1 (MQ=255)
                                                                                                                                          ccAGGCGGCGGTCCACGG‑CCT‑GGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCATTTGCAGt  >  2:38399/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
GCCGGGCTGTTGCTGCCAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCG‑GGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGG‑CCT‑GGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCATTTGCAGT  >  NC_002947/5684661‑5684946

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