Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 646,209 A→G R124R (CGT→CGC citT ← citrate:succinate antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913646,2090AG94.7% 115.5 / ‑1.7 38R124R (CGT→CGCcitTcitrate:succinate antiporter
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (17/19);  total (18/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.46e-01

ACCGGAAAAACCGTACCCCCGGTACGCGCGGTGTTGGAAGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGAATTTCACCAGGAAAAGGGCAATGCGACGACCTAACCCGGAAACTTCATACCCTAATGCAAAAATAAATGCGCCAAATACCAGCCATACCGTGGTGCTGGAAAAACCAGCCAGGCCCCATTTCAGCGCCTGTTTTTGCGCATTAAACGCTGGGTCAGCTAAT  >  NC_000913/646072‑646344
                                                                                                                                         |                                                                                                                                       
aCCGGAAAAACCGTACCCCCGGTACGCGCGGTGTTGGAAGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGAATTTCACCAGGAAAAGGGCAATgcggcg                                                                                                                                       >  1:74350/1‑140 (MQ=255)
 ccGGAAAAACCGTACCCCCGGTACGCGCGGTGTTGGAAGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGAATTTCACCAGGAAAAGGGCAATGcggcga                                                                                                                                      >  1:79785/1‑140 (MQ=255)
            gTACCCCCGGTACGCGCGGTGTTGGAAGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGAATTTCACCAGGAAAAGGGCAATGCGGCGACCTAACCCGGa                                                                                                                           <  2:358896/140‑1 (MQ=255)
            gTACCCCCGGTACGCGCGGTGTTGGAAGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGAATTTCACCAGGAAAAGGGCAATGCGGCGACCTAACCCGGa                                                                                                                           <  1:725748/140‑1 (MQ=255)
            gTACCCCCGGTACGCGCGGTGTTGGAAGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGAATTTCACCAGGAAAAGGGCAATGCGGCGACCTAACCCGGa                                                                                                                           <  1:714901/140‑1 (MQ=255)
              aCCCCCGGTACGCGCGGTGTTGGAAGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGAATTTCACCAGGAAAAGGGCAATGCGGCGACCTAACCCGGaaa                                                                                                                         <  1:193152/140‑1 (MQ=255)
                   cGGTACGCGCGGTGTTGGAAGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGAATTTCACCAGGAAAAGGGCAATGCGGCGACCTAACCCGGAAACTTCa                                                                                                                    <  1:470489/140‑1 (MQ=255)
                    ggTACGCGCGGTGTTGGAAGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGAATTTCACCAGGAAAAGGGCAATGCGACGACCTAACCCGGAAACTTCAt                                                                                                                   >  2:19379/1‑140 (MQ=255)
                     gTACGCGCGGTGTTGGAAGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGAATTTCACCAGGAAAAGGGCAATGCGGCGACCTAACCCGGAAACTTCATa                                                                                                                  >  1:529952/1‑140 (MQ=255)
                               tgtTGGAAGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGAATTTCACCAGGAAAAGGGCAATGCGGCGACCTAACCCGGAAACTTCATACCCTAATGCa                                                                                                        >  2:659295/1‑140 (MQ=255)
                                         tgtAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGAATTTCACCAGGAAAAGGGCAATGCGGCGACCTAACCCGGAAACTTCATACCCTAATGCAAAAATAAATg                                                                                              >  1:404892/1‑140 (MQ=255)
                                             aaCGGTGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGAATTTCACCAGGAAAAGGGCAATGCGGCGACCTAACCCGGAAACTTCATACCCTAATGCAAAAATAAATGCGcc                                                                                          <  1:16671/140‑1 (MQ=255)
                                             aaCGGTGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGAATTTCACCAGGAAAAGGGCAATGCGGCGACCTAACCCGGAAACTTCATACCCTAATGCAAAAATAAATGCGcc                                                                                          <  2:74350/140‑1 (MQ=255)
                                             aaCGGTGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGAATTTCACCAGGAAAAGGGCAATGCGGCGACCTAACCCGGAAACTTCATACCCTAATGCAAAAATAAATGCGcc                                                                                          <  1:472356/140‑1 (MQ=255)
                                                  tGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGAATTTCACCAGGAAAAGGGCAATGCGGCGACCTAACCCGGAAACTTCATACCCTAATGCAAAAATAAATGCGCCAAATa                                                                                     <  2:79785/140‑1 (MQ=255)
                                                   gCCAGCAGAATGTCGATAATGATAATCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGAATTTCACCAGGAAAAGGGCAATGCGGCGACCTAACCCGGAAACTTCATACCCTAATGCAAAAATAAATGCGCCAAATAc                                                                                    <  1:587492/140‑1 (MQ=255)
                                                   gCCAGCAGAATGTCGATAATGACAATCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGAATTTCACCAGGAAAAGGGCAATGCGGCGACCTAACCCGGAAACTTCATACCCTAATGCAAAAATAAATGCGCCAAATAc                                                                                    >  1:105049/1‑140 (MQ=255)
                                                      agcagAATGTCGATAATGACAATCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGAATTTCACCAGGAAAAGGGCAATGCGGCGACCTAACCCGGAAACTTCATACCCTAATGCAAAAATAAATGCGCCAAATACCAg                                                                                 <  2:529952/140‑1 (MQ=255)
                                                        cagAATGTCGATAATGACAATCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGAATTTCACCAGGAAAAGGGCAATGCGACGACCTAACCCGGAAACTTCATACCCTAATGCAAAAATAAATGCGCCAAATACCAGcc                                                                               <  1:19379/140‑1 (MQ=255)
                                                              gTCGATAATGACAATCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGAATTTCACCAGGAAAAGGGCAATGCGGCGACCTAACCCGGAAACTTCATACCCTAATGCAAAAATAAATGCGCCAAATACCAGCCATACCg                                                                         >  1:229881/1‑140 (MQ=255)
                                                                   tAATGACAATCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGAATTTCACCAGGAAAAGGGCAATGCGGCGACCTAACCCGGAAACTTCATACCCTAATGCAAAAATAAATGCGCCAAATACCAGCCATACCgtggtg                                                                    >  1:772747/1‑140 (MQ=255)
                                                                   tAATGACAATCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGAATTTCACCAGGAAAAGGGCAATGCGGCGACCTAACCCGGAAACTTCATACCCTAATGCAAAAATAAATGCGCCAAATACCAGCCATACCgtggtg                                                                    <  2:404892/140‑1 (MQ=255)
                                                                   tAATGACAATCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGAATTTCACCAGGAAAAGGGCAATGCGGCGACCTAACCCGGAAACTTCATACCCTAATGCAAAAATAAATGCGCCAAATACCAGCCATACCgtggtg                                                                    >  1:437504/1‑140 (MQ=255)
                                                                     aTGACAATCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGAATTTCACCAGGAAAAGGGCAATGCGGCGACCTAACCCGGAAACTTCATACCCTAATGCAAAAATAAATGCGCCAAATACCAGCCATACCGTGGTGCt                                                                  >  2:799554/1‑140 (MQ=255)
                                                                            tCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGAATTTCACCAGGAAAAGGGCAATGCGGCGACCTAACCCGGAAACTTCATACCCTAATGCAAAAATAAATGCGCCAAATACCAGCCATACCGTGGTGCTGGaaaaa                                                           <  2:772747/140‑1 (MQ=255)
                                                                            tCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGAATTTCACCAGGAAAAGGGCAATGCGGCGACCTAACCCGGAAACTTCATACCCTAATGCAAAAATAAATGCGCCAAATACCAGCCATACCGTGGTGCTGGaaaaa                                                           <  2:437504/140‑1 (MQ=255)
                                                                                 tAACCCAACGTCAGCGTGCGTTTGCCCATGAATTTCACCAGGAAAAGGGCAATGCGGCGACCTAACCCGGAAACTTCATACCCTAATGCAAAAATAAATGCGCCAAATACCAGCCATACCGTGGTGCTGGAAAAAccagc                                                      >  1:597838/1‑140 (MQ=255)
                                                                                       aaCGTCAGCGTGCGTTTGCCCATGAATTTCACCAGGAAAAGGGCAATGCGGCGACCTAACCCGGAAACTTCATACCCTAATGCAAAAATAAATGCGCCAAATACCAGCAATACCGTGGTGCTGGAAAAACCAGCCAGGcc                                                >  1:90962/1‑140 (MQ=255)
                                                                                        aCGTCAGCGTGCGTTTGCCCATGAATTTCACCAGGAAAAGGGCAATGCGGCGACCTAACCCGGAAACTTCATACCCTAATGCAAAAATAAATGCGCCAAATACCAGCCATACCGTGGTGCTGGAAAAACCAGCCAGGccc                                               >  2:566658/1‑140 (MQ=255)
                                                                                              gcgtgcgtTTGCCCATGAATTTCACCAGGAAAAGGGCAATGCGGCGACCTAACCCGGAAACTTCATACCCTAATGCAAAAATAAATGCGCCAAATACCAGCCATACCGTGGTGCTGGAAAAACCAGCCAGGCCCCATTTc                                         >  1:177998/1‑140 (MQ=255)
                                                                                               cgtgcgtTTGCCCATGAATTTCACCAGGAAAAGGGCAATGCGGCGACCTAACCCGGAAACTTCATACCCTAATGCAAAAATAAATGCGCCAAATACCAGCCATACCGTGGTGCTGGAAAAAccagcca                                                    <  1:304194/128‑1 (MQ=255)
                                                                                               cgtgcgtTTGCCCATGAATTTCACCAGGAAAAGGGCAATGCGGCGACCTAACCCGGAAACTTCATACCCTAATGCAAAAATAAATGCGCCAAATACCAGCCATACCGTGGTGCTGGAAAAAccagcca                                                    >  2:304194/1‑128 (MQ=255)
                                                                                                        gCCCATGAATTTCACCAGGAAAAGGGCAATGCGGCGACCTAACCCGGAAACTTCATACCCTAATGCAAAAATAAATGCGCCAAATACCAGCCATACCGTGGTGCTGGAAAAACCAGCCAGGCCCCATTTCAGCGCCTGtt                               >  2:793702/1‑140 (MQ=255)
                                                                                                           cATGAATTTCACCAGGAAAAGGGCAATGCGGCGACCTAACCCGGAAACTTCATACCCTAATGCAAAAATAAATGCGCCAAATACCAGCCATACCGTGGTGCTGGAAAAACCAGCCAGGCCCCATTTCAGCGCCTGttttt                            <  1:8804/140‑1 (MQ=255)
                                                                                                             tGAATTTCACCAGGAAAAGGGCAATGCGGCGACCTAACCCGGAAACTTCATACCCTAATGCAAAAATAAATGCGCCAAATACCAGCCATACCGTGGTGCTGGAAAAACCAGCCAGGCCCCATTTCAGCGCCTGTTTTTgc                          <  2:597838/140‑1 (MQ=255)
                                                                                                             tGAATTTCACCAGGAAAAGGGCAATGCGGCGACCTAACCCGGAAACTTCATACCCTAATGCAAAAATAAATGCGCCAAATACCAGCCATACCGTGGTGCTGGAAAAACCAGCCAGGCCCCATTTCAGCGCCTGTTTTTgc                          >  1:815625/1‑140 (MQ=255)
                                                                                                                            aaaGGGCAATGCGGCGACCTAACCCGGAAACTTCATACCCTAATGCAAAAATAAATGCGCCAAATACCAGCCATACCGTGGTGCTGGAAAAACCAGCCAGGCCCCATTTCAGCGCCTGTTTTTGCGCATTAAACGCTggg           <  2:429838/140‑1 (MQ=255)
                                                                                                                                  cAATGCGGCGACCTAACCCGGAAACTTCATACCCTAATGCAAAAATAAATGCGCCAAATACCAGCCATACCGTGGTGCTGGAAAAACCAGCCAGGCCCCATTTCAGCGCCTGTTTTTGCGCATTAAACGCTGGGTCAGCt     <  1:566658/140‑1 (MQ=255)
                                                                                                                                     tGCGGCGACCTAACCCGGAAACTTCATACCCTAATGCAAAAATAAATGCGCCAAATACCAGCCATACCGTGGTGCTGGAAAAACCAGCCAGGCCCCATTTCAGCGCCTGTTTTTGCGCATTAAACGCTGGGTCAGCTAAt  >  2:410306/1‑140 (MQ=255)
                                                                                                                                     tGCGGCGACCTAACCCGGAAACTTCATACCCTAATGCAAAAATAAATGCGCCAAATACCAGCCATACCGTGGTGCTGGAAAAACCAGCCAGGCCCCATTTCAGCGCCTGTTTTTGCGCATTAAACGCTGGGTCAGCTAAt  >  1:473420/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                       
ACCGGAAAAACCGTACCCCCGGTACGCGCGGTGTTGGAAGGTGTAAACGGTGCCAGCAGAATGTCGATAATGACAATCGCATAACCCAACGTCAGCGTGCGTTTGCCCATGAATTTCACCAGGAAAAGGGCAATGCGACGACCTAACCCGGAAACTTCATACCCTAATGCAAAAATAAATGCGCCAAATACCAGCCATACCGTGGTGCTGGAAAAACCAGCCAGGCCCCATTTCAGCGCCTGTTTTTGCGCATTAAACGCTGGGTCAGCTAAT  >  NC_000913/646072‑646344

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