Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009273 2,420,328 G→A 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,3280GA98.0% 137.8 / NA 51intergenic (‑61/+160)hisJ/argThistidine/lysine/arginine/ornithine transporter subunit/lysine/arginine/ornithine transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (28/22);  minor base C (1/0);  total (29/22)
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.99e-01

TTTGCGGAATCGCAGCAAACGCCGCAGTTGCGCTGGAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCAT‑‑‑‑A‑‑‑‑AAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCAT  >  CP009273/2420189‑2420463
                                                                                                                                            |                                                                                                                                               
tttGTGGAATCGCAGCACACGCCGTAGTTGCGCTGGAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCTTAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTACCATaaaaa                                                                                                                                        <  1:216742‑M1/149‑6 (MQ=255)
tttGCGGAATCGCAGCAAACGCCGCAGTTGCGCTGGAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTACCATaaaca                                                                                                                                        >  1:240605‑M1/1‑144 (MQ=255)
    cGGAATCGCAGCAAACGCCGCAGTTGCGCTGGAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTACCATaaaaaaaat                                                                                                                                    >  1:761081‑M1/1‑140 (MQ=255)
     ggAATCGCAGCAAACGCCGCAGTTGCGCTGGAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTCCCATaaaaaaaatt                                                                                                                                   >  2:68701‑M1/1‑139 (MQ=255)
     ggAATCGCAGCAAACGCCGCAGTTGCGCTGGAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTACCATaaaaaaaatt                                                                                                                                   >  1:95457‑M1/1‑139 (MQ=255)
       acacgcaGCAAACGCCGCAGTTGCGCGGGAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTACCATaaaaaaaattga                                                                                                                                 <  1:350855‑M1/146‑13 (MQ=255)
       aaTCGCAGCAAACGCCGCAGTTGCGACGGCGAAGGCCAGAACCAGAGAGAGCTATAGCATCAGTGTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTACCATaaaaaaaattga                                                                                                                                 <  2:2643‑M1/149‑13 (MQ=255)
         tCGCAGCAACCGCCGCCATTGCGCTGGAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTGTTTATAGGGAATATGTACGGGTTTGTTACACTTGTCGCACCATgaaaaaacttgacc                                                                                                                               >  2:185971‑M1/1‑135 (MQ=255)
            cagcaAACGCCGCAGTTGCGCTGGAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATCATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTACCATaaaaaaaattaaccctc                                                                                                                            >  1:688977‑M1/1‑132 (MQ=255)
             agcaAACGCCGCAGTTGCGCTGGAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTACCATaacaaaaattgaccatca                                                                                                                           >  1:36476‑M1/1‑131 (MQ=255)
               caAACGCCGCAGTTGCCCTGGAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATCATGTACTCAACTCCAACTCTTTCTATATGGAATATTTGACAATACTTTACACTTATCGTACCATAACAA‑‑‑‑AAATTGACCATCAAg                                                                                                                     <  2:917969/149‑1 (MQ=255)
                aaaCGCCGCAGTTGCGCTGGAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTGTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTAACGTACCATaaaaaaaattgaccatcaagg                                                                                                                        >  2:610806‑M1/1‑128 (MQ=255)
                     ccgcAGTTGCGCTGGAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTACCATaaaaaaaattgaccatcaaggcaaaa                                                                                                                   <  2:863721‑M1/149‑27 (MQ=255)
                     ccgcAGTTGCGCTGGACAATGCCAGAACCAGAGAGAGCGATAGCACCAGTTGTTTCACAATGTCCTCAACTGACAGTCTTTCTATAGGGAATATTTACAGGTTTGTTACACTTATCGTACCATaaaaaaaattgaccatcaaggcaaaa                                                                                                                   <  2:250378‑M1/149‑27 (MQ=255)
                      cgcAGTTGCGCTGGAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCCACTAGCAGTCTTTTTATAGGGAAGATTTCCAGGGTTGTTACACTTATCGTACCATaaaaaaaattgaacatccgagcaaaaa                                                                                                                  >  1:451165‑M1/1‑122 (MQ=255)
                               gcTGGAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTACCATaaaaaaaattgaccatcaaggcaaaatatcccgccc                                                                                                         >  1:503299‑M1/1‑113 (MQ=255)
                                      aaGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTATTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTACCATaaaaaaaattgaccatcaaggcaaaatatcccgccccgttaga                                                                                                  <  1:111660‑M1/149‑44 (MQ=255)
                                       aGGACAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTACCATaaaaaaaattgaccatcaaggcaaaatatcccgccccgttagac                                                                                                 <  1:455109‑M1/149‑45 (MQ=255)
                                         gCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTACCATaaaaaaaattgaccatcaaggcaaaaaatcccgccccgttagacgt                                                                                               >  2:368002‑M1/1‑103 (MQ=255)
                                           cAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTACCATaaaaaaaattgaccatcaaggcaaaatatcccgccccgttagacgtat                                                                                             >  1:45520‑M1/1‑101 (MQ=255)
                                             gcaccATAGAGATCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTACCATaaaaaaaattgaccatcaaggcaaaatatcccgccccgttagacgtatta                                                                                           <  2:761081‑M1/147‑51 (MQ=255)
                                              aaCCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTACCATaaaaaaaattgaccatcaaggcaaaatatcccgccccgttagacgtattat                                                                                          >  2:447759‑M1/1‑98 (MQ=255)
                                                ccAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTACCATaaaaaaaatggaccatcaaggcaatatatcccgccccgttagacacatgtttt                                                                                        >  1:537205‑M1/1‑96 (MQ=255)
                                                   gggagagcgATAGAACCAGTTTTCTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTACCATaaaaaaaattgaccatcaaggcaaaatatcccgccccgttagacgtattatttttc                                                                                     <  1:190330‑M1/147‑57 (MQ=255)
                                                   gagagcgCTATAGCGCCTGTTTTTGCAGACTGGCCTCAGCTGACAGTCTTTTTCTAGGGAATATTTACAGGTTTGTTTCACTTATCATACCATaaaaaaaattgaccatcaaggcaaaatatcccgccccgttagacgtattatttttc                                                                                     <  2:319099‑M1/149‑57 (MQ=255)
                                                   gagagagCGATCGCACCAGTTTTTTCATATTGTCCTCAACTGCCAGTCTTTATAACACGGTTAATGACAGGTTTGTTACACTTATCGTACCATaaaaaaaattgaccatcaaggcaaaatatcccgccccgttagacgtattatttttc                                                                                     <  1:938091‑M1/149‑57 (MQ=255)
                                                   gagagagCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTCCAGGTATGTTACACTTATCGTACCATaaaaaaaattgaccatcaagccaaaatatcccgccccattagacgtattatttttc                                                                                     >  1:821554‑M1/1‑93 (MQ=255)
                                                   gagagagCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTACCATaaaaaaaattgaccatcaaggcaaaatatcccgccccgttggacgtattatttttc                                                                                     >  1:214014‑M1/1‑93 (MQ=255)
                                                   gagagagCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTACCATaaaaaaaattgaccatcaaggcaaaatatcccgccccgttagacgtattatttttc                                                                                     >  1:214035‑M1/1‑93 (MQ=255)
                                                    agggggcgataGCACCAGTTTTTTCATAATGGTCGCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTACCATaaaaaaaattgaccatcaaggcaaaatatcccgccccgttagacgtattatttttca                                                                                    <  2:75768‑M1/144‑58 (MQ=255)
                                                    agagagCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTACCATaaaaaaaattgaccatcaaggcaaaatatcccgccccgttagacgtattatttttca                                                                                    <  2:548829‑M1/149‑58 (MQ=255)
                                                        agCGATAGCACCAGTTTTTTCATTTTTTTTTCAACCGAAAGTCTGCCTATAGGAAATATTTACAGGTTTGTTGCACTCACCGTACCAT‑‑‑‑AAAAGAAATTGACCATCAAACCAAAATATCCTCCCCCTTTACACGTATAATTCTTCATCaa                                                                            >  2:164812/1‑149 (MQ=255)
                                                          agagtgcatgagttttttaatATTGTCCTGAACTGACAGTCTTTTTATAGGTAATTTTTATATGTTTGTTACACTTATCGTACCATaaaaaaaattgaccatcaaggcaaaatatcccgccccgttagacgtattatttttcattaaaa                                                                              <  2:955897‑M1/139‑64 (MQ=255)
                                                                cACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGGAACATaaaaaaaattgaccatcaaggcaaaatatcccgccccgttagacgtattatttttcattaaaaatgctt                                                                        <  1:76907‑M1/149‑70 (MQ=255)
                                                                 aCCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGCTTGTAACACGTATCGTATCTGataagaaattgaccatcaaggcaaaagagcccgccccgttagacgtattatttttcattaaaaatgctta                                                                       >  1:942747‑M1/1‑79 (MQ=255)
                                                                        ttttCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTACCATaaaaaaaattgaccatcaaggcaaaatatcccgccccgttaaacgtattatttttcattaaaaatgcttacgtatga                                                                >  2:151193‑M1/1‑72 (MQ=255)
                                                                        tttggaTAAGGTCCGCAACGGACAGTCTTTACATAGGCAATATTTATAGGTCTGTTACACTTATAAAAAAAAAAAAA‑‑‑‑AAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGc                                                            <  2:666203/144‑1 (MQ=255)
                                                                         tttCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTACCATaaaaaaaattgaccatcaaggcaaaatatcccgccccgttagacgtattatttttcattaaaaatgcttaagtatgca                                                               >  1:839986‑M1/1‑71 (MQ=255)
                                                                              tAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTACCATaaaaaaaattgaccatcaaggcaaaatatcccgccccgttagacgtattatttttcattaaaaatgcttaagtatgcatttca                                                          >  1:468821‑M1/1‑66 (MQ=255)
                                                                                ctgTCCTCTACGGACGGTCTTTGTATAGGGAATATTTTCAGGTTTGTTACACTTATCCTACCATaaaaaaaattgaccttcaaggcaaaatatcccgccccggtagacgtattatttttcattaaaaatgcttaagtatgcatttcaat                                                        <  2:480523‑M1/148‑86 (MQ=255)
                                                                                         aCTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTACCATaaaaaaaattgaccatcaaggcaaaatatcccgccccgttagacgtattatttttcattaaaaatgcttaagtatgcatttcaatgcaccacga                                               >  1:33432‑M1/1‑55 (MQ=255)
                                                                                              cAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTACCATaaaaaaaattgaccatcaaggcaaaatatcccgccccgttagacgtattatttttcattaaaaatgcttaagtatgcatttcaatgcaccacgatagtg                                          <  2:942747‑M1/149‑100 (MQ=255)
                                                                                                 tCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTACCATaaaaaaaattgaccatcaaggcaaaatatcccgccccgttagacgtattatttttcattaaaaatgcttaagtatgcatttcaatgcaccacgatagtgcaa                                       >  2:548966‑M1/1‑47 (MQ=255)
                                                                                                  cTTTTTATAGGGAATATTTACAGGGTTGTTACACTTATCGTACCATAAAAA‑‑‑‑AAATTGACCATAAAGGGAAAATATCCCGCCCCGGTGGACGTATTATTTTTCATTAAGAATGCTTAAATATGCATTTCAAAGCACCACGAAAGaacaac                                  >  1:91187/1‑149 (MQ=255)
                                                                                                           gggAATAATTACAGGTTGGTTACACTTATCGTACCATaaaaaaaattgaccatcaaggcaaaatatcccgccccgttagacgtattatttttcattaaaaatgcttaagtaagcatttcaatgcgccacgatagtgcaacaaacctggc                             <  1:617991‑M1/149‑113 (MQ=255)
                                                                                                                           tttGTTACACTTATCGTACCATaaaaaaaattgaccatcaaggcaaaatatcccgccccgttagacgtattatttttcattaaaaatgcttaagtatgcatttcaatgcaccacgatagtgcaacaaacctggcgtgcaccatgatgat              >  1:244604‑M1/1‑22 (MQ=255)
                                                                                                                             tGTTACACTTATCGTACCGTaaaaaaaattgaccatcaaggcaaaatatcccgccccgttagacgtattatttttcattaaaaatgcttaagtatgcatttcaatgcaccacgatcgtgcaacaaccctggcgtgcaccatgatgatca            <  1:610806‑M1/149‑130 (MQ=255)
                                                                                                                             tGTTACACTTATCGTACCATaaaaaaaattgaccatcaaggcaaaatatcccgccccgttagacgtattatttttcattaaaaatgcttaagtatgcatttcaatgcaccacgatagtgcaacaaacctggcgtgcaccatgatgatca            <  2:503299‑M1/149‑130 (MQ=255)
                                                                                                                               ttACACTTATCGTACCATaaaaaaaattgaccatcaaggcaaaatatcccgccccgttagacgtattatttttcattaaaaatgctaaagtatgcatttcaatgcaccacgatagtgcaacaaacctggcgtgcaccatgatgatcatt          >  1:208821‑M1/1‑18 (MQ=255)
                                                                                                                                   acTTATCGTACCATaaaaaaaattgaccatcaaggcaaaatatcccgccccgttagacgtattatttttcattaaaaatgcttaagtatgcatttcaatgcaccacgatagggcaacaaacctggcgtgcaccgtgatgatcattccat      >  1:410893‑M1/1‑14 (MQ=255)
                                                                                                                                       aTCGTACCATaaaaaaaattgaccatcaaggcaaaatatcccgccccgttagacgtattatttttcattaaaaatgcttaagtatgcatttcaatgcaccacgatagtgcaacaaacctggcgtgcaccatgatgatcattccatcagg  >  1:450644‑M1/1‑10 (MQ=255)
                                                                                                                                            |                                                                                                                                               
TTTGCGGAATCGCAGCAAACGCCGCAGTTGCGCTGGAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCAT‑‑‑‑A‑‑‑‑AAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCAT  >  CP009273/2420189‑2420463

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: