Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,076,069 A→T I431F (ATC→TTC)  yeeR → CP4‑44 prophage; inner membrane protein YeeR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,076,0690AT94.4% 53.9 / NA 18I431F (ATC→TTC) yeeRCP4‑44 prophage; inner membrane protein YeeR
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (9/8);  minor base G (1/0);  total (10/8)
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

GGTCTGATTCTCATGAAACTGCCGGAAACCTGGAGTACAAATGAGAAAATGTTTGCTTCCGGAGGGCAGGGACATGGGTTTGCCGCTGAACGGGGAAACCATATTGTCGACAGAGTCCGTCTGAAAAACGCACGGATCCTCGGTGATAATAATGCCAGAAATGGAGCAGACAGACTGGTCAGCGGAACAGAAATCCAGACGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACGGACAGTATCG  >  NC_000913/2075934‑2076199
                                                                                                                                       |                                                                                                                                  
ggTCTGATTCTCATGAAACTGCCGGAAACCTGGAGTACAAATGAGAAAATGTTTGCTTCCGGAGGGCAGGGACATGGGTTTGCCGCTGAACGGGGAAACCATATTGTCGACAGAGTCCGTCTGAAAAACGCACGGTTCCt                                                                                                                                <  1:127876/140‑1 (MQ=255)
ggTATGATTCTCATGAAACTGCCGGAAACCTGGATTACAAATGAGAAAATGTTTGCTTCCGGAGGGCAGGGACATGGGTTTGCCGCTGAACGGGGAAACCATATTGTCGACAGAGTCCGTCTGAAAAACGCACGGTTCCt                                                                                                                                <  2:170859/140‑1 (MQ=255)
 gTCTGATTCTCATGAAACTGCCGGAAACCTGGAGTACAAATGAGAAAATGTTTGCTTCCGGAGGGCAGGGACATGGGTTTGCCGCTGAACGGGGAAACCATATTGTCGACAGAGTCCGTCTGAAAAACGCACGGTTCCTc                                                                                                                               <  2:158818/140‑1 (MQ=255)
 gTCTGATTCTCATGAAACTGCCGGAAACCTGGAGTACAAATGAGAAAATGTTTGCTTCCGGAGGGCAGGGACATGGGTTTGCCGCTGAACGGGGAAACCATATTGTCGACAGAGTCCGTCTGAAAAACGCACGGTTCCTc                                                                                                                               <  1:90071/140‑1 (MQ=255)
            aTGAAACTGCCGGAAACCTGGAGTACAAATGAGAAAATGTTTGCTTCCGGAGGGCAGGGACATGGGTTTGCCGCTGAACGGGGAAACCATATTGTCGACAGAGTCCGTCTGAAAAACGCACGGTTCCTCGGTGAtaataa                                                                                                                    >  2:284277/1‑140 (MQ=255)
                                       aaTGAGAAAATGTTTGCTTCCGGAGGGCAGGGACATGGGTTTGCCGCTGAACGGGGAAACCATATTGTCGACAGAGTCCGTCTGAAAAACGCACGGTTCCTCGGTGATAATAATGCCAGAAATGGAGCAGACAGACTGGt                                                                                         <  2:155373/140‑1 (MQ=255)
                                                          ccGGAGGGCAGGGACATGGGTTTGCCGCTGAACGGGGAAACCATATTGTCGACAGAGTCCGTCTGAAAAACGCACGGGTCCTCGGAGATAATAATGCCAGAAATGGAGCAGACAGACTGGTCAGCGGAACAGAAATCCAg                                                                      >  1:119415/1‑140 (MQ=255)
                                                                               ttGCCGCTGAACGGGGAAACCATATTGTCGACAGAGTCCGTCTGAAAAACGCACGGTTCCTCGGTGATAATAATGCCAGAAATGGAGCAGACAGACTGGTCAGCGGAACAGAAATCCAGACGAAATATTGTTCAACTGCa                                                 >  2:76875/1‑140 (MQ=255)
                                                                                tGCCGCTGAACGGGGAAACCATATTGTCGACAGAGTCCGTCTGAAAAACGCACGGTTCCTCGGTGATAATAATGCCAGAAATGGAGCAGACAGACTGGTCAGCGGAACAGAAATCCAGACGAAATATTGTTCAACTGCAg                                                >  2:147750/1‑140 (MQ=255)
                                                                                        aaCGGGGAAACCATATTGTCGACAGAGTCCGTCTGAAAAACGCACGGTTCCTCGGTGATAATAATGCCAGAAATGGAGCAGACAGACTGGTCAGCGGAACAGAAATCCAGACGAAATATTGTTCAACTGCAGCCCGTAGc                                        >  2:144701/1‑140 (MQ=255)
                                                                                                          tCGACAGAGTCCGTCTGAAAAACGCACGGTTCCTCGGTGATAATAATGCCAGAAATTGAGCAGACAGACTGGTCAGCGGAACAGAAATCCAGACGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGAc                      <  1:286063/140‑1 (MQ=255)
                                                                                                            gACAGAGTCCGTCTGAAAAACGCACGGTTCCTCGGTGATAATAATGCCAGAAATGGAGCAGACAGACTGGTCAGCGGAACAGAAATCCAGACGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACgg                    >  2:171239/1‑140 (MQ=255)
                                                                                                            gACAGAGTCCGTCTGAAAAACGCACGGTTCCTCGGTGATAATAATGCCAGAAATGGAGCAGACAGACTGGTCAGCGGAACAGAAATCCAGACGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACgg                    >  2:197689/1‑140 (MQ=255)
                                                                                                             aCAGAGTCCGTCTGAAAAACGCACGGTTCCTCGGTGATAATAATGCCAGAAATGGAGCAGACAGACTGGTCAGCGGAACAGAAATCCAGACGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGa                   >  2:284260/1‑140 (MQ=255)
                                                                                                                 agTCCGTCTGAAAAACGCACGGTTCCTCGGTGATAATAATGCCAGAAATGGAGCAGACAGACTGGTCAGCGGAACAGAAATCCAGACGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACaaa               >  1:17577/1‑140 (MQ=255)
                                                                                                                    ccGTCTGAAAAACGCACGGTTCCTCGGTGATAATAATGCCAGAAATGGAGCAGACAGACTGGTCAGCGGAACAGAAATCCAGACGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACg            >  2:109144/1‑140 (MQ=255)
                                                                                                                       tCTGAAAAACGCACGGTTCCTCGGTGATAATAATGCCAGAAATGGAGCAGACAGACTGGTCAGCGGAACAGAAATCCAGACGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACGGAc         <  1:171784/140‑1 (MQ=255)
                                                                                                                              aaCGCACGGTTCCTCGGTGATAATAATGCCAGAAATGGAGCAGACAGACTGGTCAGCGGAACAGAAATCCAGACGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACGGACAGTATCg  <  1:24102/140‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                  
GGTCTGATTCTCATGAAACTGCCGGAAACCTGGAGTACAAATGAGAAAATGTTTGCTTCCGGAGGGCAGGGACATGGGTTTGCCGCTGAACGGGGAAACCATATTGTCGACAGAGTCCGTCTGAAAAACGCACGGATCCTCGGTGATAATAATGCCAGAAATGGAGCAGACAGACTGGTCAGCGGAACAGAAATCCAGACGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACGGACAGTATCG  >  NC_000913/2075934‑2076199

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: