Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,221,347 T→C I381V (ATT→GTT)  arpA ← regulator of acetyl CoA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,221,3470TC94.4% 51.3 / ‑3.9 18I381V (ATT→GTT) arpAregulator of acetyl CoA synthetase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (7/10);  total (7/11)
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.93e-01

CTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAATTTATAGCCAAATATAAACCTGGCAGGTTAATACTATTATAGGCAGTCAACATTTCTGTTCTATGTAGTTCTGGCATTTTCTGGAGTTGAATCATCAGTTG  >  NC_000913/4221269‑4221416
                                                                              |                                                                     
cTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAATTTATAGCCAAATATAAACCTGGCAGGTTAacac                                                                     >  1:861978/1‑81 (MQ=255)
  gACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAATTTATAGCCAAATATAAACCTGGCAGGTTAACACta                                                                   >  1:861979/1‑81 (MQ=255)
         aaTTGAAAATAGTCTCAACGATATCCGCATTACCATAATTTATAGCCAAATATAAACCTGGCAGGTTAACACTATTATAg                                                             <  1:861949/80‑1 (MQ=255)
           ttGAAAATAGTCTCAACGATATCCGCATTACCATAATTTATAGCCAAATATAAACCTGGCAGGTTAACACTATTATAGGCa                                                          <  2:861950/81‑1 (MQ=255)
                      ctcAACGATATCCGCATTACCATAATTTATAGCCAAATATAAACCTGGCAGGTTAACACTATTATAGGCAGTCAACATTTc                                               >  1:861982/1‑81 (MQ=255)
                         aaCGATATCCGCATTACCATAATTTATAGCCAAATATAAACCTGGCAGGTTAACACTATTATAGGCAGTCAACATTTc                                               <  2:861951/78‑1 (MQ=255)
                           cGATATCCGCATTACCATAATTTATAGCCAAATATAAACCTGGCAAGTTAATACTATTATAGGCAGTCAACATTTCTGTTc                                          <  2:861952/81‑1 (MQ=255)
                            gATATCCGCATTACCATAATTTATAGCCAAATATAAACCTGGCAGGTTAACACTATTATAGGCAGTCAACATTTCTGTTCt                                         <  2:861953/81‑1 (MQ=255)
                                tCCGCATTACCATAATTTATAGCCAAATATAAACCTGGCAGGTTAACACTATTATAGGCAGTCAACATTTCTGTTCTATGt                                     <  1:861954/81‑1 (MQ=255)
                                      ttACCATAATTTATAGCCAAATATAAACCTGGCAGGTTAACACTATTATAGGCAGTCAACATTTCTGTTCTATGTAGTTCt                               <  2:861955/81‑1 (MQ=255)
                                              aTTTATAGCCAAATATAAACCTGGCAGGTTAACACTATTATAGGCAGTCAACATTTCTGTTCTATGTAGTTCTGGCAtttt                       <  2:861956/81‑1 (MQ=255)
                                                 tataGCCAAATATAAACCTGGCAGGTTAACACTATTATAGGCAGTCAACATTTCTGTTCTATGTAGTTCTGGCATTTTCTg                    >  2:861984/1‑81 (MQ=255)
                                                     gCCAAATATAAACCTGGCAGGTTAACACTATTATAGGCAGTCAACATTTCTGTTCTATGTAGTTCTGGCATTTTCTGGAGt                <  2:861957/81‑1 (MQ=255)
                                                       cAAATATAAACCTGGCAGGTTAACACTATTATAGGCAGTCAACATTTCTGTTCTATGTAGTTCTGGCATTTTCTGGAGtt               >  2:861990/1‑80 (MQ=255)
                                                       cAAATATAAACCTGGCAGGTTAACACTATTATAGGCAGTCAACATTTCTGTTCTATGTAGTTCTGGCATTTTCTGGAGTTg              <  2:861958/81‑1 (MQ=255)
                                                        aaaTATAAACCTGGCAGGTTAACACTATTATAGGCAGTCAACATTTCTGTTCTATGTAGTTCTGGCATTTTCTGGAGTTg              >  1:861989/1‑80 (MQ=255)
                                                            ataAACCTGGCAGGTTAACACTATTATAGGCAGTCAACATTTCTGTTCTATGTAGTTCTGGCATTTTCTGGAGTTGAatca         >  2:861988/1‑81 (MQ=255)
                                                              aaaCCTGGCAGGTTAACACTATTATAGGCAGTCAACATTTCTGTTCTATGTAGTTCTGGCATTTTCTGGAGTTGAatcat        >  1:861980/1‑80 (MQ=255)
                                                                     gCAGGTTAACACTATTATAGGCAGTCAACATTTCTGTTCTATGTAGTTCTGGCATTTTCTGGAGTTGAATCATCAGTTg  <  2:861959/79‑1 (MQ=255)
                                                                              |                                                                     
CTGACAATGAATTGAAAATAGTCTCAACGATATCCGCATTACCATAATTTATAGCCAAATATAAACCTGGCAGGTTAATACTATTATAGGCAGTCAACATTTCTGTTCTATGTAGTTCTGGCATTTTCTGGAGTTGAATCATCAGTTG  >  NC_000913/4221269‑4221416

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