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,1080GA95.0% 62.0 / ‑5.9 20intergenic (‑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/10);  minor base C (0/1);  total (9/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GACGATTATTAACGGCGGAAAGAAAACCTTCGGTTCTATTCCCGTTTTTCTCTTTTCATTCTTCATGAGTTAATTCCACTGTGAAAACGAATATTTATTTTGCGTTCCCGTTTGTTTTATTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAGGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATCCATAAAACAATATTGAAAATTTCTTTTTGCTACGCCGTGTTTTCAATATTGGTGAGGAACTTAACAATATTGAAAGTTGGATTTATCTGCGTGTGACATTTTCAATATTGGTGATTAAAGTTTTATTTCAAAATTAAAGGGC  >  NC_000913/41866‑42339
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                       
gacgaTTATTAACGGCGGAATGAAAACCTTCGGTTCTATTCCCGTTTTTCTCTTTTCATTCTTCATGAGTTAATTCCACTGTGAAAACGAATATTTAATTTGCGTTCCCGTTTGTTTTATTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGGGTAGTCAAGAAGCTGTATTCTGTGATTGGTAGCACATTTTTGTTTCGGGTGAATAGACGGCGttta                                                                                                                                                                                                                                 <  2:112450/251‑2 (MQ=255)
             ggcggAAAGAAAACCTTCGGTTCTATTCCCGTTTTTCTCTTTTCATTCTTCATGAGTTAATTCCACTGTGAAAACGAATATTTATTTTGCGTTCCCGTTTGTTTTATTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTg                                                                                                                                                                                                                    >  1:423158/1‑251 (MQ=255)
                                gTTCTATTCCCGTTTTTCTCTTTTCATTCTTCATGAGTTAATTCCACTGTGAAAACGAATATTTATTTTGCGTTCCCGTTTGTTTTATTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGtt                                                                                                                                                                                                  >  1:234843/1‑250 (MQ=255)
                                                           tcttcATGAGTTAATTCCACTTTGAAAACGAATATTTATTTTGCGTTCCCGTTTGTTTTATTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGGTTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCTGGTGAATAGAATGCGTTTTGTCGTTAATTTTGATTAATAATCAGTTTGTTTTGCTCTGTTGTGAGTAAATAATAACa                                                                                                                                                                       <  2:38567/250‑1 (MQ=255)
                                                                       aaTTCCACTGGGAAAACGAATATTTATTTTGCGTTCCCGTTTGTTTTATTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAATTCATAATGGGTATTCAATATGCTGTATTCTGTGCTTGGTATCACATTTTTGTTTCGGGAGACTAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAAAATATGACTTTCAat                                                                                                                                                          >  2:368927/1‑251 (MQ=255)
                                                                            cACTGTGAAAACGAATATTTATTTTTAGTTCACGTTTGTTTTAGTTTTGTTAACATTTAATATAATTATTACTAACCTCGTGGACGCGGTAATGGCTATCTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCTTTTTTTCGTTAATTTTGATTTATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTAAATATTgg                                                                                                                                                     <  2:349360/251‑1 (MQ=255)
                                                                                            atttatttTGCGTTCCCGTTTGTTTTATTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATTGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTGCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATCCATAAAACAat                                                                                                                                     <  1:306973/251‑1 (MQ=255)
                                                                                            atttatttTGCGTTCCCGTTTGTTTTATTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATCCATAAAACAat                                                                                                                                     >  1:91569/1‑251 (MQ=255)
                                                                                                tttttGCGTTACCTTTTGTTTTATTTTTGTTAACATTTAATATATTTATTATTACCATCGTGGACGCGTTAATATCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATAACATTTTTGTTACGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTTATTTATAATCAGTTTGTTATGCGCGGTTGTGAGTTAAAATTAATATATGACTTTCATTATTGTTGCTCCATCAAACAATATTg                                                                                                                                 >  2:268135/2‑251 (MQ=255)
                                                                                                             gtttgtttTATTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTTATCCATAAAACAATATTGAAAATTTCttttt                                                                                                                    <  2:351448/251‑1 (MQ=255)
                                                                                                                tgtttTATTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAATTAATAATGTGTTTTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGTTTAATTATCAGTTTGTTATGCTCTGTTGTGAGTTAAAAATAACATATGACTTTAATTATTGGTGTTCTATAAAACAATATTGAAAATTTCTTTTTGCt                                                                                                                 >  2:258943/1‑251 (MQ=255)
                                                                                                                          ttGTTAACATTTAATATAATTATTATTAACCTCGTGTACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATCCATAAAACAATATTGAAAATTTCTTTTTGCTACGCCgtgt                                                                                                        <  2:33009/250‑1 (MQ=255)
                                                                                                                              tAACATTTAATATAATTATTATTAACCTCGTGGACGCTTTAATGGCTAACTCATAATGGGTATTCAATAATCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATCCATAAAACAATATTGAAAATTTCTTTTTGCTACGCCGTGTTTTc                                                                                                    <  2:123382/250‑1 (MQ=255)
                                                                                                                                         ataattactatcaccctaggggaagcgTTAAAGGCTAACTCCTACTGGGTATTCACTAAGCCGCATTCTGTGATTAGTATATCATTTTCGTTTCGGGTGAATAGAAGGCGTTTTTTCGTCAATTTTGATTAATAATCAGTTTGTTCTGCTCTGTTGAGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATCCATAAAACAATATTGAAAATTTATTTTTGCTACGCCGTGTTTTCAATATTGGTg                                                                                          >  2:36252/25‑249 (MQ=255)
                                                                                                                                           aaTTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATCCATAAAACAATATTGAAAATTTCTTTTTGCTACGCCGTGTTTTCAATATTGGTGAGGa                                                                                      <  1:59868/251‑1 (MQ=255)
                                                                                                                                                                    ttAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATCCATAAAACAATATTGAAAATTTCTTTTTGCTACGCCGTGTTTTCAATATTGGTGAGGAACTTAACAATATTGAAAGTTGGAtt                                                             >  2:108465/1‑251 (MQ=255)
                                                                                                                                                                                                          ttCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATCCATAAAACAATATTGAAAATTTCTTTTTGCTACGCCGTGTTTTCAATATTGGTGAGGAACTTAACAATATTGAAAGTTGGATTTATCTGCGTGTGACATTTTCAATATTGGTGATTAAAGt                       >  1:172976/1‑251 (MQ=255)
                                                                                                                                                                                                              gtgATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCGCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATCCATAAAACAATATTGAAAATTTCTTTTTGCTACGCCGTGTTTTCAATATTGGTGAGGAACTTAACAATATTGAAAGTTGGATTTATCTGCGTGTGACATTTTCAATATTGGTGATTAAAGtttt                    <  1:301531/250‑1 (MQ=255)
                                                                                                                                                                                                                   tGGTATCACATTTTTTTTTCGGGTGAATAGAAGGCGTTGTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATCCATAAAACAATATTGAAAATTTCTTTTTGCTACGCCGTGTTTGCAATATTGGTGAGGAACTTAACAATATTGAAAGTTGGATTTATCTGCGTGTGACATGTTCAATATTGGTGATTAAAGTTTTAttta               <  2:302968/250‑2 (MQ=255)
                                                                                                                                                                                                                                ttGTGTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATCCATAAAACAATATTGAAAATTTCTTTTTGCTACGCCGTGTTTTCAATATTGGTGAGGAACTTAACAATATTGAAAGTTGGATTTATCTGCGTGTGACATTTTCAATATTGGTGATTAAAGTTTTATTTCAAAATTAAAGGGc  <  2:250300/250‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                       
GACGATTATTAACGGCGGAAAGAAAACCTTCGGTTCTATTCCCGTTTTTCTCTTTTCATTCTTCATGAGTTAATTCCACTGTGAAAACGAATATTTATTTTGCGTTCCCGTTTGTTTTATTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAGGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATCCATAAAACAATATTGAAAATTTCTTTTTGCTACGCCGTGTTTTCAATATTGGTGAGGAACTTAACAATATTGAAAGTTGGATTTATCTGCGTGTGACATTTTCAATATTGGTGATTAAAGTTTTATTTCAAAATTAAAGGGC  >  NC_000913/41866‑42339

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