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,1930CA98.0% 168.2 / NA 50intergenic (+333/‑91)metK/galPS‑adenosylmethionine synthetase/D‑galactose transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (25/24);  minor base T (0/1);  total (25/25)
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

CATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTC  >  NC_000913/3088052‑3088337
                                                                                                                                             |                                                                                                                                                
cATAGTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTaca                                                                                                                                           <  2:798311/149‑1 (MQ=255)
          ttCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGAt                                                                                                                                 >  2:468724/1‑149 (MQ=255)
          ttCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGAAGTGAt                                                                                                                                 >  2:358391/1‑149 (MQ=255)
           tCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGAtt                                                                                                                                <  2:278109/149‑1 (MQ=255)
              gTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGc                                                                                                                             <  1:178202/149‑1 (MQ=255)
                 ttCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGCTATAACTAAATTATTTTCATGCACTTCAATCATAACTAATATAAATGTTTGTGTAAGAGGCAGCACTGCTGTGATTTGCTTc                                                                                                                          >  1:434363/1‑149 (MQ=255)
                   cacaGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTcac                                                                                                                        <  2:695822/149‑1 (MQ=255)
                      aGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATc                                                                                                                     >  2:414289/1‑149 (MQ=255)
                           ttACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCttttt                                                                                                                >  1:443531/1‑149 (MQ=255)
                                   cTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTa                                                                                                        <  2:651064/149‑1 (MQ=255)
                                         aGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCAcc                                                                                                  >  2:868568/1‑149 (MQ=255)
                                         aGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATAATCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGACTTTACACTGATGTGATTTACTTCACATCTTTTTACGTTGTACTCAcc                                                                                                  >  2:821757/1‑149 (MQ=255)
                                             aaGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATc                                                                                              <  1:658976/149‑1 (MQ=255)
                                              aGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCt                                                                                             >  2:621065/1‑149 (MQ=255)
                                                 cTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTaa                                                                                          <  2:691202/149‑1 (MQ=255)
                                                 cTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTaa                                                                                          <  1:880278/149‑1 (MQ=255)
                                                 cTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTaa                                                                                          <  1:880262/149‑1 (MQ=255)
                                                  ttAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACGCCACTATTTTatt                                                                                         >  2:199469/1‑147 (MQ=255)
                                                      ttGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTcac                                                                                     >  2:861752/1‑149 (MQ=255)
                                                      ttGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCAATTAAATCATAACTAAGATCAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCTTCCTAACCTATCTTAGTCcac                                                                                     >  2:634259/1‑149 (MQ=255)
                                                                  gCGTTTAATCTATGATGATATAACTCAATTACTTTCATGCACTTAAATCATAGCTAACATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTTCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATaa                                                                         <  1:824212/149‑1 (MQ=255)
                                                                            tATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGa                                                               <  2:355837/149‑1 (MQ=255)
                                                                             atgatgATATAACTCAATTATTTTTATGCATTTAAATCATAACTCATATAAATGTTAGTGTAAGAGATTACACTGATGTGACTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGTAg                                                              <  2:433283/149‑1 (MQ=255)
                                                                               gatgatATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCt                                                                                             <  2:951715/116‑1 (MQ=255)
                                                                               gatgatATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCt                                                                                             >  1:951715/1‑116 (MQ=255)
                                                                               gatgatATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAggg                                                            >  2:595897/1‑149 (MQ=255)
                                                                                 tgatATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGca                                                          <  1:597040/149‑1 (MQ=255)
                                                                                    tataACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAACAAATAACCATATTGGAGGGcctca                                                       >  1:515009/1‑149 (MQ=255)
                                                                                     ataACTCAATTATTTTCATGCACTTAAATCATGACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGcatcat                                                      <  1:881229/149‑1 (MQ=255)
                                                                                       aaCTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGc                                                    <  2:443066/149‑1 (MQ=255)
                                                                                        aCTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGcc                                                   >  2:99327/1‑149 (MQ=255)
                                                                                              ttattTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGCCGc                                             >  1:73693/1‑149 (MQ=255)
                                                                                               tattTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCCCATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATTAAAAATAACCATATTGGAGGTCATCATACCTGACGCt                                            >  2:937694/1‑149 (MQ=255)
                                                                                                 ttttCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTaa                                          <  2:790939/149‑1 (MQ=255)
                                                                                                       tGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTcac                                                                                     >  2:455238/1‑100 (MQ=255)
                                                                                                       tGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTcac                                                                                     <  1:455238/100‑1 (MQ=255)
                                                                                                           cTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACAATATTGGAGGGCGTCATGCCTGACGCTAAAAAACAgggg                                >  1:596515/1‑149 (MQ=255)
                                                                                                            ttAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGc                               >  1:777116/1‑149 (MQ=255)
                                                                                                            ttAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGc                               >  1:777092/1‑149 (MQ=255)
                                                                                                                       cTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAg                    <  1:468724/149‑1 (MQ=255)
                                                                                                                        tAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAgg                   <  1:621065/149‑1 (MQ=255)
                                                                                                                           gATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCaa                >  2:39821/1‑149 (MQ=255)
                                                                                                                            aTAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTa                                                                                                        >  2:395501/1‑60 (MQ=255)
                                                                                                                            aTAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTa                                                                                                        <  1:395501/60‑1 (MQ=255)
                                                                                                                               aaTGTTAGTTTAAGTGATTACACTGATGTGATTTGATTAACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGAc            <  1:699229/149‑1 (MQ=255)
                                                                                                                                aTGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACg           <  1:68230/149‑1 (MQ=255)
                                                                                                                                aTGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACg           <  1:170048/149‑1 (MQ=255)
                                                                                                                                aTGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACg           >  2:738140/1‑149 (MQ=255)
                                                                                                                                       tgtAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCg    >  1:444625/1‑149 (MQ=255)
                                                                                                                                         tAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTc  <  1:738140/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                
CATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTC  >  NC_000913/3088052‑3088337

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