Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,931,183 (G)7→8 intergenic (‑90/‑133) ravA ← / → kup regulatory ATPase RavA/K(+):H(+) symporter Kup

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,931,1761.G96.7% 90.4 / ‑3.0 30intergenic (‑83/‑140)ravA/kupregulatory ATPase RavA/K(+):H(+) symporter Kup
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (15/14);  total (15/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.18e-01

AATTCTTTCCGCTAATAAATGAGGGTGAGCCATAATGAAGTGGCGTCCTTTCGTCAAAAGTTCTGCGTAAATTGCGAGTATAGACGTTTCTTGCTGGTGGCTAAAATAGTCTCAAA‑‑GGGGGGGTATTTTTCTTTGAGCCAGGTTAATGTGGCCGCATTTAGGAGTACGATTTTGCCGTTAATCGTGCATACTGTGCGCTTTTTTGTGGGCCAAGGGACTAAGCACACATTTCATATTTCAACG  >  NC_000913/3931061‑3931303
                                                                                                                    |                                                                                                                                
aaTTCTTTACGCTAATAATTGATGGTGAGCCATAATGAAGTGGCGTCCTTTCGTCAAAAGTTCTGCGTAAATTGCGAGTATAGACGTTTCTTGCTGGTGGCTAAAATAGTCTCAAAGGGGGGGGTATTTTTCt                                                                                                                 <  1:134615/133‑1 (MQ=255)
 aTTCTTTCCGCTAATAAATGAGGGTGAGCCATAATGAAGTGGCGTCCTTTCGTCAAAAGTTCTGCGTAAATTGCGAGTATAGACGTTTCTTGCTGGTGGCTAAAATAGTCTCAAAGGGGGGGGTATTTTTCtt                                                                                                                <  1:498930/133‑1 (MQ=255)
     tttCCGCTAATAAATGAGGGTGAGCCATAATGAAGTGGCGTCCTTTCGTCAAAAGTTCTGCGTAAATTGCGAGTATAGACGTTTCTTGCTGGTGGCTAAAATAGTCTCAAAGGGGGGGGTATTTTTCTTTGAg                                                                                                            <  2:7279/133‑1 (MQ=255)
                        gTGAGCCATAATGAAGTGGCGTCCTTTCGTCAAAAGTTCTGCGTAAATTGCGAGTATAGACGTTTCTTGCTGGTGGCTAAAATAGTCTCAAAGGGGGGGGTATTTTTCTTTGAGCCAGGTTAATGTgg                                                                                              >  1:128636/1‑128 (MQ=255)
                        gTGAGCCATAATGAAGTGGCGTCCTTTCGTCAAAAGTTCTGCGTAAATTGCGAGTATAGACGTTTCTTGCTGGTGGCTAAAATAGTCTCAAAGGGGGGGGTATTTTTCTTTGAGCCAGGTTAATGTgg                                                                                              <  2:128636/128‑1 (MQ=255)
                        gTGAGCCATAATGAAGTGGCGTCCTTTCGTCAAAAGTTCTGCGTAAATTGCGAGTATAGACGTTTCTTGCTGGTGGCTAAAATAGTCTCAAAGGGGGGGGTATTTTTCTTTGAGCCAGGTTAATGTGGCCGCa                                                                                         >  2:305796/1‑133 (MQ=255)
                           aGCCATAATGAAGTGGCGTCCTTTCGTCAAAAGTTCTGCGTAAATTGCGAGTATAGACGTTTCTTGCTGGTGGCTAAAATAGTCTCAAAGGGGGGGGTATTTTTCTTTGAGCCAGGTTAATGTGGCCGCAttt                                                                                      >  2:25247/1‑133 (MQ=255)
                           aGCCATAATGAAGTGGCGTCCTTTCGTCAAAAGTTCTGCGTAAATTGCGAGTATAGACGTTTCTTGCTGGTGGCTAAAATAGTCTCAAAGGGGGGGGTATTTGTCTTTGAGCCAGGTTAATGTGCCCGCAttt                                                                                      >  1:167263/1‑133 (MQ=255)
                                 aaTGAAGTGGCGTCCTTTCGTCAAAAGTTCTGCGTAAATTGCGAGTATAGACGTTTCTTGCTGGTGGCTAAAATAGTCTCAAAGGGGGGGGTATTTTTCTTTGAGCCAGGTTAATGTGGCCGCATTTAGGAGt                                                                                >  2:77700/1‑133 (MQ=255)
                                     aaGTGGCGTCCTTTCGTCAAAAGTTCTGCGGAAATTGCGAGTATAGACGTTTCTTGCTGGTGGCTAAAATAGTCTCAAAGGGGGGGGTATTTTTCTTTGAGCCAGGTTAATGTGGCCGCATTTAGGAGTACGa                                                                            >  2:95366/1‑133 (MQ=255)
                                         ggCGTCCTTTCGTCAAAAGTTCTGCGTAAATTGCGAGTATAGACGTTTCTTGCTGGTGGCTAAAATAGTCTCAAAGGGGGGGGTATTTTTCTTTGAGCCAGGTTAATGTGGCCGCATTTAGGAGTACGAtttt                                                                        <  1:77700/133‑1 (MQ=255)
                                           cGTCCTTTCGTCAAAAGTTCTGCGTAAATTGCGAGTATAGACGTTTCTTGCTGGTGGCTAAAATAGTCTCAAAGGGGGGGGTATTTTTCTTTGAGCCAGGTTAATGTGGCCGCATTTAGGAGTACGATTTTGc                                                                      >  1:461903/1‑133 (MQ=255)
                                           cGTCCTTTCGTCAAAAGTTCTGCGTAAATTGCGAGTATAGACGTTTCTTGCTGGTGGCTAAAATAGTCTCAAAGGGGGGGGTATTTTTCTTTGAGCCAGGCTAATGTGTCCGCATTTAGGAGTACGATTTTGc                                                                      >  2:217692/1‑133 (MQ=255)
                                             tCCTTTCGTCAAAAGTTCTGCGTAAATTGCGAGTATAGACGTTTCTTGCTGGTGGCTAAAATAGTCTCAAAGGGGGGGGTATTTTTCTTTGAGCCAGGTTAATGTGGCCGCATTTAGGAGTACGATTTTGCCg                                                                    >  2:112690/1‑133 (MQ=255)
                                              ccTTTCGTCAAAAGTTCTGCGTAAATTGCGAGTATAGACGTTTCTTGCTGGTGGCTAAAATAGTCTCAAAGGGGGGGGTATTTTTCTTTGAGCCAGGTTAATGTGGCCGCATTTAGGAGTACGATTTTGCCGt                                                                   <  1:305796/133‑1 (MQ=255)
                                                tttCGTCAAAAGTTCTGCGTAAATTGCGAGTATAGACGTTTCTTGCTGGTGGCTAAAATAGTCTCAAAGGGGGGGGTATTTTTCTTTGAGCCAGGTTAATGTGGCCGCATTTAGGAGTACGATTTTGCCGTTa                                                                 <  2:458500/133‑1 (MQ=255)
                                                tttCGTCAAAAGTTCTGCGTAAATTGCGAGTATAGACGTTTCTTGCTGGTGGCTAAAATAGTCTCAAAGGGGGGGGTATTTTTCTTTGAGCCAGGTTAATGTGGCCGCATTTAGGAGTACGATTTTGCCGTTa                                                                 <  2:469740/133‑1 (MQ=255)
                                                 ttCGTCAAAACTTCTACGTAAATGACGGATATCGCCGTTTCTTGCTGGTGGCTAAAATAGTCTCAAAGGGGGGGGTTTTTTTCTTTGAGCCAGGTTAATGTGGCCGCATTTAGAGTACGATTTTGCCGTTAAt                                                               <  2:141202/133‑1 (MQ=255)
                                                              cTGCGTAAATTGCGAGTATAGACGTTTCTTGCTGGTGGCTAAAATAGTCTCAAAGGGGGGGGTATTTTTCTTTGAGCCAGGTTAATGTGGCCGCATTTAGGAGTACGATTTTGCAGTGAATCGTACATACtgt                                                   >  1:238879/1‑133 (MQ=255)
                                                                      aTTGCGAGTATAGACGTTTCTTGCTGGTGGCTAAAATAGTCTCAAAGGGGGGGGTATTTTTCTTTGAGCCAGGTTAATGTGGCCGCATTTAGGAGTACGATTTTGCCGTTAATCGTGCATACTGTGCGCtttt                                           >  1:25453/1‑133 (MQ=255)
                                                                           gAGTATAGAGCTTTCCTGCTGGTTTCGAGAATAGTCTCAAAGGGGGGGGTATTTTTCTTTGAGCCAGGTTAATGTGGCCGCATTTAGGAGTACGATTTTGCCGTTAATCGTGCATACTGTGCGCTTTTTtatg                                      <  2:238879/133‑1 (MQ=255)
                                                                            aGTATAGACGTTTCTTGGTGGTGGCTAAAATAGTCTCAAAGGGGGGGGGTATTTTTCTTTGAGCCAGGTTAATGTGGCCGCATTTAGGAGTACGATTTTGCCATTAATCGGGCATACGGTGCGCCTTTTtgtg                                      >  1:234461/1‑133 (MQ=255)
                                                                               ataGACGTTTCTTGCTGGTGGCTAAAATAGTCTCAAA‑‑GGGGGGGTATTTTTCTTTGAGCCAGGTTAATGTGGCCGCATTTAGGAGTACGATTTTGCCGTTAATCGTGCATACTGTGCGCTTTTTTGTGGGCCa                                 <  1:404502/133‑1 (MQ=255)
                                                                                taGACGTTTCTTGCTGGTGGCTAAAATAGTCTCAAAGGGGGGGGTATTTTTCTTTGAGCCAGGTTAATGTGGCCGCATTTAGGAGTACGATTTTGCCGTTAATCGTGCATACTGTGCGCTTTTTTGTGGGCCa                                 >  2:94983/1‑133 (MQ=255)
                                                                                     gTTTCTTGCTGGTGGCTAAAATAGTCTCAAAGGGGGGGGTATTTTTCTTTGAGCCAGGTTAATGTGGCCGCATTTAGGAGTACGATTTTGCCGTTAATCGTGCATACTGTGCGCTTTTTTGTGGGCCAAGGGa                            <  1:25247/133‑1 (MQ=255)
                                                                                            gCTGGTGGCTAAAATAGTCTCAAAGGGGGGGGTATTTTTCTTTGAGCCAGGTTAATGTGGCCGCATTTAGGAGTACGATTTTGCCGTTAATCGTGCATACTGTGCGCTTTTTTGTGGGCCAAGGGACTAAGca                     <  1:112690/133‑1 (MQ=255)
                                                                                                          tAGTCTCAAAGGGGGGGGTATTTTTCTTTGAGCCAGGTTAATGTGGCCGCATTTAGGAGTACGATTTTGCCGTTAATCGTGCATACTGTGCGCTTTTTTGTGGGCCAAGGGACTAAGCACACATTTCATAttt       <  1:312194/133‑1 (MQ=255)
                                                                                                             tctcAAAGGGGGGGGTATTTTTCTTTGAGCCAGGTTAATGTGGCCGCATTTAGGAGTACGATTTTGCCGTTAATCGTGCATACTGTGCGCTTTTTTGTGGGCCAAGGGACTAAGGACGCATTTCTTATTTCaa    >  1:201875/1‑133 (MQ=255)
                                                                                                             tctcAAAGGGGGGGGCATTTTTCTTTGAGCCAGGTTAATGTGGCCGCATTTAGGAGTACGATTTTGCCGTTAATCGTGCATACTGTGCGCTTTTTTGTGGGCCAAGGGACTAAGCACACATTTCATATTTCaa    >  1:19397/1‑133 (MQ=255)
                                                                                                               tcAAAGGGGGGGGTATTTTTCTTTGAGCCAGGTTAATGTGGCCGCATTTAGGAGTACGATTTTGCCGTTAATCGTGCATACTGTGCGCTTTTTTGTGGGCCAAGGGACTAAGCACACATTTCATATTTCAACg  <  1:79850/133‑1 (MQ=255)
                                                                                                                    |                                                                                                                                
AATTCTTTCCGCTAATAAATGAGGGTGAGCCATAATGAAGTGGCGTCCTTTCGTCAAAAGTTCTGCGTAAATTGCGAGTATAGACGTTTCTTGCTGGTGGCTAAAATAGTCTCAAA‑‑GGGGGGGTATTTTTCTTTGAGCCAGGTTAATGTGGCCGCATTTAGGAGTACGATTTTGCCGTTAATCGTGCATACTGTGCGCTTTTTTGTGGGCCAAGGGACTAAGCACACATTTCATATTTCAACG  >  NC_000913/3931061‑3931303

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: