Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,123,823 2 bp→AG coding (231‑232/648 nt) grxB ← reduced glutaredoxin 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,123,8230GA87.5% 52.2 / ‑0.5 24L78L (CTG→TTG) grxBreduced glutaredoxin 2
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (8/13);  total (10/14)
Fisher's exact test for biased strand distribution p-value = 5.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01
*NC_0009131,123,8240TG90.9% 58.6 / ‑4.4 24L77F (TTA→TTC) grxBreduced glutaredoxin 2
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base G (8/12);  total (10/14)
Fisher's exact test for biased strand distribution p-value = 5.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

CAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAGAGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AGTAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGG‑GTGTTTCTGCGTCGTCGTTGAGCAGAACATGTAATTCGACGGGGATATTTTTCAGGCCGAAAATCATGCGGGCTTTGAGGCAGTAAGGGCAGTGATCGTAAATGTATAGCTTCACGT  >  NC_000913/1123584‑1124056
                                                                                                                                                                                                                                                ||                                                                                                                                                                                                                                         
cAGTTTGCCCAGCGCACGTAAATCATCGCTGATATTTTTAATCAGCCCGTCAGAGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGc                                                                                                                                                                                                                                 <  1:269678/251‑1 (MQ=255)
          cgactctctccctatctgCTGATATGTTTAACCAGACCGTCAGTGTGGGCCAGAAGGTCGGCAAAATTCCCCGCGTTGTCCTCTTTCGGGTCGACGAAATATTTGCGCGGGGCGGGAATAGAAGACCCATCAAATGCCGATTTGGCAAACCGCGGCAACAGCAGGTTGTGGGCGTAGCAATAGACCTGGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTACCGTCGAGttt                                                                                                                                                                                                                       <  2:17249/234‑1 (MQ=255)
                       catcGCTGATATTCTTAATCAGACCGTCAGAGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGa                                                                                                                                                                                                          <  1:95468/251‑1 (MQ=255)
                           gCTGATATTCTTAATCAGACCGTCAGAGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCCTCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGttt                                                                                                                                                                                                                       >  1:17249/1‑234 (MQ=255)
                             tGATATTCTTAATCAGACCGTCAGAGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTTCCGTCGAGTTTATCGACATAGTGAACGATg                                                                                                                                                                                                    >  1:171667/1‑251 (MQ=255)
                                tattcttaatatgtcatgatgggTGGGCCAGAAGGGCGGCAAAATTACCCGCGCTGGCTTTTTTCGTGTCGACAAAATATTTGCGCGCGGCGGGAGTAGAAAACTCAGCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTGGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTc                                                                                                                                                                                                  <  2:64404/228‑1 (MQ=255)
                                    cTTAATCAGACCGTCAGAGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTCCAG‑AACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGc                                                                                                                                                                                             >  1:156213/1‑251 (MQ=255)
                                                       ggggCCAGCAGGTCGGCAAATTTCCCTGCGCTGGCTTCTTTCTTGTCGACGAAATATGTGGGCGGGGGGGGAGGAGAAAACTCATCAAAGGCGGATTGGGAAAAAGGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AGTAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCt                                                                                                                                                                          <  2:143264/250‑1 (MQ=255)
                                                              gcagGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCC‑TCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTc                                                                                                                                                                                         <  1:223081/228‑1 (MQ=255)
                                                              gcagGTCGGCAAAATTACCCGAGCTGGCCTCTTTATTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCAGATTTGGCAAAACGCGGAAACAGAAGTTTGTTGGCGTAACCATTGACCGTGCGACGCCACTCTTCAAAGGCAGGGGGACGTTTGCCGGTCCAG‑AACGGTTTGCC‑CCGAGTTTAAAGACATAGTGAAAGGTGTTCATGCTTTcagatcggaagatggttgggaggg                                                                                                                                                                  >  2:223081/1‑228 (MQ=255)
                                                                   tCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCtttttt                                                                                                                                                              >  1:62123/1‑250 (MQ=255)
                                                                        aaaaTTGCCCGCGCTGTCCTCTTTCTTGTCGACGAAATATTGGCGCGCGGGGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTCGCCATTGACCTTGCGCCGCCACTCTTCAATGGCAGGGGAATGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCGTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGaa                                                                                                                                                         <  1:29108/251‑1 (MQ=255)
                                                                             tACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAACCGCGGCAACAGCAGTTTTTGGGCGTAGCCATTGGCCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATggg                                                                                                                                                     <  1:45442/250‑1 (MQ=255)
                                                                               cccGAGCTGGCCACTTTCTGGCCGACGAATTATTAGCGCGCGGCAGGAGTTGAAAACTCATCAAATGCCGTTTCGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAAc                                                                                                                                                  <  2:117582/251‑1 (MQ=255)
                                                                                                                           gggAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCCTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAAAGGTGTCCATGCTTTCTGGCATATAGCGGCTGTCCTCCTTTTGCTGAATGGGAACCTGTTTTTGACCGACCCTCCGGGTGG‑GTGTTTCTGcgtcgttgt                                                                                                     >  1:247871/1‑248 (MQ=255)
                                                                                                                                                 gCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGG‑GTGTTTCTGCGTCGTCGTTGAGCAGAACATGTAATTCGAc                                                                               <  1:135677/251‑1 (MQ=255)
                                                                                                                                                   cGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGG‑GTGTTTCTGCGTCGTCGTTGAGCAGAACATGTAATTCGaagg                                                                             >  1:30300/1‑248 (MQ=255)
                                                                                                                                                    gATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGTACCTGTTTTTGACAGACCATCCGGGTGG‑GTGTTTCTGCGTCGTCGTTGAGCAGAACATGTAATTCGTCggg                                                                            >  1:19640/1‑251 (MQ=255)
                                                                                                                                                            aaaaCGTGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGG‑GTGTTTCTGCGTCGTCGTTGAGCAGAACATGTAATTCGACGGGGAGAtttt                                                                    >  1:92567/1‑251 (MQ=255)
                                                                                                                                                               aCGCGGAAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGTAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGTGAACCTGTTTTTGACCGACCATCCGGGTGG‑GTGTTTCTGCGTCGTCGTTGAGCAGAACATGTAATTCGACGGGGATATTTTTCa                                                                 <  1:36705/251‑1 (MQ=255)
                                                                                                                                                                                  gtggtgttgctgttgacttcgcgcAGCCATTCTCCAAGGGCATGGGAACTGTTGCCGGTT‑AATAACGTTTGGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGGCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGTCCATCCGGGTGG‑GTGTTTCTGCGTCGTCGTTGAGCAGAACATGTAATGAGACGGGGATATTTTTCAGGCCGAAAATCATGCgtgc                                              <  2:120403/231‑4 (MQ=255)
                                                                                                                                                                                        tAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGG‑GTGTTTCTGCGTCGTCGTTGAGCAGAAAATGTAATTCGACGGGGATATTTTTCAGGCCGAAAATCATGCGGGCTTTg                                          >  1:16742/1‑249 (MQ=255)
                                                                                                                                                                                                      gcgtggcccctctTCAATGGCAGGGGTGCGTGTGCATGTC‑AAGATCGGTATGTCGTCTTGTTTATCGTCTGAGTTAACGATGTCCATGCTTCCTGGCAGATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACTGACACTCCGG‑TGGTGTGTTTCTGCGTCGTTGTTGAGCATATCGTGTAAGTCGACGGTGATAGTTTTTAGGGCGAAAATCAGGTTGGCTTTGAGGCAGTAATGGCAGt                          <  2:184980/242‑1 (MQ=255)
                                                                                                                                                                                                                              ggAACGTTTGCCGGTC‑AAGAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGG‑GTGTTTCTGCGTCGTCGTTGAGCAGAACATGTAATTCGACGGGGATATTTTTCAGGCCGAAAATCATGCGGGCTTTGAGGCAGTAAGGGCAGTGATCGTAAATGTATAGCTTCACGt  <  1:331803/251‑1 (MQ=255)
                                                                                                                                                                                                                                                ||                                                                                                                                                                                                                                         
CAGTTTGTCCAGCGCACGTAAATCATCGCTGATATTCTTAATCAGACCGTCAGAGTGGGCCAGCAGGTCGGCAAAATTACCCGCGCTGGCCTCTTTCTTGTCGACGAAATATTTGCGCGCGGCGGGAGTAGAAAACTCATCAAATGCCGATTTGGCAAAACGCGGCAACAGCAGTTTGTTGGCGTAGCCATTGACCTTGCGCAGCCACTCTTCAATGGCAGGGGAACGTTTGCCGGTC‑AGTAACGGTTTGCCGTCGAGTTTATCGACATAGTGAACGATGTCCATGCTTTCTGGCATATAGCGGCTGTCATCTTTTTGCAGAATGGGAACCTGTTTTTGACCGACCATCCGGGTGG‑GTGTTTCTGCGTCGTCGTTGAGCAGAACATGTAATTCGACGGGGATATTTTTCAGGCCGAAAATCATGCGGGCTTTGAGGCAGTAAGGGCAGTGATCGTAAATGTATAGCTTCACGT  >  NC_000913/1123584‑1124056

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