Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 645,062 G→A 11.8% intergenic (‑59/+55) rna ← / ← citT ribonuclease I/citrate/succinate antiporter; citrate carrier

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913645,0620GA11.8% 89.0 / 5.0 34intergenic (‑59/+55)rna/citTribonuclease I/citrate/succinate antiporter; citrate carrier
Reads supporting (aligned to +/- strand):  ref base G (13/17);  new base A (2/2);  total (15/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.81e-01

CTGTTTTGCCTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAGTAATCTTTT  >  NC_000913/644914‑645210
                                                                                                                                                    |                                                                                                                                                    
cTGTTTTGCCTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCg                                                                                                                                                      <  1:82062/149‑1 (MQ=255)
   ttttGCCTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCa                                                                                                                                                   >  4:250835/1‑149 (MQ=255)
         cTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTggg                                                                                                                                             <  3:1098104/149‑1 (MQ=255)
            cAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGtaagccc                                                                                                                                          <  2:1349081/149‑7 (MQ=255)
              aCGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTAc                                                                                                                                        <  4:778581/149‑1 (MQ=255)
               cGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTAcc                                                                                                                                       <  2:59/149‑1 (MQ=255)
                gCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAg                                                                                                                                                  >  4:684763/1‑137 (MQ=255)
                gCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAg                                                                                                                                                  <  3:684763/137‑1 (MQ=255)
                gCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCt                                                                                                                                      <  4:1346209/149‑1 (MQ=255)
                 cTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTg                                                                                                                                     <  3:250835/149‑1 (MQ=255)
                                       gggAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAg                                                                                                               >  2:299816/1‑149 (MQ=255)
                                              agagAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTgg                                                                                                        <  1:299816/149‑1 (MQ=255)
                                               gagaAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGt                                                                                                                                                <  2:227002/108‑1 (MQ=255)
                                               gagaAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGt                                                                                                                                                >  1:227002/1‑108 (MQ=255)
                                                       gcgAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACCCTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGCgg                                                                                                                                              <  2:956568/102‑1 (MQ=255)
                                                       gcgAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTgg                                                                                                                                              >  1:956568/1‑102 (MQ=255)
                                                                  ggcgTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTg                                                                                                                                               >  4:1150057/1‑90 (MQ=255)
                                                                  ggcgTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTg                                                                                                                                               <  3:1150057/90‑1 (MQ=255)
                                                                             aGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCtt                                                                                                                                          <  3:636233/84‑1 (MQ=255)
                                                                             aGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCtt                                                                                                                                          >  4:636233/1‑84 (MQ=255)
                                                                                    tttCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCACCTGGCGGGCttttt                                                                                                                                       >  4:482700/1‑77 (MQ=38)
                                                                                    tttCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCACCTGGCGGGCttttt                                                                                                                                       <  3:482700/80‑4 (MQ=38)
                                                                                    tttCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCACCTGGCGGGCtttt                                                                                                                                        >  1:224582/1‑77 (MQ=38)
                                                                                    tttCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCACCTGGCGGGCtttt                                                                                                                                        <  2:224582/79‑3 (MQ=38)
                                                                                        aTAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGtgatga                                                                                                                        >  4:29930/1‑91 (MQ=255)
                                                                                        aTAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGtgatga                                                                                                                        <  3:29930/91‑1 (MQ=255)
                                                                                                      aCTCGTACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATaga                                                >  4:320732/1‑149 (MQ=255)
                                                                                                                                 ttAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCga                                                                               >  2:526368/1‑91 (MQ=255)
                                                                                                                                 ttAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCga                                                                               <  1:526368/91‑1 (MQ=255)
                                                                                                                                      tCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGAc                <  3:320732/149‑1 (MQ=255)
                                                                                                                                       cATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACg               >  4:950776/1‑149 (MQ=255)
                                                                                                                                                  ccGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAGTaa        >  2:834313/1‑145 (MQ=255)
                                                                                                                                                  ccGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAGTaa        <  1:834313/145‑1 (MQ=255)
                                                                                                                                                  ccGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAGTAATCtt    >  4:972531/1‑149 (MQ=255)
                                                                                                                                                   cGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAGTAATCttt   <  3:972531/149‑1 (MQ=255)
                                                                                                                                                    gCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGa                 <  3:810988/134‑1 (MQ=255)
                                                                                                                                                    gCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGa                 >  4:810988/1‑134 (MQ=255)
                                                                                                                                                    gCCAGGTAAGCCCACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAGTAATCtttt  >  3:1439572/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTGTTTTGCCTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGGGTGTGGAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCAACCAACAGCAACATAGAGATGTAAATCACCCCGAAGATTGCGCCAAGACGCCAGTAATCTTTT  >  NC_000913/644914‑645210

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