Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 4,635,590 T→C 40.0% K56E (AAA→GAA)  uvrY ← two‑componenent system BarA/UvrY ‑ regulatory subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154514,635,5900TC40.0% 21.7 / 32.1 30K56E (AAA→GAA) uvrYtwo‑componenent system BarA/UvrY ‑ regulatory subunit
Reads supporting (aligned to +/- strand):  ref base T (10/8);  new base C (5/7);  total (15/15)
Fisher's exact test for biased strand distribution p-value = 7.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATCTGCGGGCTGATATAGCGCTGGCCGGCAAAGGCCAGACGAATGGCCTGGACCATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGCCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATACCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTAATCAAG  >  AE015451/4635369‑4635756
                                                                                                                                                                                                                             |                                                                                                                                                                      
gATCTGCGGGCTGATATAGCGCTGGCCGGCAAAGGCCAGACGAATGGCCTGGACCATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGCCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGca                                                                                                                                                          >  1:462269/1‑236 (MQ=255)
                                     aGACGAATGGCCTGGACCATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGCCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGc                                                                                                                     >  1:394247/1‑236 (MQ=255)
                                              gCCTGGACCATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGCCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTCCACATCCATCAg                                                                                                                                                            <  1:456325/188‑1 (MQ=255)
                                               ccTGGACCATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGCCGGCAGCGCGGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCg                                                                                                           >  1:495730/1‑236 (MQ=255)
                                               ccTGGACCATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGCCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTCCACATCCATCAGcac                                                                                                                                                         >  1:44213/1‑190 (MQ=255)
                                                   gcccATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGCCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTCCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCg                                                                                                                                    <  1:240339/205‑1 (MQ=255)
                                                         ttCGTCAAGGCCTGCACCCTTGGTCAGGTAGCCGGCCGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACATCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTCCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCt                                                                                                 <  1:486257/236‑1 (MQ=255)
                                                          tCGTCAAGGCCTGCACCCTTGGTCAGGTAGCCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGc                                                                                                                                         <  1:242446/195‑1 (MQ=255)
                                                           cGTCAAGGCCTGCACCCTTGGTCAGGTAGCCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCAGGGCATTTCCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGtt                                                                                                                            >  1:260306/1‑207 (MQ=255)
                                                             tCAAGGCCTGCACCCTTGGTCAGGTAGCCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTCCACATCCATCAGCACAACGTCCGGtcg                                                                                                                                             >  2:486253/1‑186 (MQ=255)
                                                                                       gCCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAg                                                                                                                              >  2:197322/1‑177 (MQ=255)
                                                                                           gCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCaa                                                               <  1:197326/236‑1 (MQ=255)
                                                                                                      cTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATcc                                                          <  2:462265/230‑1 (MQ=255)
                                                                                                                cGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATACCGGTTc                                           <  1:273578/235‑1 (MQ=255)
                                                                                                                                           acCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCt                                                                                     <  2:394243/166‑1 (MQ=255)
                                                                                                                                                cACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATACCGGTTCGCACCAGATCGTGacca                          >  1:117639/1‑220 (MQ=255)
                                                                                                                                                cACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTCCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGcc                                                                                                            <  1:259172/138‑1 (MQ=255)
                                                                                                                                                  ccGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGca                                                                                                                                                          <  1:79775/90‑1 (MQ=255)
                                                                                                                                                      cACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTCCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATACCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTAATc     <  1:478069/235‑1 (MQ=255)
                                                                                                                                                        cTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATACCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTAATCAAg  >  1:127369/1‑236 (MQ=255)
                                                                                                                                                          aCCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTa                                                     >  2:177941/1‑183 (MQ=255)
                                                                                                                                                          aCCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTa                                                     <  1:177945/183‑1 (MQ=255)
                                                                                                                                                            cTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGACtctc                                                                                                    >  2:2481/1‑134 (MQ=255)
                                                                                                                                                                   tCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATACCGGTTCGTACCAGATCGTGATCATCGAcc                    >  2:339027/1‑207 (MQ=255)
                                                                                                                                                                   tCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATACCGGTTCGTACCAGATCGTGATCATCGAcc                    <  1:339031/207‑1 (MQ=255)
                                                                                                                                                                        gTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTCCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATACCGGTTCGTACCAGATCGTGATCATCGACCa                   >  2:159225/1‑203 (MQ=255)
                                                                                                                                                                                                    tCCAGGCCGCCGATGCCGGGCATTTTAAcatccat                                                                                                                                                               >  2:51810/1‑35 (MQ=255)
                                                                                                                                                                                                    tCCAGGCCGCCGATGCCGGGCATTTCCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGTAATCTCCCTCACCCTCCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATAcc                                                <  1:320646/146‑1 (MQ=255)
                                                                                                                                                                                                    tCCAGGCCGCCGATGCCGGGCATTTCCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATAcc                                                >  2:320642/1‑146 (MQ=255)
                                                                                                                                                                                                             ccgATGCCGGGCATTTCCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAg                                                                                                                              <  2:44210/59‑1 (MQ=255)
                                                                                                                                                                                                                             |                                                                                                                                                                      
GATCTGCGGGCTGATATAGCGCTGGCCGGCAAAGGCCAGACGAATGGCCTGGACCATTTCGTCAAGGCCTGCACCCTTGGTCAGGTAGCCGGCAGCGCCGGCCTGCAACAAGCGTGTGGGGAACGGGTCTTCCTCGCACACCGTCACCGCCACTACCTTGATGTCAGGGTGGCTGCGCAGCAGTTTGCGGGTGGCCTCCAGGCCGCCGATGCCGGGCATTTTCACATCCATCAGCACAACGTCCGGTTTCAGCTCCCGGGCCAGTTTCAGTGCCGACTCGCCGGAATCTCCCTCACCCACCACCTGCAGGCCGTCGATGTCGGCCAACATCCGTGTAATACCGGTTCGTACCAGATCGTGATCATCGACCACTAAGACCCTAATCAAG  >  AE015451/4635369‑4635756

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