Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,420,353 A→G 64.4% 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,3530AG64.4% 22.8 / 23.4 31intergenic (‑86/+135)hisJ/argThistidine/lysine/arginine/ornithine transporter subunit/lysine/arginine/ornithine transporter subunit
Reads supporting (aligned to +/- strand):  ref base A (6/5);  new base G (12/8);  total (18/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

CGCCGCAGTTGCGCTGGAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAGGTACAGCTTCCCAGCGACGTATCAGTCAC  >  CP009273/2420208‑2420495
                                                                                                                                                 |                                                                                                                                              
cgccgcAGTGGCGCTGGAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCAGAAAATTGACCATCAAGGCAAGata                                                                                                                                             <  1:273172/149‑1 (MQ=255)
        ttGCGCTGGAGAGGGCCAGAACCACAGAGCGCGATAGCCCTAGTCTTTTCGGAATGTCCTCAACTGACAGTCTTTTTATAGGGAATCTTCACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATc                                                                                                                                           <  2:231723/143‑1 (MQ=255)
           cgcTGGAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTTTTTTTTTTTTCAACTGAAAGTCTCTTTCTAGGGTATATTTACAGGTTGGGTATACTTAGCGTGCCATTTCCTTTACCCTCAAGGCACGatattcgc                                                                                                                                        >  2:210557/1‑139 (MQ=255)
                    aGGCCAGAACCAGAGAGAGCGATAGCACCAATTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTAAACTTATCGTGCCATAAAATTGCCCATCAAGGCAAGAAATCCCTccccagtag                                                                                                                               >  2:226731/1‑138 (MQ=255)
                        cAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTGGACGt                                                                                                                           >  2:282030/1‑143 (MQ=255)
                                    gagCAATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTATAATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAAACGTATTATTTTTCAt                                                                                                               >  2:35474/1‑143 (MQ=255)
                                                   ttttttCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCAt                                                                                          >  1:172031/1‑149 (MQ=255)
                                                              tGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGcacc                                                                               >  1:296856/1‑149 (MQ=255)
                                                                 ccTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACTCTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACg                                                                            <  1:164853/149‑1 (MQ=255)
                                                                 ccTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACg                                                                            <  1:133854/149‑1 (MQ=255)
                                                                  cTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGa                                                                           >  1:96693/1‑149 (MQ=255)
                                                                    cAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGcacc                                                                               <  2:239948/143‑1 (MQ=255)
                                                                      aCTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGt                                                                       <  1:218760/149‑1 (MQ=255)
                                                                          aCAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGcaa                                                                   <  1:313556/149‑1 (MQ=255)
                                                                             gTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGcaacaa                                                                >  1:131171/1‑149 (MQ=255)
                                                                                tttttATAGGGAATTTTTACAGTTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGcaa                                                                   <  2:77614/143‑1 (MQ=255)
                                                                                tttttATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAAcc                                                             >  1:204263/1‑149 (MQ=255)
                                                                                    tataGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCCCGATAGTGCAACAAACCTGGc                                                         >  1:163482/1‑149 (MQ=255)
                                                                                    tataGGAAATATTTACAGTTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTGCCAATTTATGCTTAAGTATGTATTTCAATGCACCACGATAGTGCAACAAACCTGGc                                                         <  1:254174/149‑1 (MQ=255)
                                                                                       aGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCt                                                            <  2:151422/143‑1 (MQ=255)
                                                                                         ggaatcaatACAGTTCTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTgg                                                          >  2:321433/9‑143 (MQ=255)
                                                                                                ttACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCatgat                                             >  1:202147/1‑149 (MQ=255)
                                                                                                       tttGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCAtt                                      >  1:241685/1‑149 (MQ=255)
                                                                                                              cacTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTc                                     <  2:156933/143‑1 (MQ=255)
                                                                                                              cacTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAg                               <  1:305204/149‑1 (MQ=255)
                                                                                                                cTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTACAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCCGGt                             >  1:17029/1‑149 (MQ=255)
                                                                                                                   aTCGCGCCATAAATTTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCa                                <  2:316053/143‑1 (MQ=255)
                                                                                                                    tCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAg                               >  2:309858/1‑143 (MQ=255)
                                                                                                                         ccATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAGGTACa                          >  2:253888/1‑143 (MQ=255)
                                                                                                                                   gACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAGGTACAGCTTCCCAGc                >  2:222719/1‑143 (MQ=255)
                                                                                                                                       aTCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAGGTACAGCTTCCCAGCGACg            >  2:178307/1‑143 (MQ=255)
                                                                                                                                                 aaTATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAGGTACAGCTTCCCAGCGACGTATCAGTCAc  <  2:116274/143‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
CGCCGCAGTTGCGCTGGAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAGGTACAGCTTCCCAGCGACGTATCAGTCAC  >  CP009273/2420208‑2420495

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