Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,684,804 G→C 8.5% R697G (CGC→GGC)  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,8040GC 8.5% 169.7 / 7.3 59R697G (CGC→GGC) PP_4988signal transduction histidine kinase CheA/CheY‑like receiver
Reads supporting (aligned to +/- strand):  ref base G (21/33);  new base C (3/2);  total (24/35)
Fisher's exact test for biased strand distribution p-value = 3.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.48e-01

GCTGGCCGGGCTGTTGCTGCCAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCG‑GGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGG‑CCT‑GGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCATTTGCAGT  >  NC_002947/5684657‑5684946
                                                                                                                                                    |                                                                                                                                                
gCTGGCCGGGCTGTTGCTGCCAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCG‑GGTCTTCATGCAGGCTGGCCAggcgg                                                                                                                                                 <  3:329933/149‑1 (MQ=255)
    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                                                                                                                                <  3:209098/148‑1 (MQ=255)
                ctgcCAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCG‑GGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGG‑CCT‑gg                                                                                                                                <  5:157950/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                                                                                                 <  2:43439/149‑1 (MQ=255)
                                              gAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCG‑GGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGG‑CCT‑GGCCAACCTGCTCGATCAGTTCCGCGGCACCgg                                                                                                 >  8:59236/1‑149 (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                                                                             <  6:169479/148‑1 (MQ=255)
                                                                   ggCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCG‑GGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGG‑CCT‑GGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTg                                                                             <  2:133080/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                                            <  4:212327/133‑1 (MQ=255)
                                                                                                                   aTGGTGCG‑GGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGG‑CCT‑GGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGaacaac                                            >  3:212327/1‑133 (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                                         <  2:180068/113‑1 (MQ=37)
                                                                                                                                        ggttGGCCAGGCCGTGGACCGCCG‑CCT‑GGCCAGCCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGAc                                         >  1:180068/4‑116 (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)
                                                                                                                                                    |                                                                                                                                                
GCTGGCCGGGCTGTTGCTGCCAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCG‑GGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGG‑CCT‑GGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCATTTGCAGT  >  NC_002947/5684657‑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: ‑