Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,459,849 A→C 10.8% I79S (ATT→AGT)  bfr ← bacterioferritin, iron storage and detoxification protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,459,8490AC10.8% 97.9 / 8.2 37I79S (ATT→AGT) bfrbacterioferritin, iron storage and detoxification protein
Reads supporting (aligned to +/- strand):  ref base A (15/18);  new base C (2/2);  total (17/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.20e-01

TCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCACTCTTTGAA  >  CP009273/3459615‑3460114
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                         
tcttcGCGGATCTGTGCTTGCAGATGATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGACATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACCCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCATCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTc                                                                                                                                                                                                                                   <  1:676239/275‑1 (MQ=255)
  ttcGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACCTCTTCCCCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAACCCTTcc                                                                                                                                                                                                                                  >  1:30123/1‑274 (MQ=255)
            tgtgCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTGTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCATGTAAGTTTGGAAGACCTTCCAGAAAAAGaa                                                                                                                                                                                                                        >  1:149716/1‑274 (MQ=255)
                 ttgacgttatgtttGAAGGCCCATCTTCTGAATCAAATCAAGTTCCGTTTCCAGCCAGTCGATATTGCCTTCTTCATCACGCAAAATTTGTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCt                                                                                                                                                                                                                  <  2:149715/264‑1 (MQ=255)
                 ttgaccaTAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCt                                                                                                                                                                                                                  <  2:30123/269‑1 (MQ=255)
                                             tGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCACTg                                                                                                                                                                                                                                                                         >  1:182113/1‑192 (MQ=255)
                                             tGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCACTg                                                                                                                                                                                                                                                                         >  1:606949/1‑192 (MQ=255)
                                               aaTCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCANATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTtcatcccat                                                                                                                                                                                    >  1:600881/1‑275 (MQ=255)
                                                 tCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCACTGtgaa                                                                                                                                                                                                                                                                     <  2:606940/192‑4 (MQ=255)
                                                 tCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCACTGtgaa                                                                                                                                                                                                                                                                     <  2:182110/192‑4 (MQ=255)
                                                                          gTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGaa                                                                                                                                                                                                                        >  1:206889/1‑212 (MQ=255)
                                                                              atatGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATgcg                                                                                                                                                                                                                    <  2:206885/212‑1 (MQ=255)
                                                                                     cttcttCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAg                                                                                                                                                                                                                                >  2:517075/1‑193 (MQ=255)
                                                                                         ttcATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGaaaa                                                                                                                                                                                                                            <  1:517082/193‑1 (MQ=255)
                                                                                            aaaaagCAAAATTTATATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGTTTGAGAcccc                                                                                                                                       <  2:600872/270‑1 (MQ=255)
                                                                                             tcaCGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATg                                                                                                                                                                                                                      >  2:318841/1‑195 (MQ=255)
                                                                                                 gCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCt                                                                                                                                                                                                                  <  1:318847/195‑1 (MQ=255)
                                                                                                                 catATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTc                                                                                                                                                                                                                                                                      >  2:628797/1‑127 (MQ=255)
                                                                                                                     tCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGtt                                                                                                                                                                                                                                                                  <  1:628806/127‑1 (MQ=255)
                                                                                                                                        gAACGCTATCGGCATTACCAAGTGTCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATg                                                                                                                                                                                                                      >  1:492653/1‑152 (MQ=255)
                                                                                                                                         aaCGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAAt                                                                                                                                                                                >  1:167313/1‑189 (MQ=255)
                                                                                                                                          aCGCTATCGGCATAACCAATTTCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTNAACATCTTNACCAATGTTCATATTGCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACCTTCGGGCCTGGAGAAAGTACTGATTGATTGCGAc                                                                                         <  1:464300/275‑1 (MQ=255)
                                                                                                                                           cGCTATCGGCATAACCAATTGCCGCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATTCGCTCAATATAACGATCGGCGTGTCTCATCTCATNAATGGATTCATGATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGTGCATGGAGAAAGTACTGATTGATTGCGAc                                                                                         >  1:487930/1‑274 (MQ=255)
                                                                                                                                            gCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCt                                                                                                                                                                                                                  <  2:492646/152‑1 (MQ=255)
                                                                                                                                             cTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGAt                                                                                                                                                                            <  2:167310/189‑1 (MQ=255)
                                                                                                                                                  ggCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATg                                                                                                                                                                                                                                                                         <  2:135646/91‑1 (MQ=255)
                                                                                                                                                            aTTGCCGCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTCCCTCAACATCTTCACCAATGTTCAgat                                                                                                                                                                                                                                                                  >  2:138798/1‑86 (MQ=255)
                                                                                                                                                                ccTCACGCAAATCCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACCTCTTCACCAATGTTCAGTTTGcc                                                                                                                                                                                                                                                              <  1:138799/88‑1 (MQ=255)
                                                                                                                                                                  tCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTttgt                                                                  >  1:691632/1‑274 (MQ=255)
                                                                                                                                                                      gCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATATTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTgaga                                                             <  2:691623/275‑1 (MQ=255)
                                                                                                                                                                          aTTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGtt                                                                                                                                                                                                                                                                  >  1:210993/1‑74 (MQ=255)
                                                                                                                                                                              ttcGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGcc                                                                                                                                                                                                                                                              <  2:210989/74‑1 (MQ=255)
                                                                                                                                                                                                                  agcaTTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATg                                                                                                                                                                                                                      >  1:249370/1‑78 (MQ=255)
                                                                                                                                                                                                                      tttCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCt                                                                                                                                                                                                                  <  2:249365/78‑1 (MQ=255)
                                                                                                                                                                                                                       ttCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCt                                                                                                                                                                                                                                                      >  1:506451/1‑41 (MQ=255)
                                                                                                                                                                                                                           tCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTaa                                                                                                                                                                                                                                                  <  2:506444/41‑1 (MQ=255)
                                                                                                                                                                                                                                   gtCACAAATGTTTAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTACAGAAATAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGGTTCATGATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACCGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCACTCTTTGaa  <  2:487923/272‑1 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                         
TCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCAAGTGCCAGATCAGAACGCAGCATTTCCTCAACATCTTCACCAATGTTCAGTTTGCCCAGGTCCTGTAAGTTTGGAAGACCTTCCAGAAAAAGAATGCGCTCAATATAACGATCGGCGTGTTTCATCTCATCAATGGATTCATGATACTCCACATCATTGAGACGTTTGAGACCCCAGTTTTTAAACATTCGGGCATGGAGAAAGTACTGATTGATTGCGACAAGCTCATTTCCCAACAGTTTGTTGAGATAATTTATAACTTTAGTATCACCTTTCATTTTTTGACTCCTTCGCTTCCACTCTTTGAA  >  CP009273/3459615‑3460114

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: