Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 30,615 C→T P322L (CCT→CTT)  carA → carbamoyl phosphate synthetase small subunit, glutamine amidotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091330,6150CT79.2% 39.8 / 9.6 24P322L (CCT→CTT) carAcarbamoyl phosphate synthetase small subunit, glutamine amidotransferase
Reads supporting (aligned to +/- strand):  ref base C (3/2);  new base T (11/8);  total (14/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 = 3.41e-01

GGCGAGCGGTGCGAAGACTGTCAAAATGAAATTTGGTCACCACGGCGGCAACCATCCGGTTAAAGATGTGGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCTGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGCACCCTGAAGCCAGCCCTGGTCCACACGACGCCGCGCCGTTGTTC  >  NC_000913/30481‑30751
                                                                                                                                      |                                                                                                                                        
ggcgAGCGGTGCGAAGACTGTCAAAATGAAATTTGGTCACCACGGCGGCAACCATCCGGTTAAAGATGTGGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACTTGCa                                                                                                                                      <  2:472354/139‑1 (MQ=255)
                      aaaaTGAAATTTGGTCACCACGGCGGCAACCATCCGGTTAAAGATGTGGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACTTGCAAACCTGCGTGTCACGCATAAAt                                                                                                                <  1:550955/139‑1 (MQ=255)
                              aTTTGGTCACCACGGCGGCAACCATCCGGTTAAAGATGTGGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACTTGCAAACCTGCGTGTCACGCATAAATCCCTGTTc                                                                                                        >  1:247967/1‑139 (MQ=255)
                                ttGGTCACCACGGCGGCAACCATCCGGTTAAAGATGTGGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACTTGCAAACCTGCGTGTCACACATAAATCCCTGTTCGa                                                                                                      <  2:115594/139‑1 (MQ=255)
                                     caccacGGCGGCAACCATCCGGTTAAAGATGTGGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCTGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTa                                                                                                 <  1:166913/139‑1 (MQ=255)
                                                              aaGATGTGGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACTTGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGa                                                                        >  1:319949/1‑139 (MQ=255)
                                                                         aaaaaCGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCTGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCAt                                                             >  2:105657/1‑139 (MQ=255)
                                                                           aaaCGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACTTGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTc                                                           >  1:110350/1‑139 (MQ=255)
                                                                                          cACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACTTGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGCAc                                            >  2:129373/1‑139 (MQ=255)
                                                                                            ccccccAGAACAACGGTTTTGCGGTGGACTAAGCAACATTACTTGCAAACCTGCGTGTCACGCATAAATTCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGCAccc                                          >  2:415932/4‑139 (MQ=255)
                                                                                             cGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACTTGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGa                                                                                                      <  2:434538/78‑1 (MQ=255)
                                                                                             cGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACTTGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGCACCCt                                         >  1:368882/1‑139 (MQ=255)
                                                                                             cGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACTTGCAAACCTGCGTGTCACGCATAAATCCCTCTTCGa                                                                                                      >  1:434538/1‑78 (MQ=255)
                                                                                               cccAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCTGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGTATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGCACCCTGa                                       >  1:251731/1‑139 (MQ=255)
                                                                                                 cAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACTTGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTa                                                                                                 <  2:252300/79‑1 (MQ=255)
                                                                                                 cAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACTTGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTa                                                                                                 >  1:252300/1‑79 (MQ=255)
                                                                                                  aGAACCACGGTTTTGCGGTGGACGAAGCAACATTACTTGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGCACCCTGAagc                                    <  2:371552/139‑1 (MQ=255)
                                                                                                  aGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCTGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGCACCCTGAagc                                    >  1:198023/1‑139 (MQ=255)
                                                                                                    aaCCACGGTTTTGCGGTGGACGAAGCAACATTACTTGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGa                                                                                                      >  1:585895/1‑71 (MQ=255)
                                                                                                    aaCCACGGTTTTGCGGTGGACGAAGCAACATTACTTGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGa                                                                                                      <  2:585895/71‑1 (MQ=255)
                                                                                                       cACGGTTTTGCGGTGGACGAAGCAACATTACTTGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGCACCCTGAagccagcc                               <  2:247967/139‑1 (MQ=255)
                                                                                                         cGGTTTTGCGGTGGACGAAGCAACATTACTTGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGCACCCTGAAGCCAGCCCt                             >  1:502466/1‑139 (MQ=255)
                                                                                                              ttGCGGTGGACGAAGCAACATTACTTGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGCACCCTGAAGCCAGCCCTGGTcc                        >  1:220172/1‑139 (MQ=255)
                                                                                                                                    aCCTGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGCACCCTGAAGCCAGCCCTGGTCCACACGACGCCGCGCCGTTGTTc  <  1:105657/139‑1 (MQ=255)
                                                                                                                                     ccTGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGCACCCTGAAGCCAGCCCTGGTCCACacgacg                <  1:168093/124‑1 (MQ=255)
                                                                                                                                     ccTGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGCACCCTGAAGCCAGCCCTGGTCCACacgacg                >  2:168093/1‑124 (MQ=255)
                                                                                                                                      |                                                                                                                                        
GGCGAGCGGTGCGAAGACTGTCAAAATGAAATTTGGTCACCACGGCGGCAACCATCCGGTTAAAGATGTGGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCTGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGCACCCTGAAGCCAGCCCTGGTCCACACGACGCCGCGCCGTTGTTC  >  NC_000913/30481‑30751

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