Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,420,353 A→C 22.7% intergenic (‑86/+135) 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,3530AC22.7% 42.8 / 11.6 31intergenic (‑86/+135)hisJ/argThistidine/lysine/arginine/ornithine transporter subunit/lysine/arginine/ornithine transporter subunit
Reads supporting (aligned to +/- strand):  ref base A (15/8);  new base C (4/3);  total (20/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 = 5.50e-01

AACGCCGCAGTTGCGCTGGAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAGGTACAGCTTCCCAGCGACGTATCAGTCA  >  CP009273/2420206‑2420494
                                                                                                                                                   |                                                                                                                                             
aaCGCCGCAGTTGCGCTGGAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCaaaa                                                                                                                                              <  2:577414/149‑1 (MQ=255)
          ttGCGCTGGAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAGAATTGACCATCAAGGCAAAATATCCCGccc                                                                                                                                    >  2:213006/1‑149 (MQ=255)
          ttGCGCTGGAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATAGTGCCATAAAATTGACCATCCAGGCAAAATATc                                                                                                                                          >  1:288468/1‑143 (MQ=255)
               cTGGAGAAGGCCAGAACCAGAGAGAGTGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTa                                                                                                                               <  2:98513/149‑1 (MQ=255)
                         ccAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAACGTTTACAGGTTTGTTAAACTTATCGTGCCATAAAATTGACCACCAAGACAAAATATCCCGCCCCGCTAGACg                                                                                                                           >  1:423280/1‑143 (MQ=255)
                          cAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACCCTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGCTAGACGt                                                                                                                          >  1:55966/1‑143 (MQ=255)
                                cAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTattatt                                                                                                                    >  1:153986/1‑143 (MQ=255)
                                   ggagcgCGATAGCATCAGTTTTTGCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTAttttt                                                                                                                 <  1:248655/142‑1 (MQ=255)
                                   agagagCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTaa                                                                                                           >  2:554036/1‑149 (MQ=255)
                                      gagCGATAGCACCAGGTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGCCCATCAAGACAAAATATCCCGCCCCGTTAAACGTCTTATTTTTCAt                                                                                                              >  1:332224/1‑143 (MQ=255)
                                         cGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTCGAGGTTTGTTAAGCTTATCGTGACATACACCCGAAAATCCAGGAAGGATATACCGCTCCGGCAGACGCATTATTTTTCACTaa                                                                                                           >  1:190477/1‑143 (MQ=255)
                                          gATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAATCATTTTATACGGAATATTTACATGTTTGTTACACTTACCGTGCCATAAAAGTGACCATCAAGGCAAAATCTCCCGCCCCGTTTGACGTATTATTTTTCATTAAATATGCt                                                                                                    <  2:252491/149‑1 (MQ=255)
                                          gATAGCACCAGTTTTTTCATAATGTCCTCAACTAACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCt                                                                                                    >  2:329760/1‑149 (MQ=255)
                                           aTAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCtt                                                                                                   >  2:543919/1‑149 (MQ=255)
                                                gccAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAACATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATTCt                                                                                                    >  1:300534/2‑143 (MQ=255)
                                                  cAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAAAGCTTAATTATg                                                                                            >  2:108907/1‑149 (MQ=255)
                                                                       aaCTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAACATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATTCTTAAGTATGCATTTCAATGcacca                                                                             >  1:70984/1‑143 (MQ=255)
                                                                                tCTTTTTATAGGGAATATTTACCGGTTTGTTACACTTATCGTGCCATAAAATTGACCATAAAGGCAAAATATCCCGCCCCGATAGTCGTATTATTTTTCATTAAAAATGCTTAAGTATGCCGTTCAATGCCCCACGATAGTGc                                                                    >  1:122570/1‑143 (MQ=255)
                                                                                         aGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAACATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATTCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTg                                                     <  2:214068/149‑1 (MQ=255)
                                                                                             aaTATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCAcc                                                 >  2:255677/1‑149 (MQ=255)
                                                                                               tatTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGCCGTg                                                     >  1:133707/1‑143 (MQ=255)
                                                                                                 tttACGGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCatga                                             <  2:635195/149‑1 (MQ=255)
                                                                                                 tttACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCatga                                             <  2:635194/149‑1 (MQ=255)
                                                                                                      aGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATc                                        <  2:346539/149‑1 (MQ=255)
                                                                                                      aGGTTTGTTACACTTATCGTGCAATAAAATTGACCATCAAGGCAAAATATCCCGCCTCGTTAGACGTATTATTCTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATc                                        <  2:252649/149‑1 (MQ=255)
                                                                                                              tACACTTATCGTGCCATAAAATTGACCATCAAGGCAACATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATTCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATc                                <  2:111209/149‑1 (MQ=255)
                                                                                                                 acTTATCGTGCCATAAAATTGACCATCAAGGCAACATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATTCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAgg                             >  2:12180/1‑149 (MQ=255)
                                                                                                                          gCCATAAAATTGACCATCAAGGCAAAATAGCCCGCCCCGTTAGACGTCTTATTTTTCATTAAAAATGCTTAAGTATGCATTTCGATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATTATCATTCCATCAGGTAc                          <  1:422272/143‑1 (MQ=255)
                                                                                                                              tAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATAAGGTACAGCt                      >  1:26965/1‑143 (MQ=255)
                                                                                                                                     gACAATTAAGGAAACATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATTCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAGGTACAGCTTCCCAGCGACGTa         <  2:267590/149‑1 (MQ=255)
                                                                                                                                      aCCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTAGCGTGAAGTGTGGGGATCATTCCATCAGGTACAGCTTCGCAGCGACGGAt        >  2:276761/1‑149 (MQ=255)
                                                                                                                                            aaGGCAACATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATTCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAGGTACAGCTTCCCAGCGACGTATCAGTCa  >  2:445448/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
AACGCCGCAGTTGCGCTGGAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAGGTACAGCTTCCCAGCGACGTATCAGTCA  >  CP009273/2420206‑2420494

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: