Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,684,812 A→G 7.9% V694A (GTG→GCG)  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,8120AG 7.9% 177.1 / 4.8 63V694A (GTG→GCG) PP_4988signal transduction histidine kinase CheA/CheY‑like receiver
Reads supporting (aligned to +/- strand):  ref base A (24/34);  new base G (3/2);  total (27/36)
Fisher's exact test for biased strand distribution p-value = 6.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.85e-01

GGGCTGTTGCTGCCAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCG‑GGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGG‑CCT‑GGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCATTTGCAGTGGCAGTTGATGGCC  >  NC_002947/5684664‑5684960
                                                                                                                                                     |                                                                                                                                                      
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                                                                                                                                    >  6:351797/1‑142 (MQ=255)
                         aacGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCG‑GGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGG‑CCT‑GGCCAACCTGCt                                                                                                                                    <  5:351797/142‑1 (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                                                                                           <  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                                                       >  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)
                                                                                                                                                 gTCCACGG‑CCT‑GGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCATTTGCAGTGGCAGTTGAt      >  2:385489/1‑149 (MQ=255)
                                                                                                                                                  tCCACGG‑CCT‑GGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCATTTGCAGTGGCAGTTGATg     <  1:37350/149‑1 (MQ=255)
                                                                                                                                                    cACGG‑CCT‑GGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGcca                                      >  7:198024/1‑114 (MQ=255)
                                                                                                                                                    cACGG‑CCT‑GGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGcca                                      <  8:198024/114‑1 (MQ=255)
                                                                                                                                                    cACGG‑CCT‑GGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCATTTGCAGTGGCAGTTGATGGc   >  6:251307/1‑149 (MQ=255)
                                                                                                                                                     aCGG‑CCT‑GGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCATTTGCAGTGGCAGTTGATGGcc  <  1:246214/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GGGCTGTTGCTGCCAGCGCGACAACAACGTCTCTGCATCGAGCAGGATGTCCATGGCATGGGCCAGGAAGCTGGCGACCAACTGCGGATCACGCTTGCTGCGGCCGGCATGGTGCG‑GGTCTTCATGCAGGCTGGCCAGGCGGCGGTCCACGG‑CCT‑GGCCAACCTGCTCGATCAGTTCCGCGGCACCGGGGATGGCGGCCAGGGGCGTGCTATCGAGTTGCGCCAAGCCCAGATGGAACAACGACCGCGCAGCCTCCAGCCATTCCATTTCGGCCATTTGCAGTGGCAGTTGATGGCC  >  NC_002947/5684664‑5684960

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