Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,002,709 A→G 8.9% W31R (TGG→CGG)  yigF ← DUF2628 family putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,002,7090AG 8.9% 141.6 / 4.4 45W31R (TGG→CGG) yigFDUF2628 family putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (20/21);  new base G (2/2);  total (22/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.92e-01

TACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAGATAGAATACCAAGTAATAGAACTATAATTGCTTTTTTCCATAACCCCAGGACAAATAAATATATCCAGGAACAGAAGAAAGCAATGAAATTCATTTGAATAGTTAAGCGTTGTCTTACTTTTAATGCTTTAAATGCTGCTTTATATTCTGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAACTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATAGAG  >  NC_000913/4002497‑4002939
                                                                                                                                                                                                                    |                                                                                                                                                                                                                                      
tACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAGATAGAATACCAAGTAATAGAACTATAATTGCTTTTTTCCATAACCCCAGGACAAATAAATATATCCAGGAACAGAAGAAAGCAATGAAATTCATTTGAATAGTTAAGCGTTGTCTTACTTTTAATGCTTTAAATGCTGCTTTATATTCTGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATg                                                                                                                                                                          >  4:56375/1‑275 (MQ=255)
          cccAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAGATAGAATACCAAGTAATAGAACTATAATTGCTTTTTTCCATAACCCCAGGACAAATAAATATATCCAGGAACAGAAGAAAGCAATGAAATTCATTTGAATAGTTAAGCGTTGTCTTACTTTTAATGCTTTAAATGCTGCTTTATATTCTGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCg                                                                                                                                                                  <  2:56376/273‑1 (MQ=255)
                                                                                   ttgtCCATANCCCCAGGACAAATAAATATATCCAGGAACAGAAGAAAGCAATGAAATTCATTTGAATAGTTAAGCGTTGTCTTACTTTTAATGCTTTAAATGCTGCTTTATATTCTGGTGTTGCCCCCCAGAATCCAGGAAAACCATg                                                                                                                                                                                                                      <  3:137103/145‑1 (MQ=255)
                                                                                                          aaataTATCCAGGAACAGAAGAAAGCAATGAAATTCATTTGAATAGTTAAGCGTTGTCTTACTTTTAATGCTTTAAATGCTGCTTTATATTCTGGTGTTGCCCCCCAGAATCCa                                                                                                                                                                                                                                 >  8:442106/1‑114 (MQ=255)
                                                                                                          aaataTATCCAGGAACAGAAGAAAGCAATGAAATTCATTTGAATAGTTAAGCGTTGTCTTACTTTTAATGCTTTAAATGCTGCTTTATATTCTGGTGTTGCCCCCCAGAATCCa                                                                                                                                                                                                                                 <  7:442106/114‑1 (MQ=255)
                                                                                                            atatatCCAGGAACAGAAGAAAGCAATGAAATTCATTTGAATAGTTAAGCGTTGTCTTACTTTTAATGCTTTAAATGCTGCTTTATATTCTGGTGTTGCCCCCCAGAATCCa                                                                                                                                                                                                                                 >  6:820842/1‑112 (MQ=255)
                                                                                                            atatatCCAGGAACAGAAGAAAGCAATGAAATTCATTTGAATAGTTAAGCGTTGTCTTACTTTTAATGCTTTAAATGCTGCTTTATATTCTGGTGTTGCCCCCCAGAATCCa                                                                                                                                                                                                                                 <  5:820842/112‑1 (MQ=255)
                                                                                                                     ggAACAGAAGAAAGCAATGAAATTCATTTGAATAGTTAAGCGTTGTCTTACTTTTAATGCTTTAAATGCTGCTTTATATTCTGGTGTTGCCCCCCAGAATcc                                                                                                                                                                                                                                  >  8:420930/1‑102 (MQ=255)
                                                                                                                     ggAACAGAAGAAAGCAATGAAATTCATTTGAATAGTTAAGCGTTGTCTTACTTTTAATGCTTTAAATGCTGCTTTATATTCTGGTGTTGCCCCCCAGAATcc                                                                                                                                                                                                                                  <  7:420930/102‑1 (MQ=255)
                                                                                                                      gAACAGAAGAAAGCAATGAAATTCATTTGAATAGTTAAGCGTTGTCTTACTTTTAATGCTTTAAATGCTGCTTTATATTCTGGTGTTGCCCCCCAGAATCCAGGAAAACCa                                                                                                                                                                                                                        >  2:208377/1‑111 (MQ=255)
                                                                                                                          agaagaAAGCAATGAAATTCATTTGAATAGTTAAGCGTTGTCTTACTTTTAATGCTTTAAATGCTGCTTTATATTCTGGTGTTGCCCCCCAGAATCCAGGAAAACCATGtt                                                                                                                                                                                                                    <  4:208376/111‑1 (MQ=255)
                                                                                                                                                    aTAGTTAAGCGTTGTCTTACTTTTAATGCTTTAAATGCTGCTTTATATTCTGGTGTTGCCCCCCAGaa                                                                                                                                                                                                                                     <  3:138448/68‑1 (MQ=255)
                                                                                                                                                                tGTCTTACTTTTAATGCTTTAAATGCTGCTTTATATTCTGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAACTATTGTTTTGATTAATCAGTAGGTGCAAGCATTCGCATACTGAAAGTGGAAAACTACCAACTCAAAAAAATCATCGAACCATTGCCTGCACAGGCAAAATCTTCGTCTAACATTGtga          >  4:194608/1‑275 (MQ=255)
                                                                                                                                                                   cTTACTTTTAATGCTTTAAATGCTGCTTTATATTCTGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCAt                                                                                                                                                                                       <  5:927638/99‑1 (MQ=255)
                                                                                                                                                                   cTTACTTTTAATGCTTTAAATGCTGCTTTATATTCTGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCAt                                                                                                                                                                                       >  6:927638/1‑99 (MQ=255)
                                                                                                                                                                        ttttAATGCTTTAAATGCTGCTTTATATTCTGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAACTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATagag  <  1:97860/275‑1 (MQ=255)
                                                                                                                                                                         tttAATGCTTTAAATGCTGCTTTATATTCTGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATaaaaa                                                                                                                                                                                                         <  7:979398/75‑1 (MQ=255)
                                                                                                                                                                         tttAATGCTTTAAATGCTGCTTTATATTCTGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATaaaaa                                                                                                                                                                                                         >  8:979398/1‑75 (MQ=255)
                                                                                                                                                                                 tttAAATGCTGCTTTATATTCTGGTGTTGCCCCCCGGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCAt                                                                                                                                                                    <  6:1080556/104‑1 (MQ=255)
                                                                                                                                                                                 tttAAATGCTGCTTTATATTCTGGTGTTGCCCCCCGGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCAt                                                                                                                                                                    <  8:894266/104‑1 (MQ=255)
                                                                                                                                                                                 tttAAATGCTGCTTTATATTCTGGTGTTGCCCCCCGGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCAt                                                                                                                                                                    >  5:1080556/1‑104 (MQ=255)
                                                                                                                                                                                 tttAAATGCTGCTTTATATTCTGGTGTTGCCCCCCGGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCAt                                                                                                                                                                    >  7:894266/1‑104 (MQ=255)
                                                                                                                                                                                 tttAAATGCTGCTTTATATTCTGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCg                                                                                                                                                                  >  8:146290/1‑106 (MQ=255)
                                                                                                                                                                                 tttAAATGCTGCTTTATATTCTGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCg                                                                                                                                                                  <  7:146290/106‑1 (MQ=255)
                                                                                                                                                                                    aaaTGCTGCTTTATATTCTGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGaa                                                                                                                                                                        <  5:962086/97‑1 (MQ=255)
                                                                                                                                                                                    aaaTGCTGCTTTATATTCTGGTGTTGCCCCCCAGAATCCAGGAAAACCATATTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGaa                                                                                                                                                                        >  6:962086/1‑97 (MQ=255)
                                                                                                                                                                                      aTGCTGCTTTATATTCTGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCtt                                                                                                                                                    <  5:862950/115‑1 (MQ=255)
                                                                                                                                                                                      aTGCTGCTTTATATTCTGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCtt                                                                                                                                                    >  6:862950/1‑115 (MQ=255)
                                                                                                                                                                                         ctgctTTATATTCTGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATAtt                                                                                                                                                        <  8:1039548/108‑1 (MQ=255)
                                                                                                                                                                                         ctgctTTATATTCTGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATAtt                                                                                                                                                        >  7:1039548/1‑108 (MQ=255)
                                                                                                                                                                                                    tCTGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATaaaaa                                                                                                                                                                                                         >  6:129700/1‑48 (MQ=255)
                                                                                                                                                                                                    tCTGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATaaaaa                                                                                                                                                                                                         <  5:129700/48‑1 (MQ=255)
                                                                                                                                                                                                    tCTGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAACTATTGttt                                                                                                                        <  5:419313/129‑1 (MQ=255)
                                                                                                                                                                                                    tCTGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAACTATTGttt                                                                                                                        >  6:419313/1‑129 (MQ=255)
                                                                                                                                                                                                            tGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAACTATTGtt                                                                                                                         >  5:214683/1‑120 (MQ=255)
                                                                                                                                                                                                            tGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAACTATTGtt                                                                                                                         <  6:214683/120‑1 (MQ=255)
                                                                                                                                                                                                               cccccAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCAt                                                                                                                                              <  5:140692/96‑1 (MQ=255)
                                                                                                                                                                                                               cccccAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCAt                                                                                                                                              >  6:140692/1‑96 (MQ=255)
                                                                                                                                                                                                                 cccAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATa                                                                                                                                  <  5:649311/106‑1 (MQ=255)
                                                                                                                                                                                                                 cccAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATa                                                                                                                                  >  6:649311/1‑106 (MQ=255)
                                                                                                                                                                                                                  ccAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACtt                                                                                                                                                                              >  1:29532/1‑61 (MQ=255)
                                                                                                                                                                                                                  ccAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACt                                                                                                                                       <  7:965377/100‑1 (MQ=255)
                                                                                                                                                                                                                  ccAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACt                                                                                                                                       >  8:965377/1‑100 (MQ=255)
                                                                                                                                                                                                                  ccAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATa                                                                                                                                  <  7:851290/105‑1 (MQ=255)
                                                                                                                                                                                                                  ccAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATa                                                                                                                                  >  8:851290/1‑105 (MQ=255)
                                                                                                                                                                                                                    |                                                                                                                                                                                                                                      
TACGGCAATTCCCAGAATGTTAACGCCAATTAGTGCACCAACAAATAGAGATAGAATACCAAGTAATAGAACTATAATTGCTTTTTTCCATAACCCCAGGACAAATAAATATATCCAGGAACAGAAGAAAGCAATGAAATTCATTTGAATAGTTAAGCGTTGTCTTACTTTTAATGCTTTAAATGCTGCTTTATATTCTGGTGTTGCCCCCCAGAATCCAGGAAAACCATGTTGATCATAAAAATTAAATCGGTATTTCCATTTTTCACTTAATGAACCATCGTTCATATATTCCTTACTCATAAATACTCCATAACTATTGTTTTGATGAATCAGTAGGTGCAAGCATTAGCATACTGAAAGTGGAAAAATAACAAATCAAAAAAATCATCGAACCATTGCCTGAACAGGCAAAATCTTCGGCTATCATTGTGATGATAGAG  >  NC_000913/4002497‑4002939

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