Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,518,660 T→C T63A (ACG→GCG)  mnaT ← L‑amino acid N‑acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,518,6600TC97.2% 112.9 / ‑4.5 37T63A (ACG→GCG) mnaTL‑amino acid N‑acyltransferase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (25/10);  total (27/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

TCGGCTTAACAATTTACGCCCCAGACCTTTGCCCTGATGATCGGGATGGACATAAACCGAATGTTCCACGGTATGGCGAAAACCATCGAAACTACGCCAGTCGCCAAACGAGGCATATCCCGTCACTACGCCGTTTTCCTCGCTCACCAGCACTGGATAACCTGCTAAAGTCCGCGCTTCAAACCAGGCAATGCGGTTATCAGCATCCACCGTTTGGTCATTCCAGATAGCCGCCGTATACAACACGGCGT  >  NC_000913/1518538‑1518788
                                                                                                                           |                                                                                                                                
tcGGCTTAACAATTTACGCCCCAGCCCTTTGCCCTGATGAGCGGGAGGGACATAAACCGAATGTCCCACGGTATGGCGAAACCCATCGAAACTCCGCCAGTCGCCAAACGAGGCATATCCCGCCACTACGCCGTTTTCCt                                                                                                                 <  2:512981/140‑1 (MQ=255)
tcGGCTTAACAATTTACGCCCCAGACCTTTGCCCTGATGATCGGGATGGACATAAACCGAATGTTCCACGGTATGGCGAAAACCATCGAAATTACGCCAGCCGCTAAACGAGGCATACCCCGGCACTACGCCGTTTTCCt                                                                                                                 >  2:374823/1‑140 (MQ=255)
    cTTAACAATTTACGCCCCAGACCTTTGCCCTGATGATCGGGATGGACATAAACCGAATGTTCCACGGTATGGCGAAAACCATCGAAACTACGCCAGTCGCCAGACGAGGCGTATCCCGCCACTACGCCGTTTTCCTCGCt                                                                                                             >  1:299728/1‑140 (MQ=255)
       aaCAATTTACGCCCCAGACCTTTGCCCTGATGATCGGGATGGACATAAACCGAATGTTCCACGGTATGGCGAAAACCATCGAAACTACGCCAGTCGCCAAACGAGGCATATCCCGCCACTACGCCGTTTTCCTCGCTCAc                                                                                                          >  1:11877/1‑140 (MQ=255)
        aCAATTTACGCCCCAGACCTTTGCCCTGATGATCGGGATGGACATAAACCGAATGTTCCACGGTATGGCGAAAACCATCGAAACTACGCCAGTCGCCAAACGTGGCATATACCGCCAATACGCCGTTTTCCTCGCTCAcc                                                                                                         >  1:282253/1‑140 (MQ=255)
                   cccAGCCCTTTCCCCTGATGATCGGGATGGAAATAAACCGAATGTTCCAGGGTAGGGCGAAAACCATCGAAACTCCGCCGGCCGCAAAACGAGGCATATCCCGCCACTACGCCGTTTTCCTCGCTCACCAGCACTGGATa                                                                                              <  2:783841/140‑1 (MQ=255)
                             tGCCCTGATGATCGGGATGGACATAAACCGAATGTTCCACGGTATGGCGAAAACCATCGAAACTACGCCAGTCGCCAAACGAGGCATATCCCGTCACTACGCCGTTTTCCTCGCTCACCAGCACTGGATAACCTGCTaaa                                                                                    >  2:766427/1‑140 (MQ=255)
                             tGCCCTGATGATCGGGATGGACATAAACCGAATGTTCCACGGTATGGCGAAAACCATCGAAACTACGCCAGTCGCCAAACGAGGCATATCCCGCCACTACGCCGTTTTCCTCGCTCACCAGCACTGGATAACCTGCTaaa                                                                                    >  1:155840/1‑140 (MQ=255)
                             tGCCCTGATGATCGGGATGGACATAAACCGAATGTTCCACGGTATGGCGAAAACCATCGAAACTACGCCAGTCGCCAAACGAGGCATATCCCGCCACTACGCCGTTTTCCTCGCTCACCAGCACTGGATAACCTGCTaaa                                                                                    >  1:682297/1‑140 (MQ=255)
                             tGCCCTGATGATCGGGATGGACATAAACCGAATGTTCCACGGTATGGCGAAAACCATCGAAACTACGCCAGTCGCCAAACGAGGCATATCCCGCCACTACGCCGTTTTCCTCGCTCACCAGCACTGGATAACCTGCTaaa                                                                                    >  1:529928/1‑140 (MQ=255)
                             tGCCCTGATGATCGGGATGGACATAAACCGAATGTTCCACGGTATGGCGAAAACCATCGAAACTACGCCAGTCGCCAAACGAGGCATATCCCGCCACTACGCCGTTTTCCTCGCTCACCAGCACTGGATAACCTGCTaaa                                                                                    >  1:449319/1‑140 (MQ=255)
                                    atgatCGGGATGGACATAAACCGAATGTTCCACGGTATGGCGAAAACCATCGAAACTACGCCAGTCGCCAAACGAGGCATATCCCGCCACTACGCCGTTTTCGTCGCTCACCAGCACTGGATAATCTGCTAAAGTCcgcg                                                                             >  1:318998/1‑140 (MQ=255)
                                    atgatCGGGATGGACATAAACCGAATGTTCCACGGTATGGCGAAAACCATCGAAACTACGCCAGTCGCCAAACGAGGCATATCCCGCCACTACGCCGTTTTCGTCGCTCACCAGCACAGGATAACCTGCTAAAGTCcgcg                                                                             >  1:351664/1‑140 (MQ=255)
                                    atgatCGGGATGGACATAAACCGAATGTTCCACGGTATGGCGAAAACCATCGAAACTACGCCAGTCGCCAAACGAGGCATATCCCGCCACTACGCCGTTTTCCTCGCTCACCAGCACTGGATAACCTGCTAAAGTCcgcg                                                                             >  2:61774/1‑140 (MQ=255)
                                    atgatCGGGATGGACATAAACCGAATGTTCCACGGTATGGCGAAAACCATCGAAACTACGCCAGTCGCCAAACGAGGCATATCCCGCCACTACGCCGTTTTCCTCGCTCACCAGCACTGGATAACCTGCTAAAGTCcgcg                                                                             >  2:485642/1‑140 (MQ=255)
                                         ggggATGAACATAACCCGAATTTCCCACGGTATGGCAAAACCCACCAAAACTCCGCCAGTCCCAAACCGAGGCATATCCCGCCACTACGCCGTTTTCCTCGCTCACCAGCACTGGATAACCTGCTAAAGTCCGCGCTTCa                                                                        <  2:155840/139‑1 (MQ=255)
                                              tGGACATAAACCGAATGTTCCACGGTATGGCGAAAACCATCGAAACTACGCCAGTCGCCAAACGAGGCATATCCCGCCACTACGCCGTTTTTCTCGCTCAACAGGACTAGATAACCTGCTAAAGTCCGCGCTTCAAACCa                                                                   >  2:265094/1‑140 (MQ=255)
                                              tGGACATAAACCGAATGTTCCACGGTATGGCGAAAACCATCGAAACTACGCCAGTCGCCAAACGAGGCATATCCCGCCACTACGCCGTTTTCCTCGCTCACCAGCACTGGATAACCTGCTAAAGTCCGCGCTTCAAACCa                                                                   >  1:477568/1‑140 (MQ=255)
                                              tGGACATAAACCGAATGTTCCACGGTATGGCGAAAACCATCGAAACTACGCCAGTCGCCAAACGAGGCATATCCCGCCACTACGCCGTTTTCCTCGCTCACCAGCACTGGATAACCTGCTAAAGTCCGCGCTTCAAACCa                                                                   >  2:571448/1‑140 (MQ=255)
                                              tGGACATAAACCGAATGTTCCACGGTATGGCGAAAACCATCGAAACTACGCCAGTCGCCAAACGAGGCATATCCCGCCACTACGCCGTTTTCCTCGCTCACCAGCACTGGATAACCTGCTAAAGTCCGCGCTTCAAACCa                                                                   >  2:480813/1‑140 (MQ=255)
                                                    aaaaCCGAATGTTCCCGGGTAGGGAAAAACCCATCGAACCTCCCCCGGTCGCAAAACGGGGCTTATCCCGCCACTACGCCGTTTTCCTCGCTCACCAGCACTGGATAACCTGCTAAAGTCCGCGCTTCAAACCAGGCAAt                                                             <  2:477568/139‑1 (MQ=255)
                                                     aacccGAATGTCCCACGGTATGGCAAAACCCATCAAAACTCCGCCGGTCGCCAAACGAGGCATATCCCGCCACTACGCCGTTTTCCTCGCTCACCAGCACTGGATAACCTGCTAAAGTCCGCGCTTCAAACCAGGCAATg                                                            <  2:557165/137‑1 (MQ=255)
                                                       aCCGAATGTTCCACGGTATGGCGAAAACCATCGAAACTACGCCAGTCGCCAAACGAGGCATATCCCGCCACTACGCCGTTTTCCTCGCTCACCAGCACTGGATAACCTGCTAAAGTCCGCGCTTCAAACCAGGCAATGCg                                                          >  2:175846/1‑140 (MQ=255)
                                                          gAATGTTCCACGGTATGGCGAAAACCATCGAAACTACGCCAGTCGCCAAACGAGGCATATCCCGCCACTACGCCGTTTTCCTCGCTCACCAGCACTGGATAACCTGCTAAAGTCCGCGCTTCAAACCAGGCAATGCGGtt                                                       >  1:177042/1‑140 (MQ=255)
                                                           aaTGTTCCACGGTATGGCGAAAACCATCGAAACTACGCCAGTCGCCAAACGAGGCATATCCCGCCACTACGCCGTTTTCCTCGCTCACCAGCACTGGATAACCTGCTAAAGTCCGCGCTTCAAACCAGGCAATGCGGTTa                                                      >  1:750434/1‑140 (MQ=255)
                                                           aaTGTTCCACGGTATGGCGAAAACCATCGAAACTACGCCAGTCGCCAAACGAGGCATATCCCGCCACTACGCCGTTTTCCTCGCTCACCAGCACTGGATAACCTGCTAAAGTCCGCGCTTCAAACCAGGCAATGCGGTTa                                                      >  1:557165/1‑140 (MQ=255)
                                                                   aCGGTATGGCGAAAACCATCGAAACTACGCCACTCGCCAAACGAGGCATATCCCGCCACTACGCCGTTTTCCTCGCTCACCAGCACTGGATAACCTGCTAAAGTCCGCGCTTCAAACCAGGCAATGCGGTTATCAGCATc                                              <  1:541360/140‑1 (MQ=255)
                                                                         tGGCGAAAACCATCGAAACTACGCCAGTCGCCAAACGAGGCATATCCCGCCACTACGCCGTTTTCCTCGCTCACCAGCACTGGATAACCTGCTAAAGTCCGCGCTTCAAACCAGGCAATGCGGTTATCAGCATCCACCGt                                        >  1:21907/1‑140 (MQ=255)
                                                                          ggCGAAAACCATCGAAACTACGCCAGTCGCCAAACGAGGCATATCCCGCCACTACGCCGTTTTCCTCGCTCACCAGCACTGGATAACCTGCTAAAGTCCGCGCTTCAAACCAGGCAATGCGGTTATCAGCATCCACCGtt                                       >  1:633302/1‑140 (MQ=255)
                                                                              aaacCCATCAAAACTCCGCCAGTCGCCAAACGAGGCATACCCCGCCACTACGCCGTTTTCCTCGCTCACCAGCACTGGATAACCTGCTAAAGTCCGCGCTTCAAACCAGGCAATGCGGTTATCAGCATCCACCGTTTGGt                                   <  2:682297/140‑1 (MQ=255)
                                                                                     tCGAAACTACGCCAGTCGCCAAACGAGGCATATCCCGCCACTACGCCGTTTTCCTCGCTCACCAGCACTGGATAACCTGCTAAAGTCCGCGTTTCAAACCTGGCGATGCGGTTATCACCATCGACCGTTTGGTCATTCCa                            >  2:380320/1‑140 (MQ=255)
                                                                                      cGAAACTACGCCAGTCGCCAAACGAGGCATATCCCGCCACTACGCCGTTTTCCTCGCTCACCAGCACTGGATAACCTGCTAAAGTCCGCGCGTCAAACCAGGCAGTGCAGTTATCAGAATCCGCCGTTTGGTCATTCCAg                           >  2:250191/1‑140 (MQ=255)
                                                                                             cgcccAGCCCCCAACGGAGCCATTCCCCGCCACTACGCCGTTTTCCTCGCTCACCAGCACTGGATAACCTGCTAAAGTCCGCGCTTCAAACCAGGCAATGCGGTTATCAGCATCCACCGTTTGGTCATTCCAGATAgccg                    <  1:524755/137‑1 (MQ=255)
                                                                                             aCGCCAGTCGCCAAACGAGGCATATCCCGCCACTACGCCGTTTTCCTCGCTCACCAGCACTGGATAACCTGCTAAAGTCCGCGCTTCAAACCAGGCAATGCGGTTATCAGTATCAACGGTTTGGGCATTCCAGATAgccg                    >  1:239835/1‑140 (MQ=255)
                                                                                               cccAGTCGCAAACGGAGGCATACCCCGCCACTACGCCGTTTTCCTCGCTCACCAGCACTGGATAACCTGCTAAAGTCCGCGCTTCAAACCAGGCAATGCGGTTATCAGCATCCACCGTTTGGTCATTCCAGATAgccgcc                  <  2:21907/139‑1 (MQ=255)
                                                                                                   ggccgcAAACGGAGGTATACCCTGCCACTGCGCCGTTGTCCTCGCTCACCAGCACTGGATAACCTGCTAAAGTCCGCGCTTCAAACCAGGCAATGCGGTTATCAGCATCCACCGTTTGGTCATTCCAGATAGCCGCCGta               <  1:374823/135‑1 (MQ=255)
                                                                                                   gTCGCCAAACGAGGCATATCCCGCCACTACGCCGTTTTCCTCGCTCACCAGCACTGGATAACCTGCTAAAGTCCGCGCTTCAAACCAGGCAATGCGGTTATCAGCATCCACCGTTTGGTCATTCCAGATAGCCGCCGtat              >  2:21449/1‑140 (MQ=255)
                                                                                                       caaaaCGGGGCATACCCCGCCACTACGCCGTTTTCCTCGCTCACCAGCACTGGATAACCTGCTAAAGTCCGCGCTTCAAACCAGGCAATGCGGTTATCAGCATCCACCGTTTGGTCATTCCAGATAGCCGCCGTATacaa          <  2:177042/138‑1 (MQ=255)
                                                                                                                taattctCCCGCCACGACGCCGATTGCCTTGCTCACCCGCACTGAATAACCTGCTAAAGTCCGCGCTTCAAACCAGGCAATGCGGTTATCAGCATCCACCGTTTGGTCATTCCAGATAGCCGCCGTATACAACACGGCGt  <  1:380320/134‑1 (MQ=255)
                                                                                                                           |                                                                                                                                
TCGGCTTAACAATTTACGCCCCAGACCTTTGCCCTGATGATCGGGATGGACATAAACCGAATGTTCCACGGTATGGCGAAAACCATCGAAACTACGCCAGTCGCCAAACGAGGCATATCCCGTCACTACGCCGTTTTCCTCGCTCACCAGCACTGGATAACCTGCTAAAGTCCGCGCTTCAAACCAGGCAATGCGGTTATCAGCATCCACCGTTTGGTCATTCCAGATAGCCGCCGTATACAACACGGCGT  >  NC_000913/1518538‑1518788

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