Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,420,328 G→A 70.4% intergenic (‑61/+160) 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,3280GA70.4% 24.7 / 11.6 24intergenic (‑61/+160)hisJ/argThistidine/lysine/arginine/ornithine transporter subunit/lysine/arginine/ornithine transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (2/5);  new base A (9/7);  total (12/12)
Fisher's exact test for biased strand distribution p-value = 3.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGAATCGCAGCAAACGCCGCAGTTGCGCTGGAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCA  >  CP009273/2420193‑2420465
                                                                                                                                       |                                                                                                                                         
cGGAATCGCAGCAAACGCCGCAGTTGCGCTGGAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGAc                                                                                                                              <  2:271188/149‑1 (MQ=255)
           caAACGCCGCAGTTGCGCTGGAGCAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTAACATAAaattaattgaccatcaag                                                                                                                   >  1:106457‑M1/1‑131 (MQ=255)
               cgccccATTTGCGCTGGAGAAGGCCAGAACCAGAGAGAGCGCTAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGATTATTTTCATGTTTGTTATTTTTATCGTGCCATAAAAGTGACCATCAAGGCAAAata                                                                                                               <  1:11445/149‑1 (MQ=255)
                     aggtgctCTGGAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTAACATAAaattaattgaccatcaaggcaaaatatc                                                                                                         <  2:293143‑M1/146‑29 (MQ=255)
                       ttGCGCTGGAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTAACATAAaattaattgaccctctaggcaaaatccccc                                                                                                       >  2:202456‑M1/1‑119 (MQ=255)
                       ttGCGCTGGAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTGTTTTCATAATGTCCTCAACTGACAGTCTTCTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTAACATAAaattaattgaccatcaaggcaatatatcca                                                                                                       >  2:324313‑M1/1‑119 (MQ=255)
                                     gCCAGAACCAGAGGGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGCCACATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGtatta                                                                                         >  1:224494/1‑149 (MQ=255)
                                           aCCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTAACATAAaattaattgaccatcaaggcaaaatatcccgccccgctagacgtattatt                                                                                   >  2:313390‑M1/1‑99 (MQ=255)
                                              agggagagggaGAGCACCAGTTTTTTCATAGTTTTCTCCACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTAACATAAaattaattgaccatcaaggcaaaatatcccgccccgttatacgtattattttt                                                                                <  2:145550‑M1/146‑54 (MQ=255)
                                                    agCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAAGCAAAATATCCCGCCCCGTTAGACGGTTTATTTTTGATTcaaaaa                                                                          >  2:79184/1‑148 (MQ=255)
                                                      cGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTAACATAAaattaattgaccatcaaggcaaaatatcccgccccgttagacgtattatttttcattaaaa                                                                        <  1:242759‑M1/149‑62 (MQ=255)
                                                             aCCAGTTTTTTTATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTAACATAAaattaattgaccatcaaggcaaaatatcccgccccgttagacgtattctttttcattaaaaatgctta                                                                 >  2:138416‑M1/1‑81 (MQ=255)
                                                                    ttttCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTAACATAAaattaattgaccatcaaggcaaaatatcccgccccgttagacgtattatttttcattaaaaatgcttaagtatgc                                                          >  2:260722‑M1/1‑74 (MQ=255)
                                                                     tttAATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTAACATAAaattaattgaccatcacggccaaatatcctgccccgttagacgtattatttttcattaaaaatgcttaagtatgca                                                         <  2:50120‑M1/149‑77 (MQ=255)
                                                                           aaTGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTAACATAAaattaattgaccatcaaggcaaaatatcccgccccgttagacgtattatttttcattaaaaatgcttacttatgcgattcaa                                                   >  2:27677‑M1/1‑67 (MQ=255)
                                                                                  tCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGTCTGCTagtgaattgaccatcaaggcaaaatatcccgccccgttagacgtattatttttcattaaaaatgcttaagtatgcatttcaatgcacca                                            <  1:332023‑M1/149‑90 (MQ=255)
                                                                                      cTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTAACATAAaattaattgaccatcaaggcaaaatatcccgccccgttagacgtattatttttcattaaaaatgcttaagtatgcatttcaatgcaccacgat                                        >  1:204984‑M1/1‑56 (MQ=255)
                                                                                                ttttATAGTGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCt                              <  1:301179/149‑1 (MQ=255)
                                                                                                ttttATAGGGAATATTTACTGGTTTGTTACACTTATCGTGCCATAAAATTGTCCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCt                              <  1:275528/149‑1 (MQ=255)
                                                                                                  ttATAGGGAATATTTACAGGTTTGTTACACTTATCGTAACATAAaattaattgaccatcaagacaaaatatcccgccccgttagacgtattatttttcattaaaaatgcttaagtatgcatttcaatgcaacacgatagtgcaacaaac                            >  2:95701‑M1/1‑44 (MQ=255)
                                                                                                        ggAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCTCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCa                      >  2:96591/1‑149 (MQ=255)
                                                                                                              tttACAGGTTTGTTACACTTATCGTAACATAAaattaattgaccatcaaggcaaaatatcccgccccgttagacgtattatttttcattaaaaatgcttaagtatgcatttcaatgcaccacgatagtgcaacaaacctggcgtgcacc                <  2:106457‑M1/149‑118 (MQ=255)
                                                                                                                    ggTTTGTTACACTTATCGTAACATAAaattaattgaccatcaaggccaaatatcccgccccgttagacgtattatttttcattaaaaatgcttaagtattcatttcaatgcaccacgatagtgcaacaaacctggcgtgcaccatgatg          <  1:97615‑M1/149‑124 (MQ=255)
                                                                                                                            acacTTATCGTAACATAAaattaattgaccatcaaggcaatagaccccgccccgttagacgtattatttttcattaaaaatgcttaagtatgcatttcaatgcaccacgatagtgcaacaaacctggcgtgcaccatgatgatcattcc  <  1:313390‑M1/149‑132 (MQ=255)
                                                                                                                                       |                                                                                                                                         
CGGAATCGCAGCAAACGCCGCAGTTGCGCTGGAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCA  >  CP009273/2420193‑2420465

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