Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,245,036 A→C 21.3% I10S (ATC→AGC)  aruH ← arginine‑‑pyruvate transaminase AruH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,245,0360AC21.3% 56.9 / 15.4 37I10S (ATC→AGC) aruHarginine‑‑pyruvate transaminase AruH
Reads supporting (aligned to +/- strand):  ref base A (15/14);  new base C (5/3);  total (20/17)
Fisher's exact test for biased strand distribution p-value = 7.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.14e-01

CGCGCAGGCTGTCGATGGCCGCCTCGACGATAGGCGCGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAACGCGCCGTTGATTTCTTCCATCTCGGCCTGGGTG  >  NC_002947/4244887‑4245182
                                                                                                                                                      |                                                                                                                                                  
cgcgcAGGCTGTCGATGGCCGCCTCGACGATAGGCCCGGCAGTAGCGAAGGCCTGGTCACCCACCGACAGCAGCATGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGa                                                                                                                                                    <  1:444794/150‑1 (MQ=255)
   gcAGGCTGTCGATGGCCGCCTCGACGATAGGCGCGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGc                                                                                                                                                 >  2:477636/1‑150 (MQ=255)
     aGGCTGTCGATGGCCGCCGCGGCGATAGGCGCGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGt                                                                                                                                               <  2:123190/150‑1 (MQ=255)
        cTGTCGATGGCCGCCTCGACGATAGGCGCGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTgg                                                                                                                                            <  1:477636/150‑1 (MQ=255)
            cGATGGCCGCCTCGACGATAGGCGCGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGa                                                                                                                                        <  1:400110/150‑1 (MQ=255)
                      cTCGACGATAGGCGCGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTa                                                                                                                              <  1:393333/150‑1 (MQ=255)
                             aTAGGCGCGGGAGTGGCGAAGTCCGGATCACCCACCGACAGCAGCAGTACGTCCTTGTCCGCAGCCTGCAGCGCCAGGGAGCGGTAGTGGACGTTCCCGGCCGGCTAGTTGTCCTGTAGGATGCGTTGGGTGAGCTTGGCGTAGCGCATg                                                                                                                       <  1:108502/150‑1 (MQ=255)
                                     gggAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAAtctc                                                                                                               <  1:225545/150‑1 (MQ=255)
                                          tgtCGAAGTCCGGGTCACCCACCGACCACAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTc                                                                                                          >  1:267654/1‑150 (MQ=255)
                                                      ggTCACCCACCGACAGCAGCAGGATGTCCTTGACCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATAGCGTACGACGCCAGCCAGGCCACGATGAGGCTTTGGGTGAGCTTGGCGTAGCGCATGGCACTCTAGTCTCAACGTCAATTTc                                                                                              >  1:200826/1‑150 (MQ=255)
                                                         cacccaccGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGt                                                                                           <  2:393533/150‑1 (MQ=255)
                                                          acccaccGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTc                                                                                          >  1:249900/1‑150 (MQ=255)
                                                            ccaccGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCATCGCCGGCGACGGCGCTg                                                                                                                                                  <  1:541238/92‑1 (MQ=255)
                                                            ccaccGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGGGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTcga                                                                                        >  2:377982/1‑150 (MQ=255)
                                                             caccGACAGCAGCAGGCTGCCCTTGCCCTCAGCCTGCAGAACAAGGGCGAGGCACTGGCTGTCCCCGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTAAGCTTGGCGCAGCGCATGGCAATCTCCTTTCAATGACAATTTCAGTcgac                                                                                       >  2:178831/1‑150 (MQ=255)
                                                               ccGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTcgacga                                                                                     <  1:540259/150‑1 (MQ=255)
                                                                        gcagGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAAt                                                                                                 <  1:437007/129‑1 (MQ=255)
                                                                        gcagGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAAt                                                                                                 >  2:437007/1‑129 (MQ=255)
                                                                        gcagGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCg                                                                            <  2:156633/150‑1 (MQ=255)
                                                                          agGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCggg                                                                          <  2:200826/150‑1 (MQ=255)
                                                                                         tCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCAt                                                           >  1:281971/1‑150 (MQ=255)
                                                                                                  aGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTAc                                                  >  2:405001/1‑150 (MQ=255)
                                                                                                                   gTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAAcgc                                 <  2:322486/150‑1 (MQ=255)
                                                                                                                           cccAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAGTTTCAGTAAACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAACGCGCCGTTGa                         >  1:349797/1‑150 (MQ=255)
                                                                                                                           cccAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAACGCGCCg                             <  2:383046/146‑1 (MQ=255)
                                                                                                                           cccAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTGCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAACGCGCCg                             >  1:383046/1‑146 (MQ=255)
                                                                                                                            acccaacGCATCGCCGGCGACGGCGCTGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAACGCGCCGTTGAt                        >  1:162079/7‑150 (MQ=255)
                                                                                                                             cAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAACGCGCCGTTGAtt                       <  2:281971/150‑1 (MQ=255)
                                                                                                                              ccaacGCATCGCCGGCGACGGCGCTGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGgt                                                      <  2:495997/114‑1 (MQ=255)
                                                                                                                              ccaacGCATCGCCGGCGACGGCGCTGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGgt                                                      >  1:495997/5‑118 (MQ=255)
                                                                                                                              ccaacGCATCGCCGGCGACGGCGCTGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAACGCGCCGTTGAttt                      >  1:368123/5‑150 (MQ=255)
                                                                                                                              ccaacGCATCGCCGGCGACGGCGCTGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAACGCGCCGTTGAttt                      >  1:302970/5‑150 (MQ=255)
                                                                                                                                   gCATCGCCGGCGACGGCGCTGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAACGCGCCGTTGATTTCTTcc                 >  1:613148/1‑150 (MQ=255)
                                                                                                                                   gCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAACGCGCCGTTGATTTCTTcc                 >  2:260127/1‑150 (MQ=255)
                                                                                                                                          ccggcgacGGCGCTGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAACGCGCCGTTGATTTCTTCCATCTCg           <  2:613148/143‑1 (MQ=255)
                                                                                                                                          ccGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATTTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAACGCGCCGTTGATTTCTTCCATCTCg           >  1:346827/1‑150 (MQ=255)
                                                                                                                                                 cGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAACGCGCCGTTGATTTCTTCCATCTCGGCCTggg    >  1:332445/1‑150 (MQ=255)
                                                                                                                                                   gCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAACGCGCCGTTGATTTCTTCCATCTCGGCCTGGGTg  >  1:288421/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
CGCGCAGGCTGTCGATGGCCGCCTCGACGATAGGCGCGGGAGTGTCGAAGTCCGGGTCACCCACCGACAGCAGCAGGATGTCCTTGTCCTCAGCCTGCAGCGCCAGGGCGCGGTAGTGGATGTCCCAGGCCGCAGCGCCGTCGCCGGCGATGCGTTGGGTGAGCTTGGCGTAGCGCATGGCAATCTCCTTTCAATGGCAATTTCAGTCGACGAGTATCTTCGGGTTGGAATTGCCCCATACGGTGTACAAATCGGCGAGCAACGCGCCGTTGATTTCTTCCATCTCGGCCTGGGTG  >  NC_002947/4244887‑4245182

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