Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 42,108 G→A intergenic (‑177/‑295) caiT ← / → fixA L‑carnitine:gamma‑butyrobetaine antiporter/putative electron transfer flavoprotein FixA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091342,1080GA96.2% 81.0 / NA 26intergenic (‑177/‑295)caiT/fixAL‑carnitine:gamma‑butyrobetaine antiporter/putative electron transfer flavoprotein FixA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (9/16);  minor base T (0/1);  total (9/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACGATTATTAACGGCGGAAAGAAAACCTTCGGTTCTATTCCCGTTTTTCTCTTTTCATTCTTCATGAGTTAATTCCACTGTGAAAACGAATATTTATTTTGCGTTCCCGTTTGTTTTATTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAGGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATCCATAAAACAATATTGAAAATTTCTTTTTGCTACGCCGTGTTTTCAATATTGGTGAGGAACTTAACAATATTGAAAGTTGGATTTATCTGCGTGTGACATTTTCAATATTGGTGATTAAAGTTTTATTTCAAAATTAAAGGGCGTGATATC  >  NC_000913/41866‑42347
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                               
gacgaTTATTAACGGCGGAAAGAAAACCTTCGGTTCTATTCCCGTTTTTCTCTTTTCATTCGTCATGAGTTAATTCCACTGTGAAAATGAATATTTATTTTGCGTTCCCGTTTGTTTTATTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGtttt                                                                                                                                                                                                                                         <  2:28602/251‑1 (MQ=255)
                     gaaaACATTCGGTTCTATTCCCGTTTTTCTCTTTTCATTCTTCATGAGTTAATTCCACTGTGAAAACGAATATTTATTTTGCGTTCCCGTTTGTTTTCTTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATtaata                                                                                                                                                                                                                     >  2:548842/1‑250 (MQ=255)
                      aaaaCCTTCGGTTCTATTCCCGATTTTCTCTTTTCATTCTTCATGAGTTAATTCCACTGTTAAAACGAATATTTATTTTGCGTTCCCGTTTGTTTTATTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGGGTAGTCAATAAGCGGTAGTCTGTGATTGGGATCACATTTTTGTTGCGGGGGAATAGAAGGCGTTTTTTCGTTAATCTTGATTaataat                                                                                                                                                                                                                   <  2:56835/251‑1 (MQ=255)
                        aaCCTTCGGTTCTATTCCCGTTTTTCTCTTTTCATTCTTCATGAGTTAATTCCACTGTGAAAACGAATATTTATTTTGCGTTCCCGTTTGTTTTATTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAAGTTTGATTAATAATCa                                                                                                                                                                                                                 <  2:82251/251‑1 (MQ=255)
                              cGGTTCTATTCCCGTTTTTCTCTTTTCATTCTTCATGAGTTAATTCCACTGTGAAAACGAATATTTATTTTGCGTTCCCGTTTGTTTTATTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGTGTATTCAAGAAGCTGTATTCTGTGATTGGTATCACATTTTTGGTTCGTGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTg                                                                                                                                                                                                            <  1:359478/250‑1 (MQ=255)
                                       tCCCGTTTTTCTCTTTTCATTCTTCATGAGTTAATTCCACTGTGAAAACGAAGATTTATTTTGCGTTCCCGTTTGTTTTATTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGctct                                                                                                                                                                                                   <  1:534326/250‑1 (MQ=255)
                                                   ctTTTCATTCTTCATGAGTTAATTCCACTGTGAAAACGAATATTTATTTTGCGTTCCCGTTTGTTTTATTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTaa                                                                                                                                                                                        <  1:160761/249‑1 (MQ=255)
                                                                  gAGTTAATTCCACTGTGAAAACGAATATTTATTTTGCGTTCCCGTTTGTTTTATTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACtt                                                                                                                                                                       >  1:406308/1‑251 (MQ=255)
                                                                     ttAATTCCACTGTGAAAACGAATATTTATTTTGCGTTCCCGTTTGTTTTATTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTc                                                                                                                                                                     >  1:527770/1‑250 (MQ=255)
                                                                                gtgtaaaCGAATATTTATTTTGCGTTCCCGTTTGTTTTATTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGCTTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGAt                                                                                                                                                         >  1:288809/5‑251 (MQ=255)
                                                                                 tgAAAACGAATATTAATTTTGCGTTCCCGTTTGTTTTATTTTTGTTCACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAAATCATAATTTGTATTCAATAAGCTGTATTATGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGGTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTATAAAATAACATCTGACTTTCAATATTGGTgctc                                                                                                                                                        >  2:372674/1‑248 (MQ=255)
                                                                                       cGAATATTTATTTTGCGTTCCCGTTTGTTTTATTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGGAAGCGTTTATGGCTAACTCATAATGGGTATTCAATTAGCTGTATTCTGTGATTGGTATCCCCTTTTTGATTCGGGTAAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCACTGTTGTGAGTAACAAATAACATCTGACTTTCACTATTGGCGATCCATac                                                                                                                                                   >  2:362044/1‑249 (MQ=255)
                                                                                            atttatttTGCGTTCCCGTTTGTTTTATTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATCCATAAAACAat                                                                                                                                             <  2:104420/251‑1 (MQ=255)
                                                                                              ttatttTGCGTTCCCGTTTGTTTTATTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTTAATAATAATCAGTTTGTTATGCTCTTTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATCCATAAAACAata                                                                                                                                            <  2:59339/250‑1 (MQ=255)
                                                                                                       gTTCCCGTTTGTTTTATTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATCCATAAAACAATATTGAAAAttt                                                                                                                                  <  1:481126/251‑1 (MQ=255)
                                                                                                          cccGTTTGTTTTATTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATCCATAAAACAATATTGAAAATTTCtt                                                                                                                               <  1:359474/251‑1 (MQ=255)
                                                                                                            cGTTTGTTTTATTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATCCATAAAACAATATTGAAAATTTCtttt                                                                                                                             <  1:509061/251‑1 (MQ=255)
                                                                                                                            gTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATCCATAAAACAATATTGAAAATTTCTTTTTGCTACGCCGTGtttt                                                                                                             <  1:20957/251‑1 (MQ=255)
                                                                                                                                          taATTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATCCATAAAACAATATTGAAAATTTCTTTTTGCTACGCCGTGTTTTCAATATTGGTGAgg                                                                                               >  1:38788/1‑251 (MQ=255)
                                                                                                                                                                                 cATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAAGAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATAGTGGTGATCCATAAAACAATATTGAAAATTTCTTTTTGCTACGCCGTGTTTTCAATATTGGTGAGGAACTTAACAATAGTGAAAGTTGGATTTATCTGCGtgtg                                                         <  2:518556/250‑1 (MQ=255)
                                                                                                                                                                                        ggTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATCCATAAAACAATATTGAAAATTTCTTTTTGCTACGCCGTGTTTTCAATATTGGTGAGGAACTTAACAATATTGAAAGTTGGATTTATCTGCGTGTGACATTTTc                                                 <  1:56841/251‑1 (MQ=255)
                                                                                                                                                                                                       gTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATCCATAAAACAATATTGAAAATTTCTTTTTGCTACGCCGTGTTTTCAATATTGGTGAGGAACTTAACAATATTGAAAGTTGGATTTATCTGCGTGTGACATTTTCAATATTGGTGAtt                                    <  1:237411/249‑1 (MQ=255)
                                                                                                                                                                                                         aTTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATCCATAAAACAATATTGAAAATTTCTTTTTGCTACGCCGTGTTTTCAATATTGGTGAGGAACTTAACAATATTGAAAGTTGGATTTATCTGCGTGTGACATTTTCAATATTGGTGATTaa                                  >  2:421418/1‑249 (MQ=255)
                                                                                                                                                                                                                    ggtttctcaTTTTTGTTTCAGGTGAATAGATGGCGTATTTTCTTTAAATTTGATTAATAATCAGTTTGTTATGCTCTGTTGTTAGTAAAAAATAACAGCTGACTTTCAATATTTGTGATCTATAAAATAATATTGAAAATTTCTCTTTGATACGTCGTGTTTTCAATATTTGTGAGGAACTTAACAATATTGAAAGTTGGATTTATCTGCGTGTGACACTTTCATGATGGGTGATTAAAGCTTTATTTCa                      <  2:325577/243‑1 (MQ=255)
                                                                                                                                                                                                                         cacaTTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATCCATAAAACAATATTGAAAATTTCTTTTTGCTACGCCGTGTTTTCAATATTGGTGAGGAACTTAACAATATTGAAAGTTGGATTTATCTGCGTGTGACATTTTCAATATTGGTGATTAAAGTTTTATTTCAAAAtt                 <  2:188104/250‑1 (MQ=255)
                                                                                                                                                                                                                                       gggTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATCCATAAAACAATATTGAAAATTTCTTTTTGCTACGCCGTGTTTTCAATATTGGTGAGGAACTTAACAATATTGAAAGTTGGATTTATCTGCGTGTGACATTTTCAATATTGGTGATTAAAGTTTTATTTCAAAATTAAAGGGCGTGATATc  >  2:313796/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                               
GACGATTATTAACGGCGGAAAGAAAACCTTCGGTTCTATTCCCGTTTTTCTCTTTTCATTCTTCATGAGTTAATTCCACTGTGAAAACGAATATTTATTTTGCGTTCCCGTTTGTTTTATTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAGGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATCCATAAAACAATATTGAAAATTTCTTTTTGCTACGCCGTGTTTTCAATATTGGTGAGGAACTTAACAATATTGAAAGTTGGATTTATCTGCGTGTGACATTTTCAATATTGGTGATTAAAGTTTTATTTCAAAATTAAAGGGCGTGATATC  >  NC_000913/41866‑42347

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