Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,831,911 T→C I24M (ATA→ATG astC ← succinylornithine transaminase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,831,9110TC93.0% 132.9 / ‑0.1 43I24M (ATA→ATGastCsuccinylornithine transaminase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base C (14/26);  total (15/28)
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.16e-01

ACTCGCCTGTTCGTTCAGCGCTTCACGCAGTTCCGGATGCGCATGGCCCAGCGCGTTCACCGCAATGCCACCCGCGAAGTCGATATACTCTTTCCCCTGCTGATCCCACAAGCGCGAACCTTCGCCACGTACCGGTATAAAGGGTGCCGGAGCGTAAACAGGTATCATCCATTCATCAAAGTTTTCACGCGTAATTGGCTGAGACATAGCGACCTCTACAGTAAATAATTCGTTATTTATATGTTAATAAT  >  NC_000913/1831776‑1832026
                                                                                                                                       |                                                                                                                   
aCTCGCCTGTTCGTTCAGCGCTTCACGCAGTTCCGGATGCGCATGGCCCAGCGCGTTCACCGCAATGCCACCCGCGAAGTCGATATACTCTTTCCCCTGCTGATCCCACAAGCGCGAACCTTCGCCATGTACCGGcata                                                                                                                  >  1:347876/1‑139 (MQ=255)
         ttcgttcAGCGCTTCACGCAGTTCCGGATGCGCATGGCCCAGCGCGTTCACCGCAATGCCACCCGCGAAGTCGATATACTCTTTCCCCTGCTGATCCCACAAGCGCGAACCTTCGCCACGTACCGGCATAAAGGGTGcc                                                                                                         <  1:204389/139‑1 (MQ=255)
            gttcAGCGCTTCACGCAGTTCCGGATGCGCATGGCCCAGCGCGTTCACCGCAATGCCACCCGCGAAGTCGATATACTCTTTCCCCTGCTGATCCCACAAGCGCGAACCTTCGCCACGTACCGGCATAAAGGGTGCCGGa                                                                                                      <  1:201237/139‑1 (MQ=255)
              tcAGCGCTTCACGCAGTTCCGGATGCGCATGGCCCAGCGCGTTCACCGCAATGCCACCCGCGAAGTCGATATACTCTTTCCCCTGCTGATCCCACAAGCGCGAACCTTCGCCACGTACCGGCATAAAGGGTGCCGGAGc                                                                                                    <  2:564097/139‑1 (MQ=255)
                   gcTTCACGCAGTTCCGGATGCGCATGGCCCAGCGCGTTCACCGCAATGCCACCCGCGAAGTCGATATACTCTTTCCCCTGCTGATCCCACAAGCGCGAACCTTCGCCACGTACCGGCATAAAGGGTGCCGGAGCGTaaa                                                                                               <  2:241606/139‑1 (MQ=255)
                   gcTTCACGCAGTTCCGGATGCGCATGGCCCAGCGCGTTCACCGCAATGCCACCCGCGAAGTCGATATACTCTTTCCCCTGCTGATCCCACAAGCGCGAACCTTCGCCACGTACCGGCATAAAGGGTGCCGGAGCGTaaa                                                                                               <  2:653383/139‑1 (MQ=255)
                    cTTCACGCAGTTCCGGATGCGCATGGCCCAGCGCGTTCACCGCAATGCCACCCGCGAAGTCGATATACTCTTTCCCCTGCTGATCCCACAAGCGCGAACCTTCGCCACGTACCGGCATAAAGGGTGCCGGAGCGTAAAc                                                                                              <  1:1054955/139‑1 (MQ=255)
                                                  gcgcGTTCACCGCAATGCCACCCGCGAAGTCGATATACTCTTTCCCCTGCTGATCCCACAAGCGCGAACCTTCGCCATGTACCGGCATAAAGGGTGCCGGAGCGTAAACAGGTATCATCCATTCATCAAAGTTTTCAcg                                                                <  2:347876/139‑1 (MQ=255)
                                                  gcgcGTTCACCGCAATGCCACCCGCGAAGTCGATATACTCTTTCCCCTGCTGATCCCACAAGCGCGAACCTTCGCCACGTACCGGCATAAAGGGTGCCGGAGCGTAAACAGGTATCATCCATTCATCAAAGTTTTCAcg                                                                >  1:568594/1‑139 (MQ=255)
                                                  gcgcGTTCACCGCAATGCCACCCGCGAAGTCGATATACTCTTTCCCCTGCTGATCCCACAAGCGCGAACCTTCGCCACGTACCGGCATAAAGGGTGCCGGAGCGTAAACAGGTATCATCCATTCATCAAAGTTTTCAcg                                                                >  2:223592/1‑139 (MQ=255)
                                                            cGCAATGCCACCCGCGAAGTCGATATACTCTTTCCCCTGCTGATCCCACAAGCGCGAACCTTCGCCACGTACCGGCATAAAGGGTGCCGGAGCGTAAACAGGt                                                                                          >  1:425689/1‑103 (MQ=255)
                                                            cGCAATGCCACCCGCGAAGTCGATATACTCTTTCCCCTGCTGATCCCACAAGCGCGAACCTTCGCCACGTACCGGCATAAAGGGTGCCGGAGCGTAAACAGGt                                                                                          <  2:425689/103‑1 (MQ=255)
                                                             gCAATGCCACCCGCGAAGTCGATATACTCTTTCCCCTGCTGATCCCACAAGCGCGAACCTTCGCCACGTACCGGCATAAAGGGTGCCGGAGCGTAAACAGGTATCATCCATTCATCAAAGTTTTCACGCGTAATTGGCt                                                     <  2:1069051/139‑1 (MQ=255)
                                                             gCAATGCCACCCGCGAAGTCGATATACTCTTTCCCCTGCTGATCCCACAAGCGCGAACCTTCGCCACGTACCGGCATAAAGGGTGCCGGAGCGTAAACAGGTATCATCCATTCATCAAAGTTTTCACGCGTAATTGGCt                                                     <  1:966501/139‑1 (MQ=255)
                                                                 tGCCACCCGCGAAGTCGATATACTCTTTCCCCTGCTGATCCCACAAGCGCGAACCTTCGCCACGTACCGGCATAAAGGGTGCCGGAGCGTAAACAGGTATCATCCATTCATCAAAGTTTTCACGCGt                                                             >  2:926016/1‑127 (MQ=255)
                                                                 tGCCACCCGCGAAGTCGATATACTCTTTCCCCTGCTGATCCCACAAGCGCGAACCTTCGCCACGTACCGGCATAAAGGGTGCCGGAGCGTAAACAGGTATCATCCATTCATCAAAGTTTTCACGCGt                                                             <  1:926016/127‑1 (MQ=255)
                                                                       ccGCGAAGTCGATATACTCTTTCCCCTGCTGATCCCACAAGCGCGAACCTTCGCCACGTACCGGCATAAAGGGTGCCGGAGCGTAAACAGGTATCATCCATTCATCAAAGTTTTCACGCGTAATTGGCTGAGACATAGc                                           <  1:223592/139‑1 (MQ=255)
                                                                            aaGTCGATATACTCTTTCCCCTGCTGATCCCACAAGCGCGAACCTTCGCCACGTACCGGTATAAAGGGTGCCGGAGCGTAAACAGGTATCATCCATTCATCAAAGTTTTCACGCGTAATTGGCTGAGACATAGCGACct                                      <  1:803507/139‑1 (MQ=255)
                                                                                       ctctTTCCCCTGCTGATCCCACAAGCGCGAACCTTCGCCACGTACCGGCATAAAGGGTGCCGGAGCGTAAACAGGTATCATCCATTCATCAAAGTTTTCACGCGTAATTGGCTGAGACATAGCGACCTCTACAGTAaat                           <  2:568594/139‑1 (MQ=255)
                                                                                          tttCCCCTGCTGATCCCACAAGCGCGAACCTTCGCCACGTACCGGCATAAAGGGTGCCGGAGCGTAAACAGGTATCATCCATTCATCAAAGTTTTCACGCGTAATTGGCTGAGACATAGCGACCTCTACAGTAaataat                        <  1:840585/139‑1 (MQ=255)
                                                                                            tCCCCTGCTGATCCCACAAGCGCGAACCTTCGCCACGTACCGGCATAAAGGGTGCCGGAGCGTAAACAGGTATCATCCATTCATCAAAGTTTTCACGCGTAATTGGCTGAGACATAGc                                           >  2:144309/1‑118 (MQ=255)
                                                                                            tCCCCTGCTGATCCCACAAGCGCGAACCTTCGCCACGTACCGGCATAAAGGGTGCCGGAGCGTAAACAGGTATCATCCATTCATCAAAGTTTTCACGCGTAATTGGCTGAGACATAGc                                           <  1:144309/118‑1 (MQ=255)
                                                                                             ccccTGCTGATCCCACAAGCGCGAACCTTCGCCACGTACCGGCATAAAGGGTGCCGGAGCGTAAACAGGTATCATCCATTCATCAAAGTTTTCACGCGTAATTGGCTGAGACATAGCGACCTCTACAGTAAATAATTCg                     <  1:558565/139‑1 (MQ=255)
                                                                                                 tgctgATCCCACAAGCGCGAACCTTCGCCACGTACCGGCATAAAGGGTGCCGGAGCGTAAACAGGTATCATCCATTCATCAAAGTTTTCACGCGTAATTGGCTGAGACATAGCGACCTCTACAGTAAATAATTCGttat                 >  2:799925/1‑139 (MQ=255)
                                                                                                  gctgATCCCACAAGCGCGAACCTTCGCCACGTACCGGCATAAAGGGTGCCGGAGCGTAAACAGGt                                                                                          >  1:113344/1‑65 (MQ=255)
                                                                                                  gctgATCCCACAAGCGCGAACCTTCGCCACGTACCGGCATAAAGGGTGCCGGAGCGTAAACAGGt                                                                                          <  2:113344/65‑1 (MQ=255)
                                                                                                      aTCCCACAAGCGCGAACCTTCGCCACGTACCGGCATAAAGGGTGCCGGAGCGTAAACAGGTATCATCCATTCATCAAAGTTTTCACGCGTAATTGGCTGAGACATAGCGAc                                        <  2:94162/111‑1 (MQ=255)
                                                                                                      aTCCCACAAGCGCGAACCTTCGCCACGTACCGGCATAAAGGGTGCCGGAGCGTAAACAGGTATCATCCATTCATCAAAGTTTTCACGCGTAATTGGCTGAGACATAGCGAc                                        >  1:94162/1‑111 (MQ=255)
                                                                                                      aTCCCACAAGCGCGAACCTTCGCCACGTACCGGCATAAAGGGTGCCGGAGCGTAAACAGGTATCATCCATTCATCAAAGTTTTCACGCGTAATTGGCTGAGACATAGCGACCTCTACAGTAAATAATTCg                     >  1:14248/1‑130 (MQ=255)
                                                                                                      aTCCCACAAGCGCGAACCTTCGCCACGTACCGGCATAAAGGGTGCCGGAGCGTAAACAGGTATCATCCATTCATCAAAGTTTTCACGCGTAATTGGCTGAGACATAGCGACCTCTACAGTAAATAATTCg                     <  2:14248/130‑1 (MQ=255)
                                                                                                      aTCCCACAAGCCCGAACCTTCGCCACGTACCGGCATAAAGGGTGCCAGAGCGTAAACAGGTATCATCCTTTCATCAAAGTTTTCACGCGTAATTGGCTGAGACATAGCGACCTCTACAGTAAATAATTCGTTATTtata            <  2:402681/139‑1 (MQ=255)
                                                                                                        cccACAAGCGCGAACCTTCGCCACGTACCGGCATAAAGGGTGCCGGAGCGTAAACAGGTATCATCCATTCATCAAAGTTTTCACGCGTAATTGGCTGAGACATAGCGAc                                        >  2:397108/1‑109 (MQ=255)
                                                                                                        cccACAAGCGCGAACCTTCGCCACGTACCGGCATAAAGGGTGCCGGAGCGTAAACAGGTATCATCCATTCATCAAAGTTTTCACGCGTAATTGGCTGAGACATAGCGAc                                        <  1:397108/109‑1 (MQ=255)
                                                                                                          cacaAGCGCGAACCTTCGCCACGTACCGGCATAAAGGGTGCCGGAGCGTAAACAGGTATCATCCATTCATCAAAGTTTTCACGCGTAATTGGCTGAGACATAGCGAc                                        <  2:192322/107‑1 (MQ=255)
                                                                                                          cacaAGCGCGAACCTTCGCCACGTACCGGCATAAAGGGTGCCGGAGCGTAAACAGGTATCATCCATTCATCAAAGTTTTCACGCGTAATTGGCTGAGACATAGCGAc                                        >  1:192322/1‑107 (MQ=255)
                                                                                                             aaGCGCGAACCTTCGCTACGTACCGGCATAAAGGGTGCCGGAGCGTAAACAGGTATCATCCATTCATCAAAGTTTTCACGCGTAATTGGCTGAGACATAGCGACCTCTACAGTAAATAATTCGTTATTTATATGTtaat     >  1:914487/1‑139 (MQ=255)
                                                                                                             aaGCGCGAACCTTCGCCACGTACCGGCATAAAGGGTGCCGGAGCGTAAACAGGTATCATCCATTCATCAAAGTTTTCACGCGTAATTGGCTGAGACATAGCGACCTCTACAGTAAATAATTCGTTATTTATATGTtaat     <  1:799925/139‑1 (MQ=255)
                                                                                                               tcgcgAACCTTCGCCACGTACCGGCATAAAGGGTGCCGGAGCGTAAACAGGTATCATCCATTCATCAAAGTTTTCACGCGTAATTGGCTGAGACATAGCGACCTCTACAGTAAATAATTCGTTATTTATATGTtaataa   <  2:839853/138‑1 (MQ=255)
                                                                                                               gcgcGAACCTTCGCCACGTACCGGCATAAAGGGTGCCGGAGCGTAAACAGGTATCATCCATTCATCAAAGTTTTCACGCGTAATTGGCTGAGACATAGCGACCTCTACAGTAAATAATTCGTTATTTATATGTtaataa   <  1:469380/139‑1 (MQ=255)
                                                                                                                cgcgAACCTTCGCCACGTACCGGTATAAAGGGTGCCGGAGCGTAAACAGGTATCATCCATTCATCAAAGTTTTCACGCGTAATTGGCTGAGACATa                                             <  1:937931/96‑1 (MQ=255)
                                                                                                                cgcgAACCTTCGCCACGTACCGGTATAAAGGGTGCCGGAGCGTAAACAGGTATCATCCATTCATCAAAGTTTTCACGCGTAATTGGCTGAGACATa                                             >  2:937931/1‑96 (MQ=255)
                                                                                                                cgcgAACCTTCGCCACGTACCGGCATAAAGGGTGCCGGAGCGTAAACAGGTATCATCCATTCATCAAAGTTTTCACGCGTAATTGGCTGAGACATAGCGACCTCTACAGTAAATAATTCGTTATTTATATGTtaataat  >  2:926683/1‑139 (MQ=255)
                                                                                                                         tCGCCACGTACCGGCATAAAGGGTGCCGGAGCGTAAACAGGTATCATCCATTCATCAAAGTTTTCACGCGTAATTGGCTGAGACATAGCGAc                                        <  2:351761/92‑1 (MQ=255)
                                                                                                                         tCGCCACGTACCGGCATAAAGGGTGCCGGAGCGTAAACAGGTATCATCCATTCATCAAAGTTTTCACGCGTAATTGGCTGAGACATAGCGAc                                        >  1:351761/1‑92 (MQ=255)
                                                                                                                                       |                                                                                                                   
ACTCGCCTGTTCGTTCAGCGCTTCACGCAGTTCCGGATGCGCATGGCCCAGCGCGTTCACCGCAATGCCACCCGCGAAGTCGATATACTCTTTCCCCTGCTGATCCCACAAGCGCGAACCTTCGCCACGTACCGGTATAAAGGGTGCCGGAGCGTAAACAGGTATCATCCATTCATCAAAGTTTTCACGCGTAATTGGCTGAGACATAGCGACCTCTACAGTAAATAATTCGTTATTTATATGTTAATAAT  >  NC_000913/1831776‑1832026

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