Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,416,945 C→T 8.5% intergenic (+54/‑136) pta → / → yfcC phosphate acetyltransferase/putative inner membrane transporter; C4‑dicarboxylate anaerobic carrier family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,416,9450CT 8.5% 139.9 / 3.9 47intergenic (+54/‑136)pta/yfcCphosphate acetyltransferase/putative inner membrane transporter; C4‑dicarboxylate anaerobic carrier family protein
Reads supporting (aligned to +/- strand):  ref base C (21/22);  new base T (2/2);  total (23/24)
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.26e-01

GTTAACGACCTGTCCCGTGGCGCACTGGTTGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCCCACTTTGCATACTTACCACTTTGTTTTGTGCAAGGGAATATTTGCGCTATGTCCG  >  NC_000913/2416805‑2417087
                                                                                                                                            |                                                                                                                                              
gTTAACGACCTGTCCCGTGGCGCACTGGTTGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCg                                                                                                                                        >  1:273853/1‑149 (MQ=255)
 ttAACGACCTGTCCCGTGGCGCACTGGTTGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAGTGATTATTTCCgg                                                                                                                                       <  3:1357968/149‑1 (MQ=255)
      gACCTGTCCCGTGGCGCACTGGTTGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGttttt                                                                                                                                  >  4:997390/1‑149 (MQ=255)
      gACCTGTCCCGTGGCGCACTGGTTGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGttttt                                                                                                                                  >  2:1161614/1‑149 (MQ=255)
         cTGTCCCGTGGCGCACTGGTTGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTAt                                                                                                                               <  4:762272/149‑1 (MQ=255)
               cGTGGCGCACTGGTTGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCtt                                                                                                                         <  1:1161614/149‑1 (MQ=255)
                   gcgcACTGGTTGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGt                                                                                                                                      >  2:552246/1‑132 (MQ=255)
                   gcgcACTGGTTGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGt                                                                                                                                      <  1:552246/132‑1 (MQ=255)
                            ttGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATtctc                                                                                                                           <  4:486007/134‑1 (MQ=255)
                            ttGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATtctc                                                                                                                           >  3:486007/1‑134 (MQ=255)
                                          tACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATc                                                                                              <  3:997390/149‑1 (MQ=255)
                                          tACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATc                                                                                              <  2:273853/149‑1 (MQ=255)
                                             acCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGa                                                                                                   >  1:199797/1‑141 (MQ=255)
                                             acCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGa                                                                                                   <  2:199797/141‑1 (MQ=255)
                                                aTCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTAtt                                                                                                                              <  1:1075742/111‑1 (MQ=255)
                                                aTCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTAtt                                                                                                                              >  2:1075742/1‑111 (MQ=255)
                                                           tGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGtttt                                                                                                                                   >  1:844829/1‑95 (MQ=255)
                                                           tGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGtttt                                                                                                                                   <  4:941304/95‑1 (MQ=255)
                                                           tGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGtttt                                                                                                                                   >  3:941304/1‑95 (MQ=255)
                                                           tGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGtttt                                                                                                                                   <  2:844829/95‑1 (MQ=255)
                                                              gATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCt                                                                                                     <  1:488389/122‑1 (MQ=255)
                                                              gATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCt                                                                                                     >  2:488389/1‑122 (MQ=255)
                                                                            agcagcagTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTg                                                                                                    >  1:1029245/1‑109 (MQ=255)
                                                                            agcagcagTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTg                                                                                                    <  2:1029245/109‑1 (MQ=255)
                                                                                                     gCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCtt                                                                                                        >  2:281182/1‑80 (MQ=255)
                                                                                                     gCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCtt                                                                                                        <  1:281182/80‑1 (MQ=255)
                                                                                                                  gCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATc                                                          >  1:45474/1‑113 (MQ=255)
                                                                                                                  gCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATc                                                          <  2:45474/113‑1 (MQ=255)
                                                                                                                   cGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCaa                                                                        >  2:402970/1‑98 (MQ=255)
                                                                                                                   cGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCaa                                                                        <  1:402970/98‑1 (MQ=255)
                                                                                                                     gATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCCCACTTTGCATACTTACCACTTTGTTTTGTGCAAGGGaa                   <  1:234609/149‑1 (MQ=255)
                                                                                                                         tCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCCCACTTTGCATACTTACCACtttgtttttt                            <  4:1002809/136‑3 (MQ=255)
                                                                                                                         tCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCCCACTTTGCATACTTACCACtttgttttgt                            >  3:1002809/1‑136 (MQ=255)
                                                                                                                         tCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCCCACTTTGCATACTTACCACtttgtttt                              <  2:641018/134‑1 (MQ=255)
                                                                                                                         tCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCCCACTTTGCATACTTACCACtttgtttt                              >  1:641018/1‑134 (MQ=255)
                                                                                                                          cTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCCCACTTTGCATACTTACCACtttgttttgt                            <  1:650068/135‑1 (MQ=255)
                                                                                                                          cTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCCCACTTTGCATACTTACCACtttgttttgt                            >  2:650068/1‑135 (MQ=255)
                                                                                                                           tGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCCCACTTTGCa                                                <  3:473564/114‑1 (MQ=255)
                                                                                                                           tGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCCCACTTTGCa                                                >  4:473564/1‑114 (MQ=255)
                                                                                                                             aaCCGGAAATAGTGATTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCCCACTTTGCATACTTACCACTTTGTTTTGTGCAAGGGAATATTTgcg           >  4:1204080/1‑149 (MQ=255)
                                                                                                                                cGGAAATAGTGATTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGc                                                                                          <  1:167696/67‑1 (MQ=38)
                                                                                                                                cGGAAATAGTGATTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGc                                                                                          >  2:167696/1‑67 (MQ=38)
                                                                                                                                      taatCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCCCACTTTGCATACTTACCACTTTGTTTTGTGCAAGGGAATATTTGCGCTATGTCCg  >  1:298563/1‑149 (MQ=255)
                                                                                                                                          cACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCCCACTTTGCATACTTACCACTTTGTTTTGTGCAAGGGAATATTTgc            >  3:724971/1‑135 (MQ=255)
                                                                                                                                          cACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCCCACTTTGCATACTTACCACTTTGTTTTGTGCAAGGGAATATTTgc            <  4:724971/135‑1 (MQ=255)
                                                                                                                                           aCTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCCCACTTTGCATCCTTAc                                         <  3:370102/105‑1 (MQ=255)
                                                                                                                                           aCTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCCCACTTTGCATACTTAc                                         >  4:370102/1‑105 (MQ=255)
                                                                                                                                            |                                                                                                                                              
GTTAACGACCTGTCCCGTGGCGCACTGGTTGACGATATCGTCTACACCATCGCGCTGACTGCGATTCAGTCTGCACAGCAGCAGTAATCTCGTCATCATCCGCAGCTTTGCGCTGCGGATATCTGAACCGGAAATAATCACTATTTCCGGTTTTTTATTCTCTTAATTTGCATTAATCCTTTCTGATTATCTTGCTTAACTGCGCTGCATCAATGAATTGCGCCATCCCACTTTGCATACTTACCACTTTGTTTTGTGCAAGGGAATATTTGCGCTATGTCCG  >  NC_000913/2416805‑2417087

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