Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,318,699 T→C G21G (GGT→GGC argG → argininosuccinate synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,318,6990TC96.0% 79.8 / ‑3.9 25G21G (GGT→GGCargGargininosuccinate synthetase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (10/14);  total (10/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.70e-01

TGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGCGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTGCCATGGAATACGGCG  >  NC_000913/3318570‑3318829
                                                                                                                                 |                                                                                                                                  
tGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGCGGCCTGGACACCa                                                                                                                          >  1:153676/1‑140 (MQ=255)
        aaTTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGCGGCCTGGACACCAGTGCCGCa                                                                                                                  <  2:260377/140‑1 (MQ=255)
                                       atcaGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGCGGCCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTc                                                                                   <  2:144186/140‑1 (MQ=255)
                                         caGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGCGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCt                                                                                 <  1:381860/140‑1 (MQ=255)
                                         caGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGCGGCCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCt                                                                                 <  2:75343/140‑1 (MQ=255)
                                              aattaaGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGCGGCCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGc                                                                            >  1:203610/1‑140 (MQ=255)
                                              aattaaGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGCGGCCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGc                                                                            >  1:256881/1‑140 (MQ=255)
                                              aattaaGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGCGGCCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGc                                                                            >  1:364016/1‑140 (MQ=255)
                                              aattaaGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGCGGCCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGc                                                                            <  2:324019/140‑1 (MQ=255)
                                                 taaGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGCGGCCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCata                                                                         <  1:97313/140‑1 (MQ=255)
                                                   aGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGCGGCCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCAtata                                                                       >  2:333900/1‑140 (MQ=255)
                                                     cAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGCGGCCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACt                                                                     >  2:193448/1‑140 (MQ=255)
                                                           gtTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGCGGCCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAAc                                                               <  1:333900/140‑1 (MQ=255)
                                                                 ttATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGCGGCCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGgc                                                         <  1:278607/140‑1 (MQ=255)
                                                                                      aTCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGCGGCCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTat                                    <  1:390660/140‑1 (MQ=255)
                                                                                          cccGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGCGGCCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTatgatg                                <  2:153676/140‑1 (MQ=255)
                                                                                               tAGGTCAACGTATTGGTATCGCTTTTTCTGGCGGCCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATc                           <  2:203610/140‑1 (MQ=255)
                                                                                                     aaCGTATTGGTATCGCTTTTTCTGGCGGCCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCtcgt                     >  1:86943/1‑140 (MQ=255)
                                                                                                      aCGTATTGGTATCGCTTTTTCTGGCGGCCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCtcgtc                    <  2:86943/140‑1 (MQ=255)
                                                                                                      aCGTATTGGTATCGCTTTTTCTGGCGGCCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCtcgtc                    <  2:256881/140‑1 (MQ=255)
                                                                                                       cGTATTGGTATCGCTTTTTCTGGCGGCCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCtcgtcg                   >  1:36266/1‑140 (MQ=255)
                                                                                                         tATTGGTATCGCTTTTTCTGGCGGCCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTg                 <  2:36266/140‑1 (MQ=255)
                                                                                                            tGGTATCGCTTTTTCTGGCGGCCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTGCCa              >  1:203462/1‑140 (MQ=255)
                                                                                                                   gCTTTTTCTGGCGGCCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTGCCATGGAATa       >  1:171270/1‑140 (MQ=255)
                                                                                                                        ttCTGGCGGCCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTGCCATGGAATACGgcg  <  2:171270/140‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                  
TGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGCGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTGCCATGGAATACGGCG  >  NC_000913/3318570‑3318829

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