Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,227,339 C→T T369I (ACC→ATC)  ygjI → putative transporter YgjI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,227,3390CT83.3% 36.2 / ‑2.1 19T369I (ACC→ATC) ygjIputative transporter YgjI
Reads supporting (aligned to +/- strand):  ref base C (0/3);  new base T (9/6);  total (9/10)
Fisher's exact test for biased strand distribution p-value = 2.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCG‑TCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATT‑ATTAATATGACCGCCGCAGCGTCCATG‑CTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCC‑AAATTA‑GATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGTCGCTTCGACATTCCCGACTGGCGCGAATATTCTGACCATCATTTTTTATAACGTCGGCGGTATTG  >  NC_000913/3227105‑3227563
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                   
ggtCCTTTTCCGCGTTGTTCGGTTCCCTGCTGATGGGGACCGGAACGCCGGGGAAAATGTTCTTCTCCGAAATCCCGGAAGGCATCGTTGGTAAGAAAACCGTCGAACTGAACGATAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCG‑TCATCCCGCTGATGATTATCCCGATGCGCGGTTCCAATACCGTGCAGGATCTGATGAATACTATT‑ATTAAGATGAGCGTCGCAGCGTCCATg                                                                                                                                                                                                                     <  2:394523/250‑1 (MQ=255)
      ttACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCG‑TCATCCCGCTGATGATTATCCCGATGCTCGGTTCCACTACCGTGCAGGATCTGATGAATACTATT‑ATTAATATGATCGCCGCAGCGTCCATG‑CTTccg                                                                                                                                                                                                              >  1:31270/1‑251 (MQ=255)
        aCCGCGATGTTCGGTTCCCTGCTGATGTGGAACGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGCAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCG‑TCATCACGCTGCTGATTATCCCGATGCTCGGTTCCAATAACGTGCAGGATCTGATGAATACTATT‑ATTAATATGCTCGCCGCAGCGTCCATG‑CTTccgcc                                                                                                                                                                                                            >  1:142486/1‑251 (MQ=255)
                   cGGTGCCCTGCTGATGTGGACCGCAACGCCGGTCAAATTTGTCTTCTCCGAAAACCCGGAAGGCCTCTTGGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCG‑TCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTTTT‑ATTAATATGATCGCCGCAGCGTCCATG‑CGTCCGCCGGTATtcatca                                                                                                                                                                                                 <  2:386447/251‑1 (MQ=255)
                                                        aTTTTCTTCTCCAAAAACCCGGGAGGCATCTTTTGGAAAAAAACCCGCGAACTGAAAGAAAACGGCGTGCCGGGGCGCGCAGGGGTGATCCATTTCCTGATCGGTCA‑CCCGCTGATGATTATCCCGATGCTCGGTTCCAATCCCGTGCAGGATCTGATGAATACTATT‑ATTAATATGATCGCCGCAGCGTCCATG‑CTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCC‑AAATTA‑GATCa                                                                                                                                                          <  2:13001/251‑1 (MQ=255)
                                                         ttttCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCG‑TCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATT‑ATTAATATGATCGCCGCAGCGTCCATG‑CTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCC‑AAATTA‑GATCAc                                                                                                                                                         >  2:321027/1‑251 (MQ=255)
                                                          tttCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCG‑TCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATT‑ATTAATATGATCGCCGCAGCGTCCATG‑CTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCC‑AAATTA‑GATCAc                                                                                                                                                         <  1:121302/250‑1 (MQ=255)
                                                               gttcctaaacccctgAAGCCATCTTTGGTAAGAAAACCGCCGAACGGAACGAAAACGGCGTTCCGGCGCGCGCATCGTGGATCCAGTTCCTGATCG‑TCATCCCGCTGATGATTATCCCGATTCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATT‑ATTAATATGATCGCCGCATCGTCCATG‑CCTCCGCCGTTATTCATCATGCTGGCTTACCTGAAGTTACGCGCC‑AAATTA‑GATCTCCTGc                                                                                                                                                     <  2:329517/235‑1 (MQ=255)
                                                                                       ttGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCG‑TCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATT‑ATTAATATGATCGCCGCAGCGTCCATG‑CTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCC‑AAATTA‑GATCACCTGCCACGCGATTTCCGTATGGGCTCccgc                                                                                                                           >  1:70768/1‑251 (MQ=255)
                                                                                          gTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCG‑TCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATT‑ATTAATATGATCGCCGCAGCGTCCATG‑CTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCC‑AAATTA‑GATCACCTGCCACGCGATTTCCGTATGGGCTCCcgccgc                                                                                                                        >  2:9280/1‑251 (MQ=255)
                                                                                                             gaccGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTACTGATCG‑TCATCCAGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATT‑ATTAATATGATCGCCGCAGCGTCCATG‑CTTCCGCCGTTATTCTTCATGCTGGCTTACCTGAATTTACGCGCC‑AAATTA‑GATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGACCCGGTATTATCGTTGttt                                                                                                     >  2:176455/4‑251 (MQ=255)
                                                                                                                          gTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCG‑TCATCCCGCTGCTGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATT‑ATTAATATGATCGCCGCAGCGTCCATG‑CTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCC‑AAATTA‑GATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGc                                                                                        >  2:143461/1‑251 (MQ=255)
                                                                                                                                  gcgcgcAGCGTGGATCCAGTTCCTGATCG‑TCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATCCTATT‑ATTAATATGATCGCCGCAGCGTCCATG‑CTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCC‑AAATTA‑GATCCCCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTg                                                                                <  1:321027/251‑1 (MQ=255)
                                                                                                                                         gggtggagcatgtcccTGATGG‑TCATCCCGCTTATGATTACCCCGATGCTCGGTACAAATCCCGTGCTGGATCTGA‑GAATTATATTTATAAAGAGGGTCCCGGCAGCTTCTAGC‑CTGCCTCGTTTATTCATCATGCTGGCTTCCCTAATTTTACGCGCA‑AGATTA‑GATCCCCTGCCACGCGATTTCCGTGGGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCATTGCTGATTGCGATAGTTGCCGTAgg                                                                         <  2:166034/237‑1 (MQ=255)
                                                                                                                                                 ccAGTTCCTGATCG‑TCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATT‑ATTAATATGATCGCCGCAGCGTCCATG‑CTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCC‑AAATTA‑GATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGTCg                                                                 >  1:236408/1‑251 (MQ=255)
                                                                                                                                                                                     cgaggCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATT‑ATTAATATGATCGCCGCAGCGTCCATG‑CTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCC‑AAATTA‑GATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTGTGTCGCTTCGACATTCCCGACTGGCGCGACTATTCTGAcc                              >  1:187624/5‑251 (MQ=255)
                                                                                                                                                                                            gggTTCCAATACCGTGCGGTATCGGATGAATACTTTT‑ATTAAGATGACCCCCGCAGCGTC‑ATGTCTCCCGCCTTTATTCATCATTCGGGCTTACCTTAATTTACGCGCC‑AAATTA‑GATCCCCGGCCACGCGATTCCCGTAGGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGTCGCTTCGACATTCCCGACTGGCGCGAATATTCTGACCATCAttt                       <  1:143461/250‑1 (MQ=255)
                                                                                                                                                                                                                    tatgatttttttt‑attattaTGACCGCCCCATCTTCCACG‑CTTCCCCAGTTATTCATCATGCTGGCTTACCT‑AATTT‑CTCGCCCAAAATATGATCCCCTCCCACGCAGTTTCCGTATGTGCTCCCGCCGCACCGTTATTGTCGTTGTTTCAATGCTGATTGCGATAGTTGCCGTAGGGTTGGTCGCTTCGACATTCCCGACTGGCGCGAATATTCTGACCATCATTTTTTATAACGTCGGCGGTATTg  <  2:236408/237‑1 (MQ=255)
                                                                                                                                                                                                                    gatgaATACTATT‑ATAAATATGTCCGCCGCGGCCTCAATG‑TTCCCTCCGTTATTCATCATGCTGGTTTCCATTAATTTCCGCGCC‑AATTTA‑GACCCCCTGCCCCGCGTTTTCCGTAGTGCCGCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGTCGCTTCGACATTCCCGACTGGCGCGAATATTCTGACCATCATTTTTTATAACGTCGGCGGTATTg  <  1:163972/248‑1 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                   
TGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCG‑TCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATT‑ATTAATATGACCGCCGCAGCGTCCATG‑CTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCC‑AAATTA‑GATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGTCGCTTCGACATTCCCGACTGGCGCGAATATTCTGACCATCATTTTTTATAACGTCGGCGGTATTG  >  NC_000913/3227105‑3227563

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