Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,493,666 G→A 9.1% E399E (GAG→GAA)  trg → methyl‑accepting chemotaxis protein III, ribose and galactose sensor receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,493,6660GA 9.1% 140.5 / 4.8 44E399E (GAG→GAA) trgmethyl‑accepting chemotaxis protein III, ribose and galactose sensor receptor
Reads supporting (aligned to +/- strand):  ref base G (17/23);  new base A (2/2);  total (19/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GGTGTAGTAAAAACGATGGGCGCTATCTCCACGAGTTCGAAGAAAATTTCTGAGATCACCGCCGTCATCAACAGTATTGCTTTCCAGACGAATATTCTGGCACTGAATGCTGCCGTTGAAGCCGCGCGAGCGGGTGAGCAAGGGCGTGGATTTGCCGTTGTCGCCAGCGAAGTACGGACACTCGCAAGTCGCAGCGCTCAGGCGGCGAAAGAGATTGAAGGCTTGATCAGTGAATCAGTCAGGTTAATTGACCTGGGGTCGGATGAGGTGGCAAC  >  NC_000913/1493529‑1493803
                                                                                                                                         |                                                                                                                                         
ggTGTAGTAAAAACGATGGGCGCTATCTCCACGAGTTCGAAGAAAATTTCTGAGATCACCGCCGTCATCAACAGTATTGCTTTCCAGACGAATATTCTGGCACTGAATGCTGCCGTTGAAGCCGCGCGAGCGGGTGAGCa                                                                                                                                         <  1:188849/140‑1 (MQ=255)
            aCGATGGGCGCTATCTCCACGAGTTCGAAGAAAATTTCTGAGATCACCGCCGTCATCAACAGTATTGCTTTCCAGACGAATATTCTGGCACTGAATGCTGCCGTTGAAGCCGCGCGAGCGGGTGAGCAAGGGCGTGGAtt                                                                                                                             <  2:952739/140‑1 (MQ=255)
            aCGATGGGCGCTATCTCCACGAGTTCGAAGAAAATTTCTGAGATCACCGCCGTCATCAACAGTATTGCTTTCCAGACGAATATTCTGGCACTGAATGCTGCCGTTGAAGCCGCGCGAGCGGGTGAGCAAGGGCGTGGAtt                                                                                                                             <  2:952744/140‑1 (MQ=255)
              gATGGGCGCTATCTCCACGAGTTCGAAGAAAATTTCTGAGATCACCGCCGTCATCAACAGTATTGCTTTCCAGACGAATATTCTGGCACTGAATGCTGCCGTTGAAGCCGCGCGAGCGGGTGAGCAAGGGCGTGGATTTg                                                                                                                           >  1:120778/1‑140 (MQ=255)
                 gggCGCTATCTCCACGAGTTCGAAGAAAATTTCTGAGATCACCGCCGTCATCAACAGTATTGCTTTCCAGACGAATATTCTGGCACTGAATGCTGCCGTTGAAGCCGCGCGAGCGGGTGAGCAAGGGCGTGGATTTGCCg                                                                                                                        <  1:908897/140‑1 (MQ=255)
                               cGAGTTCGATGATAATTTCTGAGATCACCGCCGTCATCAACAGTATTGCTTTCCAGACGAATATTCTGGCACTGAATGCTGCCGTTGAAGCCGCGCGAGCGGGTGAGCAAGGGCGTGGATTTGCCGTTGTCGCCAGCGaa                                                                                                          <  2:217435/140‑1 (MQ=255)
                               cGAGTTCGAAGAAAATTTCTGAGATCACCGCCGTCATCAACAGTATTGCTTTCCAGACGAATATTCTGGCACTGAATGCTGCCGTTGAAGCCGCGCGAGCGGGTGAGCAAGGGCGTGGATTTGCCGTTGTCGCCAGCGaa                                                                                                          <  2:217327/140‑1 (MQ=255)
                               cGAGTTCGAAGAAAATTTCTGAGATCACCGCCGTCATCAACAGTATTGCTTTCCAGACGAATATTCTGGCACTGAATGCTGCCGTTGAAGCCGCGCGAGCGGGTGAGCAAGGGCGTGGATTTGCCGTTGTCGCCAGCGaa                                                                                                          <  2:845694/140‑1 (MQ=255)
                               cGAGTTCGAAGAAAATTTCTGAGAACACCGCCGTCATCAACAGTATTGCTTTCCAGACGAATATTCTGGCACTGAATGCTGCCGTTGAAGCCGCGCGAGCGGGTGAGCAAGGGCGTGGATTTGCCGTTGTCGCCAGCGaa                                                                                                          <  2:217340/140‑1 (MQ=255)
                                   ttCGAAGAAAATTTCTGAGATCACCGCCGTCAACAACAGTATTGCTTTCCAGACGTATATTCTGGCACTGTATGCTGCCGTTGAAGCCGCGCGAGCGGGTGAGCTAGGGCGTGGATTTGCCGTAGACGCCAGCGAAGTAc                                                                                                      >  2:1215048/1‑140 (MQ=255)
                                     cGAAGAAAATTTCTGAGATCACCGCCGTCATCAACAGTATTGCTTTCCAGACGAATATTCTGGCACTGAATGCTGCCGTTGAAGCCGCGCGAGCGGGTGAGCAAGGGCGTGGATTTGCCGTTGTCGCCAGCGACGTACgg                                                                                                    >  2:382462/1‑140 (MQ=255)
                                     cGAAGAAAATTTCTGAGATCACCGCCGTCATCAACAGTATTGCTTTCCAGACGAATATTCTGGCACTGAATGCTGCCGTTGAAGCCGCGCGAGCGGGTGAGCAAGGGCGTGGATTTGCCGTTGTCGCCAGCGAAGTACgg                                                                                                    >  2:47778/1‑140 (MQ=255)
                                     cGAAGAAAATTTCTGAGATCACCGCCGTCATCAACAGTATTGCTTTCCAGACGAATATTCTGGCACTGAATGCTGCCGTTGAAGCCGCGCGAGCGGGTGAGCAAGGGCGTGGATTTGCCGTTGTCGCCAGCGAAGTACgg                                                                                                    >  2:47806/1‑140 (MQ=255)
                                     cGAAGAAAATTTCTGAGATCACCGCCGTCATCAACAGTATTGCTTTCCAGACGAATATTCTGGCACTGAATGCTGCCGTTGAAGCCGCGCGAGCGGGTGAACAAGGGCGTGGATTTGCCGTTGTCGCCAGCGAAGTACgg                                                                                                    >  2:1042278/1‑140 (MQ=255)
                                     cGAAGAAAATTTCTGAGATCACCGCCGTCATCAACAGTATTGCTTTCCAGACGAATATTCTGGCACTGAATGCTGCCGTTGAAGCCGCGCGAGCGGGTGAACAAGGGCGTGGATTTGCCGTTGTCGCCAGCGAAGTACgg                                                                                                    >  2:1042346/1‑140 (MQ=255)
                                                   gagaTCACCGCCGTCATCAACAGTATTGCTTTCCAGACGAATATTCTGGCACTGAATGCTGCCGTTGAAGCCGCGCGAGCGGGTGAGCAAGGGCGTGGATTTGCCGTTGTCGCCAGCGAAGTACGGACACTCGCAAGTCg                                                                                      <  2:964598/140‑1 (MQ=255)
                                                   gagaTCACCGCCGTCATCAACAGTATTGCTTTCCAGACGAATATTCTGGCACTGAATGCTGCCGTTGAAGCCGCGCGAGCGGGTGAGCAAGGGCGTGGATTTGCCGTTGTCGCCAGCGAAGTACGGACACTCGCAAGTCg                                                                                      <  2:964677/140‑1 (MQ=255)
                                                    agaTCACCGCCGTCATCAACAGTATTGCTTTCCAGACGAATATTCTGGCACTGAATGCTGCCGTTGAAGCCGCGCGAGCGGGTGAGCAAGGGCGTGGATTTGCCGTTGTCGCCAGCGAAGTACGGACACTCGCAAGTCGc                                                                                     <  1:1190079/140‑1 (MQ=255)
                                                      aTCACCGCCGTCATCAACAGTATTGCTTTCCAGACGAATATTCTGGCACTGAATGCTGCCGTTGAAGCCGCGCGAGCGGGTGAGCAAGGGCGTGGATTTGCCGTTGTCGCCAGCGAAGTACGGACACTCGCAAGTCGCAg                                                                                   >  1:253577/1‑140 (MQ=255)
                                                          ccgccgTCATCAACAGTATTGCTTTCCAGACGAATATTCTGGCACTGAATGCTGCCGTTGAAGCCGCGCGAGCGGGTGAGCAAGGGCGTGGATTTGCCGTTGTCGCCAGCGAAGTACGGACACTCGCAAGTCGCAGCGCt                                                                               <  1:815672/140‑1 (MQ=255)
                                                                     aaCAGTATTGCTTTCCAGACGAATATTCTGGCACTGAATGCTGCCGTTGAAGCCGCGCGAGCGGGTGAGCAAGGGCGTGGATTTGCCGTTGTCGCCAGCGAAGTACGGACACTCGCAAGTCGCAGCGCTCAGGCGGCGaa                                                                    >  1:335646/1‑140 (MQ=255)
                                                                       cAGTATTGCTTTCCAGACGAATATTCTGGCACTGAATGCTGCCGTTGAAGCCGCGCGAGCGGGTGAGCAAGGGCGTGGATTTGCCGTTGTCGCCAGCGAAGTACGGACACTCGCAAGTCGCAGCGCTCAGGCGGCGAAag                                                                  <  2:458032/140‑1 (MQ=255)
                                                                       cAGTATTGCTTTCCAGACGAATATTCTGGCACTGAATGCTGCCGTTGAAGCCGCGCGAGCGGGTGAGCAAGGGCGTGGATTTGCCGTTGTCGCCAGCGAAGTACGGACACTCGCAAGTCGCAGCGCTCAGGCGGCGAAag                                                                  <  2:597131/140‑1 (MQ=255)
                                                                       cAGTATTGCTTTCCAGACGAATATTCTGGCACTGAATGCTGCCGTTGAAGCCGCGCGAGCGGGTGAGCAAGGGCGTGGATTTGCCGTTGTCGCCAGCGAAGTACGGACACTCGCAAGTCGCAGCGCTCAGGCGGCGAAag                                                                  <  2:711397/140‑1 (MQ=255)
                                                                       cAGTATTGCTTTCCAGACGAATATTCTGGCACTGAATGCTGCCGTTGAAGCCGCGCGAGCGGGTGAGCAAGGGCGTGGATTTGCCGTTGTCGCCAGCGAAGTACGGACACTCGCAAGTCGCAGCGCTCAGGCGGCGAAag                                                                  <  2:711641/140‑1 (MQ=255)
                                                                       cAGTATTGCTTTCCAGACGAATATTCTGGCACTGAATGCTGCCGTTGAAGCCGCGCGAGCGGGTGAGCAAGGGCGTGGATTTGCCGTTGTCGCCAGCGAAGTACGGACACTCGCAAGTCGCAGCGCTCAGGCGGCGAAag                                                                  <  2:986446/140‑1 (MQ=255)
                                                                          tATTGCTTTCCAGACGAATATTCTGGCACTGAATGCTGCCGTTGAAGCCGCGCGAGCGGGTGAGCAAGGGCGTGGATTTGCCGTTGTCGCCAGCGAAGTACGGACACTCGCAAGTCGCAGCGCTCAGGCGGCGAAAgaga                                                               >  1:26189/1‑140 (MQ=255)
                                                                            ttGCTTTCCAGACGAATATTCTGGCACTGAATGCTGCCGTTGAAGCCGCGCGAGCGGGTGAGCAAGGGCGTGGATTTGCCGTTGTCGCCAGCGAAGTACGGACACTCGCAAGTCGCAGCGCTCAGGCGGCGAAAGAGAtt                                                             >  1:751267/1‑140 (MQ=255)
                                                                             tGCTTTCCAGACGAATATTCTGGCACTGAATGCTGCCGTTGAAGCCGCGCGAGCGGGTGAACAAGGGCGTGGATTTGCCGTTGTCGCCAGCGAAGTACGGACACTCGCAAGTCGCAGCGCTCAGGCGGCGAAAGAGATTg                                                            <  1:1042346/140‑1 (MQ=255)
                                                                             tGCTTTCCAGACGAATATTCTGGCACTGAATGCTGCCGTTGAAGCCGCGCGAGCGGGTGAACAAGGGCGTGGATTTGCCGTTGTCGCCAGCGAAGTACGGACACTCGCAAGTCGCAGCGCTCAGGCGGCGAAAGAGATTg                                                            <  1:1042278/140‑1 (MQ=255)
                                                                               cTTTCCAGACGAATATTCTGGCACTGAATGCTGCCGTTGAAGCCGCGCGAGCGGGTGAGCAAGGGCGTGGATTTGCCGTTGTCGCCAGCGAAGTACGGACACTCGCAAGTCGCAGCGCTCAGGCGGCGAAAGAGATTGaa                                                          <  2:152296/140‑1 (MQ=255)
                                                                                    cAGACGAATATTCTGGCACTGAATGCTGCCGTTGAAGCCGCGCGAGCGGGTGAGCAAGGGCGTGGATTTGCCGTTGTCGCCAGCGAAGTACGGACACTCGCAAGTCGCAGCGCTCAGGCGGCGAAAGAGATTGAAGGCtt                                                     >  1:1142638/1‑140 (MQ=255)
                                                                                      gACGAATATTCTGGCACTGAATGCTGCCGTTGAAGCCGCGCGAGCGGGTGAGCAAGGGCGTGGATTTGCCGTTGTCGCCAGCGAAGTACGGACACTCGCAAGTCGCAGCGCTCAGGCGGCGAAAGAGATTGAAGGCTTGa                                                   >  2:429142/1‑140 (MQ=255)
                                                                                      gACGAATATTCTGGCACTGAATGCTGCCGTTGAAGCCGCGCGAGCGGGTGAGCAAGGGCGTGGATTTGCCGTTGTCGCCAGCGAAGTACGGACACTCGCAAGTCGCAGCGCTCAGGCGGCGAAAGAGATTGAAGGCTTGa                                                   >  2:429126/1‑140 (MQ=255)
                                                                                              ctCTGGCACTGACTGCTGCCGTTGAAGCCGCGCGAGCGGGTGAGCAAGGGCGTGGATTTGCCGTTGTCGCCAGCGAAGTACGGACACTCGCAAGTCGCAGCGCTCAGGCGGCGAAAGAGATTGAAGGCTTGATCAGTGaa                                           <  1:1222670/139‑1 (MQ=255)
                                                                                                       tGAATGCTGCCGTTGAAGCCGCGCGAGCGGGTGAGCAAGGGCGTGGATTTGCCGTTGTCGCCAGCGAAGTACGGACACTCGCAAGTCGCAGCGCTCAGGCGGCGAAAGAGATTGAAGGCTTGATCa                                                <  1:960692/126‑1 (MQ=255)
                                                                                                       tGAATGCTGCCGTTGAAGCCGCGCGAGCGGGTGAGCAAGGGCGTGGATTTGCCGTTGTCGCCAGCGAAGTACGGACACTCGCAAGTCGCAGCGCTCAGGCGGCGAAAGAGATTGAAGGCTTGATCa                                                >  2:960692/1‑126 (MQ=255)
                                                                                                            gctgcCGTTGAAGCCGCGCGAGCGGGTGAGCAAGGGCGTGGATTTGCCGTTGTCGCCAGCGAAGTACGGACACTCGCAAGTCGCAGCGCTCAGGCGGCGAAAGAGATTGAAGGCTTGATCAGTGAATCAGTCAGGTTAAt                             <  2:861020/140‑1 (MQ=255)
                                                                                                             ctgcCGTTGAAGCCGCGCGAGCGGGTGAGCAAGGGCGTGGATTTGCCGTTGTCGCCAGCGAAGTACGGACACTCGCAAGTCGCAGCGCTCAGGCGGCGAAAGAGATTGAAGGCTTGATCAGTGAATCAGTCAGGTTAAtt                            >  2:1121144/1‑140 (MQ=255)
                                                                                                             ctgcCGTTGAAGCCGCGCGAGCGGGTGAGCAAGGGCGTGGATTTGCCGTTGTCGCCAGCGAAGTACGGACACTCGCAAGTCGCAGCGCTCAGGCGGCGAAAGAGATTGAAGGCTTGATCAGTGAATCAGTCAGGTTAAtt                            >  1:290125/1‑140 (MQ=255)
                                                                                                                ccGTTGAAGCCGCGCGAGCGGGTGAGCAAGGGCGTGGATTTGCCGTTGTCGCCAGCGAAGTACGGACACTCGCAAGTCGCAGCGCTCAGGCGGCGAAAGAGATTGAAGGCTTGATCAGTGAATCAGTCAGGTTAATTGAc                         >  2:1125562/1‑140 (MQ=255)
                                                                                                                ccGTTGAAGCCGCGCGAGCGGGTGAGCAAGGGCGTGGATTTGCCGTTGTCGCCAGCGAAGTACGGACACTCGCAAGTCGCAGCGCTCAGGCGGCGAAAGAGATTGAAGGCTTGATCAGTGAATCAGTCAGGTTAATTGAc                         <  1:905792/140‑1 (MQ=255)
                                                                                                                               gAGCGGGTGAGCAAGGGCGTGGATTTGCCGTTGTCGCCAGCGAAGTACGGACACTCGCAAGTCGCAGCGCTCAGGCGGCGAAAGAGATTGAAGGCTTGATCAGTGAATCAGTCAGGTTAATTGACCTGGGGTCGGATGAg          <  1:48764/140‑1 (MQ=255)
                                                                                                                                       gAGCAAGGGCGTGGATTTGCCGTTGTCGCCAGCGAAGTACGGACACTCGCAAGTCGCAGCGCTCAGGCGGCGAAAGAGATTGAAGGCTTGATCAGTGAATCAGTCAGGTTAATTGACCTGGGGTCGGATGAGGTGGCAAc  >  2:931442/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                         
GGTGTAGTAAAAACGATGGGCGCTATCTCCACGAGTTCGAAGAAAATTTCTGAGATCACCGCCGTCATCAACAGTATTGCTTTCCAGACGAATATTCTGGCACTGAATGCTGCCGTTGAAGCCGCGCGAGCGGGTGAGCAAGGGCGTGGATTTGCCGTTGTCGCCAGCGAAGTACGGACACTCGCAAGTCGCAGCGCTCAGGCGGCGAAAGAGATTGAAGGCTTGATCAGTGAATCAGTCAGGTTAATTGACCTGGGGTCGGATGAGGTGGCAAC  >  NC_000913/1493529‑1493803

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑