Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,712,518 (C)6→7 intergenic (+360/‑251) nth → / → dtpA endonuclease III/dipeptide/tripeptide:H(+) symporter DtpA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,712,5121.C88.5% 79.0 / 4.3 26intergenic (+354/‑257)nth/dtpAendonuclease III/dipeptide/tripeptide:H(+) symporter DtpA
Reads supporting (aligned to +/- strand):  ref base . (0/3);  new base C (13/10);  total (13/13)
Fisher's exact test for biased strand distribution p-value = 2.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCAT‑CA‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGACTTGT  >  NC_000913/1712376‑1712641
                                                                                                                                          |                                                                                                                                 
aTCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCAT‑CA‑ccc                                                                                                                                <  1:22551/140‑1 (MQ=255)
  cAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCAT‑CA‑ccccc                                                                                                                              <  1:187277/140‑1 (MQ=255)
  cAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCAT‑CA‑ccccc                                                                                                                              <  1:189712/140‑1 (MQ=255)
       tAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCAT‑CA‑ccccccctct                                                                                                                         <  1:54242/140‑5 (MQ=255)
       tAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCAT‑CA‑ccccccctct                                                                                                                         <  2:393653/140‑5 (MQ=255)
               tttAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCAT‑CA‑ccccc                                                                                                                              <  2:26971/127‑1 (MQ=255)
               tttAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCAT‑CA‑ccccc                                                                                                                              >  1:26971/1‑127 (MQ=255)
                                                      aTTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCAT‑CACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACg                                                                           <  2:177494/140‑1 (MQ=255)
                                                      aTTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCAT‑CACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACg                                                                           <  1:44938/140‑1 (MQ=255)
                                                          cacTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCAT‑CACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTgg                                                                       >  2:166237/1‑140 (MQ=255)
                                                               gCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCAT‑CACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACg                                                                           >  1:287767/1‑131 (MQ=255)
                                                               gCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCAT‑CACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACg                                                                           <  2:287767/131‑1 (MQ=255)
                                                                         aGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCAT‑CACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGa                                                        >  2:4368/1‑140 (MQ=255)
                                                                                     aTCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCAT‑CACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTc                                            >  1:4493/1‑140 (MQ=255)
                                                                                      tctcACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCAT‑CACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAAtt                                             >  1:308771/1‑138 (MQ=255)
                                                                                      tctcACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCAT‑CACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAAtt                                             <  2:308771/138‑1 (MQ=255)
                                                                                       ctcACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCAT‑CACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTccc                                          <  1:35815/140‑1 (MQ=255)
                                                                                          aCTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCAT‑CACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGc                                       <  1:166237/140‑1 (MQ=255)
                                                                                               tAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCATCCC‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAg                                  <  2:358437/140‑1 (MQ=255)
                                                                                                aaaaaCGCTAAATAGCAGAACATATCGCCGTTAAGCCAT‑CACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGt                                 >  1:252656/1‑140 (MQ=255)
                                                                                                 aaacCGCTAAATAGCAGAACATATCGCCGTTAAGCCATCCC‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGtt                                <  1:61512/140‑1 (MQ=255)
                                                                                                                 aaaCATATCGCCGTTAAGCCATCCC‑CCCCCCCCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGattat                <  1:96405/139‑1 (MQ=255)
                                                                                                                    cATATCGCCGTTAAGCCAT‑CACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTAc             >  2:8752/1‑140 (MQ=255)
                                                                                                                     atatCGCCGTTAAGCCAT‑CACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACa            <  2:252656/140‑1 (MQ=255)
                                                                                                                       atCGCCGTTAAGCCAT‑CACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTaa                        >  1:82306/1‑126 (MQ=255)
                                                                                                                       atCGCCGTTAAGCCAT‑CACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTaa                        <  2:82306/126‑1 (MQ=255)
                                                                                                                       atCGCCGTTAAGCCAT‑CACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACaaa          >  2:221874/1‑140 (MQ=255)
                                                                                                                       atCGCCGTTAAGCCAT‑CACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACaaa          >  1:12105/1‑140 (MQ=255)
                                                                                                                         cGCCGTTAAGCCAT‑CACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAgg        >  2:165915/1‑140 (MQ=255)
                                                                                                                           ccGTTAAGCCAT‑CACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGAc      <  2:12105/140‑1 (MQ=255)
                                                                                                                           ccGTTAAGCCAT‑CACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGAc      <  2:4493/140‑1 (MQ=255)
                                                                                                                            cGTTAAGCCAT‑CACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGACt     >  1:291496/1‑140 (MQ=255)
                                                                                                                               tAAGCCAT‑CACCCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGACTTGt  >  1:270626/1‑140 (MQ=255)
                                                                                                                                          |                                                                                                                                 
ATCAGGTTAAAATCATTTAAATTTACACACGCAACAAATATTGACCTACAAAACATTACACTGGCTATTTTTCAGAAACTGGACTATCTCACTAGTAAAAACGCTAAATAGCAGAACATATCGCCGTTAAGCCAT‑CA‑CCCCCCTCTCAGTGCAGTGAAAAAATCTGCCGTTACGTTTTTTGAAAAATTTAACGCTGGATAACATTTCCCGGAATGGTTGAATTCCCCGCCTCAGTTATATGTAACAGATTATTACAAAGGACTTGT  >  NC_000913/1712376‑1712641

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