Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,420,470 C→A 100% intergenic (‑203/+18) hisJ ← / ← argT histidine/lysine/arginine/ornithine transporter subunit/lysine/arginine/ornithine transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,420,4700CA72.1% 14.8 / 4.8 14intergenic (‑203/+18)hisJ/argThistidine/lysine/arginine/ornithine transporter subunit/lysine/arginine/ornithine transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (3/1);  new base A (3/7);  total (6/8)
Fisher's exact test for biased strand distribution p-value = 2.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01
Rejected as polymorphism: Not supported by required number of reads on both strands.

GCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAGGTACAGCTTCCCAGCGACGTATCAGTCACCGTAGACATTAAAGTCGAAATACTTTTTCGCCATCTTGTCGTAGGTGCCGTCCTGACGCAGCTCGCCAAGCGCCTTATTGAAGGCAGCCGTCAG  >  CP009273/2420328‑2420589
                                                                                                                                              |                                                                                                                       
cgccgcagttgcgctggagaaggccagaaccagagagagcgatagcaccagttttttcataatgtcctcaactgacagtctttttatagggaatatttacaggtttgttacacttatcgtgccataaaaATTCCATCAGGTAAAGCTTc                                                                                                                   <  1:105789‑M2/20‑1 (MQ=255)
     caggtgcgctggagtaggccagaaccagagagagtgatagcaccagttttttcataatggcctcaactgacagtctttttatagggaatatttacaggtttgttacacttatcgtgccataaaaATTCCATCAGGTAAAGCTTCCCAGc                                                                                                              <  2:268921‑M2/25‑1 (MQ=255)
       aaTTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAGGTACAGCTTCCCAGCGa                                                                                                            >  2:55984/1‑149 (MQ=255)
                    ggCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAGGTACAGCTTCCCAGCGACGTATCAGTCAcc                                                                                               >  2:192823/1‑149 (MQ=255)
                      gccagaaccagagagagcgatagcaccagttttttcataatgtcctcaactgacagtctttttatagggaatatttacaggtttgttacacttatcgtgccataaaaATTCCATCAGGTAAAGCTTCCCAGCGACGTATCAGTCACCGt                                                                                             <  1:282788‑M2/42‑1 (MQ=255)
                       aaaaTATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAGGTACAGCTTCCCAGCGACGTATCAGTCACCGTa                                                                                            <  2:279649/149‑1 (MQ=255)
                                    gagcgatagcaccagttttttcataatgtcctcaactgacagtctttttatagggaatatttacaggtttgttacacttatcgtgccataaaaATTCCATCAGGTAAAGCTTCCCAGCGACGTATCAGTCACCGTAGACATTAAAGTCg                                                                               <  2:246664‑M2/56‑1 (MQ=255)
                                             caccagttttttcataatgtcctcaactgacagtctttttatagggaatatttacatgtttgttacacttatcgtgccataaaaATTCCATCAGGTAAAGCTTCCCAGCGACGTATCAGTCACCGTAGACATTAAAGTCGAAATACttt                                                                      >  2:175762‑M2/85‑149 (MQ=255)
                                                cggttttttcataatgtgctgaactgacagtctttttatagggaatatttacaggtttgttacacttatcgtgccataaaaATTCCATCAGGTAAAGCTTCCCAGCGACGTATCAGTCACCGTAGACATTAAAGTCGAAATACTTTTTc                                                                   <  1:165288‑M2/68‑1 (MQ=255)
                                                         cataatgtcctcaactgacagtctttttatagggaatatttacaggtttgttacacttatcgtgccataaaaATTCCATCAGGTAAAGCTTCCCAGCGACGTATCAGTCACCGTAGACATTAAAGTCGAAATACTTTTTCGCCATCTTg                                                          >  2:184664‑M2/73‑149 (MQ=255)
                                                                               ctttttatagggaatatttacaggtttgttacacttatcgtgccataaacATTCCATCAGGTAAAGCTTCCCAGCGACGTATCAGTCACCGTAGACATTAAAGTCGAAATACTTTTTCGCCATCTTGTCGTAGGTGCCGTCCTGACGCa                                    <  1:267245‑M2/99‑1 (MQ=255)
                                                                                   aaTGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAGGTACAGCTTCCCAGCGACGTATCAGTCACCGTAGACATTAAAGTCGAAATACTTTTTCGCCATCTTGTCGTAGGTGCCGTCCTGACGCAGCTc                                >  2:272276/1‑149 (MQ=255)
                                                                                                     tgtttgttacacttatcgtgccataaaaATTCCATCAGGTAAAGCTTCCCAGCGACGTATCAGTCACCGTAGACATTAAAGTCGAAATACTTTTTCGCCATCTTGTCGTAGGTGCCGTCCTGACGCAGCTCGCCAAGCGCCTTATTGaa              <  1:175762‑M2/121‑1 (MQ=255)
                                                                                                                 ttatcgtgccataaaaATTCCATCAGGTAAAGCTTCCCAGCGACGTATCAGTCACCGTAGACATTAAAGTCGAAATACTTTTTCGCCATCTTGTCGTAGGTGCCGTCCTGACGCAGCTCGCCAAGCGCCTTATTGAAGGCAGCCgtccg  >  2:124171‑M2/17‑147 (MQ=255)
                                                                                                                                              |                                                                                                                       
GCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAGGTACAGCTTCCCAGCGACGTATCAGTCACCGTAGACATTAAAGTCGAAATACTTTTTCGCCATCTTGTCGTAGGTGCCGTCCTGACGCAGCTCGCCAAGCGCCTTATTGAAGGCAGCCGTCAG  >  CP009273/2420328‑2420589

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: