Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,004,180 G→A intergenic (‑2/+50) rarD ← / ← yigI putative transporter RarD/putative thioesterase YigI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,004,1800GA96.2% 71.9 / NA 26intergenic (‑2/+50)rarD/yigIputative transporter RarD/putative thioesterase YigI
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (11/14);  minor base C (1/0);  total (12/14)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATA‑TAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCAACGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCTGGCGT  >  NC_000913/4003934‑4004422
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   
aGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCTCCACCATAAAGAAAAACGACCAGAGCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATATAATAAGCGGCAAGAGCGATTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCAGGaa                                                                                                                                                                                                                                                 <  1:45895/251‑1 (MQ=255)
 gCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCAGGAAt                                                                                                                                                                                                                                                <  1:340988/251‑1 (MQ=255)
                               tCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCTGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCAGGAATACTCATTTTTAATTGGTAACAGCAGTTTa                                                                                                                                                                                                                   <  1:259981/250‑1 (MQ=255)
                                           gAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATCCCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTGTGTTTTGCATCCATCAGGAATACTCAGTTTTAATTGGTAACAGCAGTTTACCTGCTTTTAt                                                                                                                                                                                                        <  2:455864/249‑1 (MQ=255)
                                               cAGCGTTTTTAAATAGGACCACTGGCTGCAAATGCTCAGCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCGCGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGTGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCTATCAGGAATACTCATGTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCtt                                                                                                                                                                                                   <  1:176345/250‑1 (MQ=255)
                                                            tAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGGACGTAGAGAAAAAACTTTACGTAAGCCGGCGCTATACACACAATAAAAAAAGCGGCAAAAGCGGGGAAAAAGCCCTGCCGCGTTTGTTTTTCTTCCATCAGGGATAATCATTTTTTATTGGTAAAAACAGTTTACCTGCCTTTTTGTCtttaaacctacatatt                                                                                                                                                                                     >  1:427512/1‑237 (MQ=255)
                                                             agaaCCACTGGCGACAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTTCATCCATCAGGAATACTCATTGTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCATCCTACCAta‑ta                                                                                                                                                                                    <  2:413330/247‑1 (MQ=255)
                                                                                    atcaGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCAGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATA‑TAGGTGGCGGTGGCACTGGCAAtata                                                                                                                                                            <  1:181510/251‑1 (MQ=255)
                                                                                                aTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCAAGGAGAAAATAAAATTGAAAGAAGGTGGGGGTAAAACCCCAATAAAATAAAAGGCAAAAGCGAGGAAAAAGCCCCGCCACGGTTGTTTTTCAGTCAACATGGATAATCATTTTTAATTGGTAACAACAGTTTACCTTCTTTTTTGTCTTCAACTTACCCAA‑TAGGGGGGGGGGGGAATGGCAatattaaaatgtttttt                                                                                                                                                >  1:85502/1‑238 (MQ=255)
                                                                                                 tAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGGAAGCTGGCGCTATAACCCAAATAAAATAAACGGCCAGAGCGGGTAATACGCCCTGCCGCGGTTGTTTTGCATCCACCCGGGATACTCATTTTTAATTGGTAACAGCAGTTTAACTGCTTTTTTGTTTTCAACCTACCATA‑TAGGGGGCGGGGGCACTGGCCATATAAAGCGGtttttct                                                                                                                                               >  1:133718/1‑250 (MQ=255)
                                                                                                        aaaCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCAAGAACTAAATCAACTTTAAGTAAACTGGCGCTATAACCCCAATAAAATAACCCGCAAGAGAGAGGAATACGCCCCTCCCAGTTTGTTTTGCCTCCATCAGGAATACTCATTTTTTATTTGGAACAGCAGTTTACCTGCTTTTTTTTTTTCAACCTACCATAATAG‑TGGCGGGGGGCCTGGCAAAATAAAACTTTTCTTCATTttgt                                                                                                                                         >  1:265604/1‑250 (MQ=255)
                                                                                                                  aTCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCAGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCCTA‑TAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAAcgcg                                                                                                                              <  1:89913/251‑1 (MQ=255)
                                                                                                                                                   aCGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTACCGCGTTTGTTTTGCATCCATCAGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATA‑TAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCACAGCGGGCGAAGGCGGCTTTATTGCCTGCACGCAAc                                                                                               >  2:146667/1‑249 (MQ=255)
                                                                                                                                                     tgcgtaAATCAACTTGAAGTACGCTGGCGATATACCCCAAATAAAATAAGCGGCAAGAGCGAGCAATACGCACTACCGCGATTGTTTTGCATCCAACAGAATTACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCACTCTACCATA‑TAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTATTCATTGTGTAATTCAACGCGTGCGACGCCAACTTTATTGcttgaccgtaaaaggc                                                                                           >  2:93641/4‑236 (MQ=255)
                                                                                                                                                      tagtaAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCAGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATA‑TAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGct                                                                                          >  1:257143/1‑251 (MQ=255)
                                                                                                                                                          aaaTCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCAGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATA‑TAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGctacta                                                                                      <  1:166645/251‑1 (MQ=255)
                                                                                                                                                                            cTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGAAAGTAATACGCCCTGCCGCGTTTGTTTTACATCCATCAGGAATACTCATTTTTACTTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATA‑TAGGTGGCGGTGGCACTGCCAATATAAAGCAGTTCTTCATTGTGTAATTCAACGCGGGCGACGGAGACTTTATTGCCTGCACCCAACAGGCTACTAGTAGCAGTAAAACGATCg                                                                    >  1:261251/1‑251 (MQ=255)
                                                                                                                                                                                            aataaaataaGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGATTTGCATCCATCAGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATA‑TAGGCGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGgcgc                                                     >  1:260375/1‑250 (MQ=255)
                                                                                                                                                                                             ataaaataaGCGGCAAGAGCGAGGAATTCGCCCTGCCGCGTTTGTTTTGCATCCATCAGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATA‑TAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAg                                                   <  2:267485/251‑1 (MQ=255)
                                                                                                                                                                                                   taaGCGACAAGAGCGAGTAATAAGCCCTGCCGCGTTTGTTTTGCATCCATCAGGCATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATA‑TAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATc                                             >  2:136803/1‑251 (MQ=255)
                                                                                                                                                                                                         gCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCAGGAATACTCATTTTTAATTGGTAACAGCAGTTTCCCTGCTTTTATGTCTTCAACCTACCATA‑TAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCAAcgcg                                       >  1:401278/1‑251 (MQ=255)
                                                                                                                                                                                                                cGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCAGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATA‑TAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCAACGCGAAGATCa                                <  1:12033/251‑1 (MQ=255)
                                                                                                                                                                                                                                ccGCGTTTGTTTTGCATCCATCAGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATA‑TAGGTGGCGGGGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGgtcgacggc                                                                                                                     <  2:326662/150‑1 (MQ=255)
                                                                                                                                                                                                                                     tttgtttTGCATCCCACAGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATA‑TAGGTGGCGGTGGCACTGTCAATATAAAGCTGTTCTTCATTGGGTAATTCATCGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCGCGCCCCTGCCTGGGCGCAGATATTCAACGCGAAGATCAATGGTCCCCATccacgttag            <  2:7168/250‑8 (MQ=255)
                                                                                                                                                                                                                                        gtttTTCATCCATCAGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGGTTTTATGTCTTCAACCTACCATA‑TAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCAACGCGAAGAGTAATGGTCCCCACCCGCGATAGCCg         <  2:234465/250‑1 (MQ=255)
                                                                                                                                                                                                                                                tccatcAGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATA‑TAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTACTTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCAACACGAAGATCACTGGTCCCCATCCGCGCTAGCCGCTGGCGt  >  2:244769/1‑249 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   
AGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATA‑TAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCAACGCGAAGATCAATGGTCCCCATCCGCGATAGCCGCTGGCGT  >  NC_000913/4003934‑4004422

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