Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,060,772 (T)8→9 intergenic (+126/‑78) argP → / → scpA DNA‑binding transcriptional dual regulator ArgP/methylmalonyl‑CoA mutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,060,7641.T96.7% 103.2 / ‑3.2 30intergenic (+118/‑86)argP/scpADNA‑binding transcriptional dual regulator ArgP/methylmalonyl‑CoA mutase
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base T (15/14);  total (15/15)
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.91e-01

CCTCACTTCTTTTTGTATTCCTTGAATCACATCACAAAATAGACAAATCTCAGGCGGCAAAAAACGACGTCTGAATGCA‑‑TTTTTTTTGCTGGCGACAAACCCACGTAAAAAGCTCACCGTAGGCGCAAATACCCTCATTTTGATTGCGTTTTAC  >  NC_000913/3060686‑3060839
                                                                               |                                                                            
ccTCACTTCTTTTTGTATTCCTTGAATCACATCACAAAATAGACAAATCTCAGGCGGCAAAAAACGACGTCTGAATGCA‑‑tt                                                                           >  1:503833/1‑81 (MQ=255)
    aCTTCTTTTTGTATTCCTTGAATCACATCACAAAATAGACAAATCTCAGGCGGCAAAAAACGACGTCTGAATGCA‑‑tttt                                                                         <  1:503800/79‑1 (MQ=255)
           tttGTATTCCTTGAATCACATCACAAAATAGACAAATCTCAGGCGGCAAAAAACGACGTCTGAATGCA‑‑tttttttttg                                                                   <  2:503800/78‑3 (MQ=255)
                 ttCCTTGAATCACATCACAAAATAGACAAATCTCAGGCGGCAAAAAACGACGTCTGAATGCATTTTTTTTTGCTGGCGAc                                                            <  1:503802/80‑1 (MQ=255)
                 ttCCTTGAATCACATCACAAAATAGACAAATCTCAGGCGGCAAAAAACGACGTCTGAATGCATTTTTTTTTGCTGGCGAc                                                            >  2:503830/1‑80 (MQ=255)
                   ccTTGAATCACATCACAAAATAGACAAATCTCAGGCGGCAAAAAACGACGTCTGAATGCATTTTTTTTTGCTGGCGACaa                                                          <  2:503802/80‑1 (MQ=255)
                    cTTGAATCACATCACAAAATAGACAAATCTCAGGCGGCAAAAAACGACGTCTGAATGCATTTTTTTTTGCTGGCGACAAAc                                                        >  1:503845/1‑81 (MQ=255)
                      tGAATCACATCACAAAATAGACAAATCTCAGGCGGCAAAAAACGACGTCTGAATGCATTTTTTTTTGCTGGCGACAAAcc                                                       >  2:503828/1‑80 (MQ=255)
                        aaTCACATCACAAAATAGACAAATCTCAGGCGGCAAAAAACGACGTCTGAATGCATTTTTTTTTGCTGGCGACAAACCCAc                                                    <  2:503803/81‑1 (MQ=255)
                          tcacatcacaAAATAGACAAATCTCAGGCGGCAAAAAACGACGTCTGAATGCATTTTTTTTTGCTGGCGACAAACCCACGt                                                  <  2:503804/81‑1 (MQ=255)
                           cacatcacaAAATAGACAAATCTCAGGCGGCAAAAAACGACGTCTGAATGCATTTTTTTTTTGCTGGCGACAAACCCACg                                                   >  2:503839/1‑80 (MQ=255)
                           cacatcacaAAATAGACAAATCTCAGGCGGCAAAAAACGACGTCTGAATGCATTTTTTTTTGCTGGCGACAAACCCACGt                                                  >  2:503831/1‑80 (MQ=255)
                            acatcacaAAATAGACAAATCTCAGGCGGCAAAAAACGACGTCTGAATGCATTTTTTTTTGCTGGCGACAAACCCACGTaa                                                <  1:503806/81‑1 (MQ=255)
                                 acaAAATAGACAAATCTCAGGCGGCAAAAAACGACGTCTGAATGCATTTTTTTTTGCTGGCGACAAACCCACGTAAAAAg                                            <  2:503806/80‑1 (MQ=255)
                                   aaaaTAGACAAATCTCAGGCGGCAAAAAACGACGTCTGAATGCATTTTTTTTTGCTGGCGACAAACCCACGTAAAAAGCTc                                         <  1:503808/81‑1 (MQ=255)
                                         gACAAATCTCAGGCGGCAAAAAACGACGTCTGAATGCA‑‑TTTTTTTTGCTGGCGACAAACCCACGTAAAAAGCTCACCGTAg                                  <  1:503809/81‑1 (MQ=255)
                                            aaaTCTCAGGCGGCAAAAAACGACGTCTGAATGCATTTTTTTTTGCTGGCGACAAACCCACGTAAAAAGCTCACCGTAGgc                                <  1:503810/81‑1 (MQ=255)
                                              aTCTCAGGCGGCAAAAAACGACGTCTGAATGCATTTTTTTTTGCTGGCGACAAACCCACGTAAAAAGCTCACCGTAGgcgc                              <  2:503810/81‑1 (MQ=255)
                                                 tcAGGCGGCAAAAAACGACGTCTGAATGCATTTTTTTTTGCTGGCGACAAACCCACGTAAAAAGCTCACCGTAGGCGCaaa                           >  1:503843/1‑81 (MQ=255)
                                                  cAGGCGGCAAAAAACGACGTCTGAATGCATTTTTTTTTGCTGGCGACAAACCCACGTAAAAAGCTCACCGTAGGCGCaaa                           >  2:503833/1‑80 (MQ=255)
                                                   aGGCGGCAAAAAACGACGTCTGAATGCATTTTTTTTTGCTGGCGACAAACCCACGTAAAAAGCTCACCGTAGGCGCaa                            <  1:503812/78‑1 (MQ=255)
                                                        gcAAAAAACGACGTCTGAATGCATTTTTTTTTGCTGGCGACAAACCCACGTAAAAAGCTCACCGTAGGCGCAAATACCCTc                    <  1:503813/81‑1 (MQ=255)
                                                               acgacgTCTGAATGCATTTTTTTTTGCTGGCGACAAACCCACGTAAAAAGCTCACCGTAGGCGCAAATACCCTCATTTTGa             >  1:503836/1‑81 (MQ=255)
                                                                cgacgTCTGAATGCATTTTTTTTTGCTGGCGACAAACCCACGTAAAAAGCTCACCGTAGGCGCAAATACCCTCATTTTGAt            >  2:503836/1‑81 (MQ=255)
                                                                 gacgTCTGAATGCATTTTTTTTTGCTGGCGACAAACCCACGTAAAAAGCTCACCGTAGGCGCAAATACCCTCATTTTGAtt           >  2:503841/1‑81 (MQ=255)
                                                                 gacgTCTGAATGCATTTTTTTTTGCTGGCGACAAACCCACGTAAAAAGCTCACCGTAGGCGCAAATACCCTCATTTTGAtt           >  2:503846/1‑81 (MQ=255)
                                                                     tCTGAATGCATTTTTTTTTGCTGGCGACAAACCCACGTAAAAAGCTCACCGTAGGCGCAAATACCCTCATTTTGATTGCGt       <  2:503813/81‑1 (MQ=255)
                                                                     tCTGAATGCATTTTTTTTTGCTGGCGACAAACCCACGTAAAAAGCTCACCGTAGGCGCAAATACCCTCATTTTGATTGCGt       >  2:503838/1‑81 (MQ=255)
                                                                         aaTGCATTTTTTTTTGCTGGCGACAAACCCACGTAAAAAGCTCACCGTAGGCGCAAATACCCTCATTTTGATTGCGtttt    <  2:503815/80‑1 (MQ=255)
                                                                         aaTGCATTTTTTTTTGCTGGCGACAAACCCACGTAAAAAGCTCACCGTAGGCGCAAATACCCTCATTTTGATTGCGttt     <  1:503817/79‑1 (MQ=255)
                                                                         aaTGCATTTTTTTTTGCTGGCGACAAACCCACGTAAAAAGCTCACCGTAGGCGCAAATACCCTCATTTTGATTGCGTTTTa   >  1:503844/1‑81 (MQ=255)
                                                                         aaTGCATTTTTTTTTGCTGGCGACAAACCCACGTAAAAAGCTCACCGTAGGCGCAAATACCCTCATTTTGATTGCGTTTTa   >  2:503848/1‑81 (MQ=255)
                                                                          aTGCATTTTTTTTTGCTGGCGACAAACCCACGTAAAAAGCTCACCGTAGGCGCAAATACCCTCATTTTGATTGCGTTTTAc  >  1:503838/1‑81 (MQ=255)
                                                                               |                                                                            
CCTCACTTCTTTTTGTATTCCTTGAATCACATCACAAAATAGACAAATCTCAGGCGGCAAAAAACGACGTCTGAATGCA‑‑TTTTTTTTGCTGGCGACAAACCCACGTAAAAAGCTCACCGTAGGCGCAAATACCCTCATTTTGATTGCGTTTTAC  >  NC_000913/3060686‑3060839

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