Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,154,691 G→A Q53* (CAG→TAG)  argE ← acetylornithine deacetylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,154,6910GA92.6% 74.2 / ‑0.7 27Q53* (CAG→TAG) argEacetylornithine deacetylase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (9/16);  total (10/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.26e-01

GGTTTTTTCAGTTTCGTGACGTCGACATCGCGTAGCGCATCAAGGATAAACGCAAAAAAGCCTTTCATGTCGGCGGTGCCTAAGCCGTAAAGCTTGCCGTCATGCTCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCGTCATCAAATGGCACCGTATCGGT‑ATGCCCCGCCAGCAACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCAC  >  NC_000913/4154444‑4154936
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     
ggTTTTTTCAGTTTCGTGACGTCGACAACGCGTAGCGCATCAAGAATAAACGCAAAAAAGCCTTTCATGTCGGCGGTGCCTAAGCCGTAAAGCTTGCCGTCATGCTCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCGTCATCAAATGGCACCGTATCGGT‑ATGCCCCGCCAGCAACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTacac                                                                                                                                                                                                                                                    <  1:30922/251‑1 (MQ=255)
        cAGTTTCGTGACGTCGACATCGCGTAGCGCATCAAGGATAAACGCAAAAAAGCCTTTCATGTCGGCGGTGCCTAAGCCGTAAAGCTTGCCGTCATGCCCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCGTCATCAAATGGCACCGTATCGGT‑ATGCCCCGCCAGCAACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTACACTTCcaca                                                                                                                                                                                                                                             <  1:142632/250‑1 (MQ=255)
                   cGTCGACATCGCGTAGCGCATCAAGGATAAACGCAAAAAAGCCTTTCATGTCGGCGGTGCCTAAGCCGTAAAGCTTGCCGTCATGCTCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCGTCATCAAATGGCACCGTATCGGT‑ATGCCCCGCCAGCAACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCaa                                                                                                                                                                                                                                  >  1:191870/1‑250 (MQ=255)
                      cGACATCGCGTAGCGCATCAAGGATAAACGCAAAAAAGCCTTTCATGTCGGCGGTGCCTAAGCCGTAAAGCTTGCCGTCATGCTCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCGTCATCAAATGGCACCGTATCGGT‑ATGCCCCGCCAGCAACAAGCCACCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTACACTTCCACATTGAAGCCCAActct                                                                                                                                                                                                                              >  1:126893/1‑251 (MQ=255)
                         cATCGCGTAGCGCATCAAGGATAAACGCAAAAAAGCCTTTCATGTCGGCGGTGCCTAAGCCGTAAAGCTTGCCGTCATGCTCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCGTCATCAAATGGCACCGTATCGGT‑ATGCCCCGCCAGCAACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTACACTTCCACATTGAAGCCCAAATCTTTa                                                                                                                                                                                                                           <  1:149606/251‑1 (MQ=255)
                                      cacAAGGATAACCCAAAAAAAGCTTTGCATGTCCGCGGGGCCTAAGCCGAAAAGCTTGCCGTCATGCTCCGTCAGTGTAAACGATTCGCCCCTCCAGGGCCCGTCATCAAATGGCCCCGTATCGTT‑ATGCCCCGCCAGCAACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTACACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCca                                                                                                                                                                                                              <  2:199702/249‑1 (MQ=255)
                                      aacAAGGATAAACACAAAAAGGCCTTTCAGGTCGGCGGTGCCTAAGCCTTAAAGCTTGCCGTCACGCTCCGCCAGTGTAAACGGAGCGCGCGTCCAGCGACCGTCATCAAATGGCACCGTACCGGT‑ATGCCCCGCCAGCAACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATAGAATTTGTTGCGAGTTCCTGGAACAGGCTACACTTCCTCATTGAAGCCCAAATCTTTAAACCAGTACGCca                                                                                                                                                                                                              <  1:201901/249‑1 (MQ=255)
                                       tCAAGGATAAACGCAAAAAAGCCTTTCATGTCGGCGGTACCTAAGACGTAAAGCTTGACGTCATGCTCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCGGCATCAAATGGCACCGTAACGGC‑ATGCCCCGCCAGCAACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGGAGGGAGGTCCTGGAACCGGCTAAACTTCCCCATTGAGACCCAACACTTTAAACCAGtacacca                                                                                                                                                                                                              >  2:308/1‑244 (MQ=255)
                                             aaaaaCGCAAAAAAGCCTTCCATGTCGGCGGTGCCTAAGCCGTAAAGCTTGCCGTCATGCTCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCGCCACCAAATGGCCCCGTATCGGT‑ATGCCCCGCCAGCAACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTACACTTCcaca                                                                                                                                                                                                                                             <  2:92000/211‑1 (MQ=255)
                                             aTAAACGCAAAAAAGCCTTTCATGTCGGCGGTGCCTAAGCCGTAAAGCTTGCCGTCATGCTCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCGTCATCAAATGGCACCGTATCGGT‑ATGCCCCGCCAGCAACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTACACTTCcaca                                                                                                                                                                                                                                             >  1:92000/1‑213 (MQ=255)
                                               aaaCGCAAAAAAGCCTTTCATGTCGGCGGTGCCTAAGCCGTAAAGCTTGCCGTCATGCTCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCGTCATCAAATGGCACCGTATCGGT‑ATGCCCCGCCAGCCACAAGCCGCCAGACCCCTGTCCGCTACTTGCCAGCATATTGAATTTGTTTCGAGTTCCTGGAACAGGCTACCCTTCCACATTGAAGCCCAAATCTTTAAACCAGTCAGCCAGCagagcgt                                                                                                                                                                                                     >  2:102904/1‑248 (MQ=255)
                                                   gCAAAAAAGCCTTTCATGTCGGCGGTGCCTAAGCCGTAAAGCTTGCCGTCATGCTCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCGTCATCAAATGGCACCGTATCGGT‑ATGCCCCGCCAGCAACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTACACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTa                                                                                                                                                                                                  >  2:162015/1‑250 (MQ=255)
                                                               ttCATGTCGGCGGTGCCTAAGCCGTAAAGCTTGCCGTCATGCTCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCGTCATCAAATGGCACCGTATCGGT‑ATGCCCCGCCAGCAACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTACACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCt                                                                                                                                                                                     >  1:176820/1‑251 (MQ=255)
                                                                    gccggcgtTGCCTAAGCCGTAAAGCTTGCCGTCATGCTCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCTTCATCAAATGGCACCGTATCGGT‑ATGCCCCGCCAGCAACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTACACTTCCACATTGAAGCCCAAATCTTTAAACCAGGCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGAt                                                                                                                                                                                <  1:89595/249‑1 (MQ=255)
                                                                        gtgggtcttaacccttacgtcTTGCCGTCATGCTCCGTCAGTGTAAACGGACCGCGCGCCCAGCGACCGCAATCAACAGGCCCCGGA‑CGGTTATGCCCCGCCAGCAACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTACACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAg                                                                                                                                                                            <  2:91611/231‑1 (MQ=255)
                                                                          ggTGCCTAAGCCGTAAAGCTTGCCGTCATGCCCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCGTCATCAAATGGCACCGTATCGGT‑ATGCCCCGCCAGCAACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTACACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTg                                                                                                                                                                          <  1:262860/251‑1 (MQ=255)
                                                                                                                                    catcGACCGTCATCAAATGGCCCCGTATCGGT‑ATGCCCCGCCAGCAACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTACACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATaaa                                                                                                                <  1:135967/248‑1 (MQ=255)
                                                                                                                                                 aaaatggacccgattCGGT‑ATGCCCCGCCAGCAACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTACACTTCCACATTGTAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTg                                                                                                   <  2:292555/237‑1 (MQ=255)
                                                                                                                                                               ccGGT‑ATGCCCGGCCAGCAACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTACACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGa                                                                                     <  1:213248/250‑1 (MQ=255)
                                                                                                                                                                          ccGCCAGCAACAAGCCACCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTACACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTg                                                                           <  2:126893/251‑1 (MQ=255)
                                                                                                                                                                                gCAACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACCCACCTCTGGTCATg                                                                     <  2:191870/251‑1 (MQ=255)
                                                                                                                                                                                           ccAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTACACTTCCTCATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGAGTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTGCATTGTTGACACACCTCTTGTCATGATAGTATCAat                                                          <  2:97680/251‑1 (MQ=255)
                                                                                                                                                                                            cAGCCCCCCGTCAGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTACACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATTTGCATTGATTTGATCGAGTACCTCTTCCGTGGCGCTTATTGAAGGTGTGGAAATCAGAGCGCGGTACTTCTCGATAAATGGCGGTAATTTGGTTTTCATTGTTTCAACACCCCTGGGTATGATAGTATCAat                                                          >  2:179907/1‑250 (MQ=255)
                                                                                                                                                                                                ccccTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTACACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTc                                                      >  1:357468/1‑250 (MQ=255)
                                                                                                                                                                                                         gATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTACACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATAAATATTCCTGCAGTAtt                                            >  2:65383/1‑251 (MQ=255)
                                                                                                                                                                                                                                    ccGAGTTCCTGGAACATGCTACCCTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCTCCAGCAGAGTTATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGTCGCTTATTTAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGTCGGTAATTTTTTTTTCATTGTTGACACACCTCTGTTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACTTTGac                 <  2:167578/250‑2 (MQ=255)
                                                                                                                                                                                                                                              ggAACAGGCTACACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGtaataa         <  1:162015/249‑1 (MQ=255)
                                                                                                                                                                                                                                                   aGGCTACACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTGTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACccac  <  2:176820/251‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     
GGTTTTTTCAGTTTCGTGACGTCGACATCGCGTAGCGCATCAAGGATAAACGCAAAAAAGCCTTTCATGTCGGCGGTGCCTAAGCCGTAAAGCTTGCCGTCATGCTCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCGTCATCAAATGGCACCGTATCGGT‑ATGCCCCGCCAGCAACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGCTGCACTTCCACATTGAAGCCCAAATCTTTAAACCAGTCCGCCAGCAGAGTGATTAAATCTGCATTGCTTTGATCGAGTGCCTCTTCCGTGGCGCTTATTGAAGGTGTGGCAATCAGAGCGCGGTAAATCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCAC  >  NC_000913/4154444‑4154936

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