Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,420,353 A→G 68.6% 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,3530AG68.6% 47.2 / 30.6 48intergenic (‑86/+135)hisJ/argThistidine/lysine/arginine/ornithine transporter subunit/lysine/arginine/ornithine transporter subunit
Reads supporting (aligned to +/- strand):  ref base A (8/7);  new base G (17/16);  total (25/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.98e-01

GCCGCAGTTGCGCTGGAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAGGTACAGCTTCCCAGCGACGTATCAGTCA  >  CP009273/2420209‑2420494
                                                                                                                                                |                                                                                                                                             
gccgcAGTTGCGCTGGAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGatat                                                                                                                                           >  1:174499/1‑149 (MQ=255)
gccgcAGTTGCGCTGGAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAatat                                                                                                                                           >  1:345948/1‑149 (MQ=255)
          cgcTGGAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGt                                                                                                                                 <  1:425980/149‑1 (MQ=255)
           gcTGGAGAAGGCAAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGtt                                                                                                                                <  1:143826/149‑1 (MQ=255)
                  aaGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTACAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCACGATATCCCTCCCTGTTAGACGta                                                                                                                         >  2:468708/1‑149 (MQ=255)
                   aGGCCAGAACCAGGGAGAGCGATAGCACGAGCTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGtat                                                                                                                        <  1:314166/149‑1 (MQ=255)
                     gCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGtatta                                                                                                                      >  2:243111/1‑149 (MQ=255)
                       cAGAACCAGAGAGAGCGATAGCACCAGTCTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGGTAGACGTattatt                                                                                                                    >  1:305494/1‑149 (MQ=255)
                        aGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTAttt                                                                                                                   >  1:173745/1‑149 (MQ=255)
                          aaCCAGAGAGATCGATAGCATCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAAGATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTAttttt                                                                                                                 <  1:152784/149‑1 (MQ=255)
                          aaCCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTCttttt                                                                                                                 <  2:174499/149‑1 (MQ=255)
                             cAGAGAGAGCGATAGCACCAGTTTTTTTTTTTTTTTTTCAACTGAGAGTCTTTTTATCGGGTATATTTTCAGGTTTGTTACACTTATCGTGCCATAATCTTGTCCATCGAGGCAAGATATCCCGCCCAGTTATACGTTTTATTTTTCAt                                                                                                              >  2:200619/1‑149 (MQ=255)
                             cAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTGTTGGTTTTCAt                                                                                                              >  1:60669/1‑149 (MQ=255)
                                 gagagCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTaaa                                                                                                          <  1:277026/149‑1 (MQ=255)
                                        aTAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAACTTGACCATCAAGGCAAAATATCCCGCCCCGTTACACGTATTCTTTTTCATTAAAAATGAtt                                                                                                   >  2:36149/1‑149 (MQ=255)
                                            cACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGt                                                                                               >  2:288085/1‑149 (MQ=255)
                                              ccAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTACACGTATTATTTTTCATTAAAAATGCTTAAGTAt                                                                                             >  2:186070/1‑149 (MQ=255)
                                               cAGTTTTTTCATAATGTCCTCAACGTACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTGTCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATg                                                                                            <  2:224316/149‑1 (MQ=255)
                                                 gTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCa                                                                                          <  2:387131/149‑1 (MQ=255)
                                                        cATAATGTCCTCAACTAACAGTCTTTTGATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAAt                                                                                   <  2:14566/149‑1 (MQ=255)
                                                         aTAATGTCCTCAACTGACAGTAATTTTATAGGGAATATCTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATg                                                                                  <  1:385371/149‑1 (MQ=255)
                                                           aaTGTCCTCCACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGca                                                                                <  2:407039/149‑1 (MQ=255)
                                                                       tGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGc                                                                    <  2:10882/149‑1 (MQ=255)
                                                                         ataGTCTTTTTATAGAGAATTTTTACAGCTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGcaa                                                                  <  1:266417/147‑1 (MQ=255)
                                                                           aGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGcaaca                                                                >  2:90196/1‑149 (MQ=255)
                                                                           aGTCTTTTTATAGGGAATATTTACAGCTGTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACAACGATAGTGcaaca                                                                >  2:54288/1‑149 (MQ=255)
                                                                                ttttATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCt                                                           >  2:465640/1‑149 (MQ=255)
                                                                                 tttATAGGGCATATTTACAGGTGTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATTCCGCCCCGTTAGCCGTATTATTTTTCATTAAAGAGCCATAAGTATGCATTTCAATGCAACACGATAATGCAACAAACCTg                                                          >  2:228896/1‑149 (MQ=255)
                                                                                   tataGGGAATATTTACAGGTTTGTTACGCTTATCGTGCCATAAAATTGACCATAAAGGCAAAATCTCTCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGc                                                        <  2:400819/149‑1 (MQ=255)
                                                                                     taGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGt                                                      <  1:288085/149‑1 (MQ=255)
                                                                                       gggAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCTTTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTga                                                    >  2:25674/1‑148 (MQ=255)
                                                                                        ggAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCa                                                   >  1:407224/1‑149 (MQ=255)
                                                                                        ggAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCAAAAAGCTTGTGCAACAAACCTGGCGTgga                                                   >  2:447524/1‑147 (MQ=255)
                                                                                              tttACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCatga                                             >  1:116581/1‑149 (MQ=255)
                                                                                              tttACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCatga                                             >  1:271967/1‑149 (MQ=255)
                                                                                                  cAGGTTTGGTACACTTAACGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATgatgat                                         <  2:147131/149‑1 (MQ=255)
                                                                                                   tggTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATc                                        <  1:150979/148‑1 (MQ=255)
                                                                                                    ggTTTTTTACACTTATCGGGCCATAAAATTGCCCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCa                                       <  2:373725/149‑1 (MQ=255)
                                                                                                     gTTTGGTACACTTCTCGTGCCTTAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCAt                                      <  2:201257/149‑1 (MQ=255)
                                                                                                        tGTTACACTTATCGGGCCATAAAATCGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTcc                                   <  2:148818/149‑1 (MQ=255)
                                                                                                              acTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAgg                             >  1:136896/1‑149 (MQ=255)
                                                                                                                         cATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAGGTACAGCTTccc                  >  2:404565/1‑149 (MQ=255)
                                                                                                                          agaaaaTTGACCATCAAGGCAAGATATCCCCCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAGGTACAGCTTCCCa                 <  2:190982/147‑1 (MQ=255)
                                                                                                                          aaaaaaTTGACCATCAAGGCAAGAAATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAGGTACAGCTTCCCa                 <  2:324345/147‑1 (MQ=255)
                                                                                                                                      aTCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAGGTACAGCTTCCCAGCGACGTATCAg     >  2:312754/1‑149 (MQ=255)
                                                                                                                                      aTCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAGGTACAGCTTCCCAGCGACGTATCAg     >  2:308104/1‑149 (MQ=255)
                                                                                                                                             cAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAGGTACAGCTTCCCAGCGACGTATCAGTCa  >  1:43739/1‑145 (MQ=255)
                                                                                                                                             cAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAGGTACAGCTTCCCAGCGACGTATCAGTCa  <  2:43739/145‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                             
GCCGCAGTTGCGCTGGAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAGGTACAGCTTCCCAGCGACGTATCAGTCA  >  CP009273/2420209‑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: ‑