Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,088,193 C→A intergenic (+333/‑91) metK → / → galP S‑adenosylmethionine synthetase/D‑galactose transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,088,1930CA97.1% 108.5 / ‑4.2 34intergenic (+333/‑91)metK/galPS‑adenosylmethionine synthetase/D‑galactose transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (19/14);  minor base G (1/0);  total (20/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

TTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCT  >  NC_000913/3088048‑3088338
                                                                                                                                                 |                                                                                                                                                 
ttttCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTTAATTATTTTCATGCACTTAAATCATAACTAAGATAAAGCTTATTGTAAGAGAt                                                                                                                                                <  2:485794/149‑1 (MQ=255)
  catatatCATGATTCGGTTTACACAGGTGTTACATTTCTTTTCAGTAAAGTCTTAAGTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAACCATAACTAAGATAAATGTTAGTGTAAGAGATTa                                                                                                                                              <  1:283826/146‑1 (MQ=255)
           ttattCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGAtgt                                                                                                                                     <  1:383053/149‑1 (MQ=255)
                    tttCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTACTGCACGTAAATCATAACTAATATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCtt                                                                                                                            <  1:309762/149‑1 (MQ=255)
                        acaGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTcaca                                                                                                                        >  1:596095/1‑149 (MQ=255)
                          aGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATc                                                                                                                      >  1:68526/1‑149 (MQ=255)
                                    tttCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTAcgtc                                                                                                            <  2:105300/149‑1 (MQ=255)
                                        ttttCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTAc                                                                                                        <  2:345596/149‑1 (MQ=255)
                                                    tCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGCTGTGATTTGATTCACATCTTTTTACGTCGTACTCACCTATCTTa                                                                                            >  2:497937/1‑149 (MQ=255)
                                                            gCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTcgaaa                                                                                    >  1:85438/1‑145 (MQ=255)
                                                                aTAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTGAATTCACAATcaa                                                                                >  1:553627/1‑149 (MQ=255)
                                                                   aCAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCa                                                                                                     >  1:65397/1‑125 (MQ=255)
                                                                   aCAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCa                                                                                                     <  2:65397/125‑1 (MQ=255)
                                                                         tttAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGGGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCa                                                                       >  1:597862/1‑149 (MQ=255)
                                                                         tttAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACATTGAAAAATAACCa                                                                       >  1:373699/1‑149 (MQ=255)
                                                                                  tgatgaTATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAgg                                                              <  2:117861/149‑1 (MQ=255)
                                                                                           aagtCAATTATTGTCAGGCACTTAAATCATAATTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGc                                                     <  2:609013/146‑1 (MQ=255)
                                                                                           aaCTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATa                                                                           >  2:81246/1‑127 (MQ=255)
                                                                                           aaCTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATa                                                                           <  1:81246/127‑1 (MQ=255)
                                                                                           aaCTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGc                                                     >  2:610208/1‑149 (MQ=255)
                                                                                              tCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTATTAACCTATCTTAATTCGCAATAAAAAATAACCAAATTAGAGGGAATCATGCCTg                                                  >  1:180569/1‑149 (MQ=255)
                                                                                              tCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTg                                                  >  1:44486/1‑149 (MQ=255)
                                                                                              tCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTg                                                  >  1:205207/1‑149 (MQ=255)
                                                                                                aaatattTTCATGCACTTAAATCATAACTAAGATAAATGTTTGTGTAAGAAATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGAc                                                <  2:362759/146‑1 (MQ=255)
                                                                                                              aCTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAggg                                  >  2:305692/1‑149 (MQ=255)
                                                                                                               cTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAgggg                                 >  1:775519/1‑149 (MQ=255)
                                                                                                                           cTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTAACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAg                     <  2:546579/149‑1 (MQ=255)
                                                                                                                           cTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACACACCTATCTTAATTCACCATAAAAAATAACCATATTGGAGGGCATCATGCCTGCAGCTAAAAAACAGGGCCGGTCAAACacg                     >  2:483451/1‑147 (MQ=255)
                                                                                                                             aaGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTAACAATAAAAAATAACCATATTGGAGGGCATCATGACTGACGCTCAAAAACAGGGCCGGTCAAACCAGGc                   >  2:149882/1‑149 (MQ=255)
                                                                                                                                        tAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGtttt        >  1:206323/1‑149 (MQ=255)
                                                                                                                                           tgtAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTCTCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGCCGCTAAAAAACAGGGGCGGTCAAAGAAGGCAATGACATTTTTCg     >  1:497595/1‑149 (MQ=255)
                                                                                                                                              aaGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCa                  >  2:214259/1‑133 (MQ=255)
                                                                                                                                              aaGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCa                  <  1:214259/133‑1 (MQ=255)
                                                                                                                                              aaGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCt  <  2:44486/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
TTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCT  >  NC_000913/3088048‑3088338

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