Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 5,408,980 T→C 8.0% K792R (AAG→AGG)  fhuA_3 ← Ferrichrome outer membrane transporter/phagereceptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1435,408,9800TC 8.0% 156.1 / 4.9 50K792R (AAG→AGG) fhuA_3Ferrichrome outer membrane transporter/phagereceptor
Reads supporting (aligned to +/- strand):  ref base T (26/20);  new base C (2/2);  total (28/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 = 2.49e-01

TCAGTCTCTTTGAAGATACGGGGGGTTCATTGCAAGTTTCAGGATCTTGCGCACCTGTGGGCACGGGTTTACCCGCGAACAGGGGCGCAACCCGTGCCATCCCCTGGCTGCCTTACCAGCTATAACTCACCTTCGCCGTCACCTCCCGCCCCGGGTTGTAGTAATTCGCCGTCCCCGACCCCACCACATGCTGCTCGTCAAACAGGTTGCTCACATTCAGCGCCAGGTCGGTCCCTTTGGCGATCTTGTAGTTCAGCGCCGCATC  >  ACB143/5408849‑5409113
                                                                                                                                   |                                                                                                                                     
tcagtcTCTTTGAAGATACGGGGGGTTCATTGCAAGTTTCAGGATCTTGCGCACCTGTGGGCACGGGTTTACCCGCGAACAGGGGCGCAACCCGTGCCATCCCCTGGCTGCCTTACCAGCTATAACTCACCTTCGc                                                                                                                                   >  7:126368/1‑136 (MQ=255)
     ctctTTGAAGATACGGGGGGTTCATTGCAAGTTTCAGGATCTTGCGCACCTGTGGGCACGGGTTTACCCGCGAACAGGGGCGCAACCCGTGCCATCCCCTGGCTGCCTTACCAGCTATAACTCACCTTCGCCGTCa                                                                                                                              >  1:209514/1‑136 (MQ=255)
     ctctTTCAAGATACGGTGGGTTCATTGCAAGTTTCAGGATCAAGCGCACCTGTGCGCACCCGATTACCCGCGAACAGCGCCCCAACACGTCCCAACCCCTGGCTGCCTTACCAGCAAAACCTCCCATTCGCCGTcc                                                                                                                              >  4:215604/1‑135 (MQ=255)
            aaGATACGGGGGGTTCATTGCAAGTTTCAGGATCTTGCGCACCTGTGGGCACGGGTTTACCCGCGAACAGGGGCGCAACCCGTGCCATCCCCTGGCTGCCTTACCAGCTATAACTCACCTTCGCCGTCACCTCCCg                                                                                                                       <  8:76805/136‑1 (MQ=255)
            aaGATACGGGGGGTTCATTGCAAGTTTCAGGATCTTGCGCACCTGTGGGCACGGGTTTACCCGCGAACAGGGGCGCAACCCGTGCCATCCCCTGGCTGCCTTACCAGCTATAACTCACCTTCGCCGTCACCTCCCg                                                                                                                       <  5:113777/136‑1 (MQ=255)
              gATACGGGGGGTTCATTGCAAGTTTCAGGATCTTGCGCACCTGTGGGCACGGGTTTACCCGCGAACAGGGGCGCAACCCGTGCCATCCCCTGGCTGCCTTACCAGCTATAACTCACCTTCGCCGTCACCTCCCGcc                                                                                                                     <  7:231110/136‑1 (MQ=255)
                tACGGGGGGTTCATTGCAAGTTTCAGGATCTTGCGCACCTGTGGGCACGGGTTTACCCGCGAACAGGGGCGCAACCCGTGCCATCCCCTGGCTGCCTTACCAGCTATAACTCACCTTCGCCGTCACCTCCCGcccc                                                                                                                   <  3:61350/136‑1 (MQ=255)
                  cGGGGGGTTCATTGCAAGTTTCAGGATCTTGCGCACGTGTGGGCACGGGTTTACCCGCGAACCGGGGCGAAACCCGCGCCATCCCCTGGCTGCCTTACCAGCTATAACTCACCTTCGCCGTCACCTCCCGCCCCg                                                                                                                  <  4:65130/135‑1 (MQ=255)
                            aTTGCAAGTTTCAGGATCTTGCGCACCTGTGGGCACGGGTTTACCCGCGAACAGGGGCGCAACCCGTGCCATCCCCTGGCTGCCTTACCAGCTATAACTCACCTTCGCCGTCACCTCCCGCCCCGGGTTGTAGTaa                                                                                                       >  5:325446/1‑136 (MQ=255)
                            aTTGCAAGTTTCAGGATCTTGCGCACCTGTGGGCACGGGTTTACCCGCGAACAGGGGCGCAACCCGTGCCATCCACTGGCTGCCTTACCAGCTATAACTCACCTTCGCCGTCACCTACCCCCCCGGGTGGTAGTaa                                                                                                       >  7:292407/1‑136 (MQ=255)
                             ttGCAAGTTTCAGGATCTTGCGCACCTGTGGGCACGGGTTTACCCGCGAACAGGGGCGCAACCCGTGCCATCCCCTGGCTGCCTTACCAGCTATAACTCACCTTCGCCGTCACCTCCCGCCCCGGGTTGTAGTAAt                                                                                                      >  2:351297/1‑136 (MQ=255)
                               gCAAGTTTCAGGATCTTGCGCACCTGTGGGCACGGGTTTACCCGCGAACAGGGGCGCAACCCGTGCCATCCCCTGGCTGCCTTACCAGCTATAACTCACCTTCGCCGTCACCTCCCGCCCCGGGTTGTAGTAAtt                                                                                                     <  4:337039/135‑1 (MQ=255)
                                  aGTTTCAGGATCTTGCGCACCTGTGGGCACGGGTTTACCCGCGAACAGGGGCGCAACCCGTGCCATCCCCTGGCTGCCTTACCAGCTATAACTCACCTTCGCCGTCACCTCCCGCCCCGGGTTGTAGTAATTCGcc                                                                                                 <  6:165586/136‑1 (MQ=255)
                                      tCAGGATCTTGCGCACCTGTGTGCACGGGTTTACCCGCGAACAGGGGCGCAACCCGTGCCATCCCCTGGCTGCCTTACCAGCTATAACTCACCTTCGCCGTCACCTCCCGCCCCGGGTTGTAGTAATTCGCCGTc                                                                                              <  4:391276/135‑1 (MQ=255)
                                      tCAGGATCTTGCGCACCTGTGGGCACGGGTTTACCCGCGAACAGGGGCGCAACCCGTGCCATCCCCTGGCTGCCTTACCAGCTATAACTCACCTTCGCCGTCACCTCCCGCCCCGGGTTGTAGTAATTCGCCGTc                                                                                              <  4:375978/135‑1 (MQ=255)
                                                gcgcACCTGTGGGCACGGGTTTACCCGCGAACAGGGGCGCAACCCGTGCCATCCCCTGGCTGCCTTACCAGCTATAACTCACCTTCGCCGTCACCTCCCGCCCCGGGTTGTAGTAATTCGCCGTCCCCGAc                                                                                        <  5:341443/131‑1 (MQ=255)
                                                 cgcACCTGTGGGCACGGGTTTACCCGCGAACAGGGGCGCAACCCGTGCCATCCCCTGGCTGCCTTACCAGCTATAACTCACCTTCGCCGTCACCTCCCGCCCCGGGTTGTAGTAATTCGCCGTCCCCGACCccac                                                                                   >  6:375173/1‑135 (MQ=255)
                                                         tgGGCACGGGTTTACCCGCGAACAGGGGCGCAACCCGTGCCATCCCCTGGCTGCCTTACCAGCTATAACTCACCTTCGCCGTCACCTCCCGCCCCGGGTTGTAGTAATTCGCCGTCCCCGCCCCCACCACAtgctg                                                                          >  3:271644/1‑136 (MQ=255)
                                                            gCACGGGTTTACCCGCGAACAGGGGCGCAACCCGTGCCATCCCCTGGCTGCCTTACCAGCTATAACTCACCTTCGCCGTCACCTccc                                                                                                                        >  2:337267/1‑87 (MQ=255)
                                                            gCACGGGTTTACCCGCGAACAGGGGCGCAACCCGTGCCATCCCCTGGCTGCCTTACCAGCTATAACTCACCTTCGCCGTCACCTccc                                                                                                                        <  1:337267/87‑1 (MQ=255)
                                                              aCGGGTTTACCCGCGAACAGGGGCGCAACCCGTGCCATCCCCTGGCTGCCTTACCAGCTATAACTCACCTTCGCCGTCACCTCCCGCCCCGGGTTGTAGTAATTCGCCGTCCCCGCCCCCACCACATGCTGCTAGt                                                                     >  7:158188/1‑136 (MQ=255)
                                                                gggTTTACCCGCGAACAGGGGCGCAACCCGTGCCATCCCCTGGCTGCCTTACCAGCTATAACTCACCTTCGCCGTCACCTCCCGCCCCGTGTTGTAGTAATTCGCCGTCCCCGACCCCACCACATGCTGCTCGTCa                                                                   <  3:338478/136‑1 (MQ=255)
                                                                   tttACCCGCGAACAGGGGCGCAACCCGTGCCATCCCCTGGCTGCCTTACCAGCTATAACTCACCTTCGCCGTCACCTCCCGCCCCGGGTTGTAGTAAATCGCCGTCCCCGACCCCACCCCATGCTGCTCGTCaaa                                                                 >  7:52396/1‑135 (MQ=255)
                                                                        caccgAACAGGGGCGCAACACTTTAAATCAACTGGCTGCCTTACCAGCTATAACTCACCTTCGCCGTCACCTCCCGCCCCGGGTTGTAGTAATTCGCCGTCCCCGACCCCACCACATGCTGCTCGTCAAACAGGtt                                                           >  1:247733/4‑136 (MQ=255)
                                                                         cgcgAACAGGGGCGCAACCCGTGCCATCCCCTGGCTGCCTTACCAGCTATAACTCACCTTCGCCGTCACCTCCCGCCCCGGGTTGTAGTAATTCGCCGTCCCCGACCCCACCACATGCTGCTCGTCAAACAGGTTg                                                          <  4:199939/136‑1 (MQ=255)
                                                                              aCAGGGGCGCAACCCGTGCCATCCCCTGGCTGCCTTACCAGCTATAACTCACCTTCGCCGTCACCTCCCGCCCCGGGTTGTAGTAATTCGCCGTCCCCGACCCCACCACATGCTGCTCGTCAAACAGGTTGCTcac                                                     <  4:292787/136‑1 (MQ=255)
                                                                              aCAGGGGCGCAACCCGTGCCATCCCCTGGCTGCCTTACCAGCTATAACTCACCTTCGCCGTCACCTCCCGCCCCGGGTTGTAGTAATTCGCCGTCCCCGACCCCACCACATGCTGCTCGTCAAACAGGTTGCTca                                                      >  3:69752/1‑135 (MQ=255)
                                                                               cAGGGGCGCAACCCGTGCCATCCCCTGGCTGCCTTACCAGCTATAACTCACCTTCGCCGTCACCTCCCGCCCCGGGTTGTAGTAATTCGCCGTCCCCGACCCCACCACATGCTGCTCGTCAAACAGGTTGCTccca                                                    >  4:78046/1‑136 (MQ=255)
                                                                                    gcgcAACCCGTGCCATCCCCTGGCTGCCTTACCAGCTATAACTCACCTTCGCCGTCACCTCCCGCCCCGGGTTGTAGTAATTCGCCGTCCCCGACCCCACCACATGCTGCTCGTCAAACAGGTTGCTCACATTCAg                                               <  1:108260/136‑1 (MQ=255)
                                                                                     cgcAACCCGTGCCATCCCCTGGCTGCCTTACCAGCTATAACTCACCTTCGCCGTCACCTCCCGCCCCGGGTTGTAGTAATTACCAGTCACAGACCCCACCACATTCGGCACGTCAAACAGTTTGCTCAAATTCatc                                              >  6:245240/1‑134 (MQ=255)
                                                                                               gCCATCCCCTGGCTGCCTTACCAGCTATAACTCACCTTCGCCGTCACCTCCCGCCCCGGGTTGTAGTAATTCGCCGTCCCCGCCCCCACCACATGCTGCGCGTCACTCAGGTTGCGCACATTAAGCGCCAggtcgg                                    >  3:245107/1‑136 (MQ=255)
                                                                                                ccATCCCCTGGCTGCCTTACCAGCTATAACTCACCTTCGCCGTCACCTc                                                                                                                          <  3:278339/49‑1 (MQ=255)
                                                                                                ccATCCCCTGGCTGCCTTACCAGCTATAACTCACCTTCGCCGTCACCTc                                                                                                                          >  4:278339/1‑49 (MQ=255)
                                                                                                 cATCCCCTGGCTGCCTTACCAGCTATAACTCACCTTCGCCGTCACCTCCCGCCCCGGGTTGTAGTAATTCGCCGTCCCCGACCCCCCCACATGCTGCTCGGCAAACAGGTTGCTCACATTCAGCGCCAggtcggtc                                  >  1:395325/1‑136 (MQ=255)
                                                                                                  aTCCCCTGGCTGCCTTACCAGCTATAACTCACCTTCGCCGTCACCTCCCGCCCCGGGTTGTAGTAATTCGCCGTCCCCGACCCCCCCACATGCTGCTCGTCAAACAGGTTCCTCACATTCAGCGCCAGGTCGGTcc                                 >  1:319458/1‑136 (MQ=255)
                                                                                                       cTGGCTGCCTAACCAGCTATAACTCACCTTCGCCGTCACCTCCCGCCCCGGGTTGTAGTAATTCGCCGTCCCCGACCCCACCACATGCTGCTCGTCAAACAGCTGGCTCACATTCAGCGCCAGGTCGGTCCCTTTg                            >  5:67303/1‑136 (MQ=255)
                                                                                                          gCTGCCTTACCAGCTATAACTCACCTTCGCCTTCACCTCCCGCCCCGGGTTGTAGTAATTCGCCGTCCCCGACCCCACCACATGCTGCTCGTCAAACAGGTTGCTCACATTCAGCGCCAGGTCGGTCCCTTTgg                           <  7:306096/134‑1 (MQ=255)
                                                                                                          gCTGCCTTACCAGCTATAACTCACCTTCGCCGTCACCTCCCGCCCCGGGTTGTAGTAATTCGCCGTCCCCGACCCCACCACATGCTGCTCGTCAAACAGGTTGCTCACATTCAGCGCCAGGTCGGTCCCTTTGGCg                         <  3:397487/136‑1 (MQ=255)
                                                                                                             gCCTTACCAGCTATAACTCACCTTCGCCGTCACCTCCCGCCCCGGGTTGTAGTAATTCGCCGTCCCCCACCCCACCACAAGCTGCTCGTCAAAAAGGTTGCTCACATTCATCGCCAGGACGGGCCCTTTGGCGATc                      >  7:146290/1‑136 (MQ=255)
                                                                                                                 tACCAGCTATAACTCACCCTCGCCGTCACCTCCCGCCCCGGGTTGTAGTAATTCGCCGTCCCCGACCCCACCACATGCTGCTCGTCAAACAGGTTGCTCACATTCAGCGCCAGGTCGGTCCCTTTgg                           <  7:111552/127‑1 (MQ=255)
                                                                                                                 tACCAGCTATAACTCACCCTCGCCGTCACCTCCCGCCCCGGGTTGTAGTAATTCGCCGTCCCCGACCCCACCACATGCTGCTCGTCAAACAGGTTGCTCACATTCAGCGCCAGGTCGGTCCCTTTgg                           <  7:117217/127‑1 (MQ=255)
                                                                                                                 tACCAGCTATAACTCACCCTCGCCGTCACCTCCCGCCCCGGGTTGTAGTAATTCGCCGTCCCCGACCCCACCACATGCTGCTCGTCAAACAGGTTGCTCACATTCAGCGCCAGGTCGGTCCCTTTgg                           >  8:111552/1‑127 (MQ=255)
                                                                                                                 tACCAGCTATAACTCACCCTCGCCGTCACCTCCCGCCCCGGGTTGTAGTAATTCGCCGTCCCCGACCCCACCACATGCTGCTCGTCAAACAGGTTGCTCACATTCAGCGCCAGGTCGGTCCCTTTgg                           >  8:117217/1‑127 (MQ=255)
                                                                                                                  aCCAGCTATAACTCACCTTCGCCGTCACCTCCCGCCCCGGGTTGTAGTAATTCGCCGTCCCCGACCCCACCACATGCTGCTCGTCAAACAGGTTGCTCACATTCAGCGCCAGGTCGGTCCCTTTGGCGATCTTg                   >  3:22637/1‑134 (MQ=255)
                                                                                                                   ccAGCTATAACTCACCTTCGCCGTCACCTCCCGCCCCGGGTTGTAGTAATTCGCCGTCCCCGACCCCACCACATGCTGCTCGTCAAACAGGTTGCTCACATTCAGCGCCAGGTCGGTCCCTTTGGCGATCTTGTa                 <  2:247733/135‑1 (MQ=255)
                                                                                                                    cAGCTATAACTCACCTTCGCCGTCACCTCCCGCCCCGGGTTGTAGTAATTCGCCGTCCCCGACCCCACCACATGCTGCTCGTCAAACAGGTTGCTCACATTCAGCGCCAGGTCGGTCCCTTTGGCGATCTTGTAGt               >  8:274130/1‑136 (MQ=255)
                                                                                                                          taACTCACCTTCGCCGTCACCTCCCGCCCCGGGTTGTAGTAATTCGCCGTCCCCCACCCCACCACATGCTGCTCGTCAAACAGGTTGCTCACATTCAGCGCCAGGTCGGTCCCTTTGGCGATCTTGTAGTTCAgcg         >  7:90562/1‑136 (MQ=255)
                                                                                                                           aaCTCACCTTCGCCGTCACCTCCCGCCCCGGGTTGTAGTAATTCGCCGTCCCCGACCCCACCACATGCTGCTCGTCAAACAGGTTGCTCACATTCAGCGCCAGGTCGGTCCCTTTGGCGATCTTGTAGTTCAgcgc        >  5:128143/1‑136 (MQ=255)
                                                                                                                           aaCTCACCTTCGCCGTCACCTCCCGCCCCGGGTTGTAGTAATTCGCCGTCCCCGACCCCACCACATGCTGCTCGTCAAACAGGTTGCTCACATTCAGCGCCAGGTCGGTCCCTTTGGCGATCTTGTAGTTCAgcgc        >  1:105522/1‑136 (MQ=255)
                                                                                                                                 ccTTCGCCGTCACCTCCCGCCCCGGGTTGTAGTAATTCGCCGTCCCCGACCCCACCACATGCTGCTCGTCAAACAGGTTGCTCACATTCAGCGCCAGGTCGGTCCCTTTGGCGATCTTGTAGTTCAGCGCCGCATc  <  7:169274/136‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                     
TCAGTCTCTTTGAAGATACGGGGGGTTCATTGCAAGTTTCAGGATCTTGCGCACCTGTGGGCACGGGTTTACCCGCGAACAGGGGCGCAACCCGTGCCATCCCCTGGCTGCCTTACCAGCTATAACTCACCTTCGCCGTCACCTCCCGCCCCGGGTTGTAGTAATTCGCCGTCCCCGACCCCACCACATGCTGCTCGTCAAACAGGTTGCTCACATTCAGCGCCAGGTCGGTCCCTTTGGCGATCTTGTAGTTCAGCGCCGCATC  >  ACB143/5408849‑5409113

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑