Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,466,313 A→G L138P (CTT→CCT)  bfr ← bacterioferritin, iron storage and detoxification protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,466,3130AG98.2% 191.7 / ‑3.8 54L138P (CTT→CCT) bfrbacterioferritin, iron storage and detoxification protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (29/24);  minor base C (0/1);  total (29/25)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

AGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCA  >  NC_000913/3466179‑3466447
                                                                                                                                      |                                                                                                                                      
agcaTGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTc                                                                                                                                    <  1:875960/139‑1 (MQ=255)
            cTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTccagcca                                                                                                                        <  1:503201/139‑1 (MQ=255)
             tGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTccagccag                                                                                                                       <  1:467892/139‑1 (MQ=255)
               gggTATCTGGCCTGGATCTGTTATTAACCAGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCAGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTc                                                                                                                     >  2:813494/1‑139 (MQ=255)
                   aTCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGata                                                                                                                 <  2:524614/139‑1 (MQ=255)
                   aTCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGata                                                                                                                 <  2:533670/139‑1 (MQ=255)
                    tCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGatat                                                                                                                >  1:112235/1‑139 (MQ=255)
                    tCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGatat                                                                                                                >  1:579782/1‑139 (MQ=255)
                                      ttaACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCa                                                                                              <  2:420204/139‑1 (MQ=255)
                                        aaCCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCaaa                                                                                            >  1:299813/1‑139 (MQ=255)
                                         aCCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCaaaa                                                                                           <  2:1050744/139‑1 (MQ=255)
                                         aCCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCaaaa                                                                                           >  2:771425/1‑139 (MQ=255)
                                           cGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAAtt                                                                                         >  2:89236/1‑139 (MQ=255)
                                            gggAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAAttt                                                                                        >  2:165516/1‑139 (MQ=255)
                                             ggAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTccagccag                                                                                                                       >  1:1703/1‑107 (MQ=255)
                                             ggAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTccagccag                                                                                                                       <  2:1703/107‑1 (MQ=255)
                                             ggAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTccagc                                                                                                                          >  1:24137/1‑104 (MQ=255)
                                             ggAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTccagc                                                                                                                          <  2:24137/104‑1 (MQ=255)
                                             ggAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTc                                                                                       >  2:847841/1‑139 (MQ=255)
                                                   ttCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTc                                                                                       <  2:494623/133‑1 (MQ=255)
                                                   ttCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTc                                                                                       >  1:494623/1‑133 (MQ=255)
                                                   ttCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTatcat                                                                                 >  2:929343/1‑139 (MQ=255)
                                                        cctcccGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCatat                                                                            >  2:776914/1‑139 (MQ=255)
                                                            ccGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATcgcg                                                                        >  1:87161/1‑139 (MQ=255)
                                                              gACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGc                                                                      <  2:579782/139‑1 (MQ=255)
                                                                              cttcGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGtt                                                                                                                                     >  2:971010/1‑60 (MQ=255)
                                                                              cttcGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGtt                                                                                                                                     <  1:971010/60‑1 (MQ=255)
                                                                                 cgcgGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGata                                                                                                                 <  2:838130/77‑1 (MQ=255)
                                                                                 cgcgGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGata                                                                                                                 >  1:838130/1‑77 (MQ=255)
                                                                                   cgGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCACGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGCAATCATGAACGCTAt                                                 <  1:813494/139‑1 (MQ=255)
                                                                                        cTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATc                                                            <  1:1206052/123‑1 (MQ=255)
                                                                                        cTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATc                                                            >  2:1206052/1‑123 (MQ=255)
                                                                                         tgtgCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCAt                                           <  1:89236/139‑1 (MQ=255)
                                                                                                gCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTccagccag                                                                                                                       >  1:1157999/1‑56 (MQ=255)
                                                                                                gCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTccagccag                                                                                                                       <  2:1157999/56‑1 (MQ=255)
                                                                                                gCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTccagcc                                                                                                                         <  1:710933/54‑1 (MQ=255)
                                                                                                gCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTccagcc                                                                                                                         >  2:710933/1‑54 (MQ=255)
                                                                                                gCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACg                                                                                                <  1:447588/79‑1 (MQ=255)
                                                                                                gCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACg                                                                                                >  2:447588/1‑79 (MQ=255)
                                                                                                 cAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAAtt                                   >  1:1079492/1‑139 (MQ=255)
                                                                                                  aGATAATTTTGCAGGCCCGTCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCACGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTg                                  >  1:768599/1‑139 (MQ=255)
                                                                                                  aGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATgg                                                                                                              >  2:994745/1‑63 (MQ=255)
                                                                                                  aGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATgg                                                                                                              <  1:994745/63‑1 (MQ=255)
                                                                                                  aGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTg                                  >  1:749353/1‑139 (MQ=255)
                                                                                                   gATAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGc                                 >  2:1110652/1‑139 (MQ=255)
                                                                                                     tAATTTTGCAGGCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCt                               <  2:112235/139‑1 (MQ=255)
                                                                                                                gCCCATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAAttc                    >  1:1170313/1‑139 (MQ=255)
                                                                                                                 cccATCTTCTGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAAttct                   <  1:847841/139‑1 (MQ=255)
                                                                                                                      cttctGAATCAGATCAGGTTCCGTTTTCAGCCAGTCGCTATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTcgcg              <  1:776914/139‑1 (MQ=255)
                                                                                                                      cttctGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAATTTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTcgcg              <  1:771425/139‑1 (MQ=255)
                                                                                                                          tGAATCAGATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCAt          >  1:295336/1‑139 (MQ=255)
                                                                                                                                 gATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGt                                                                >  2:348564/1‑78 (MQ=255)
                                                                                                                                 gATCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGt                                                                <  1:348564/78‑1 (MQ=255)
                                                                                                                                  aTCAGGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCa  >  1:72686/1‑139 (MQ=255)
                                                                                                                                  aTCAGGGTCCCTTTCCAAGCAGTCGATTTGGCCCTCCTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGGAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCa  <  2:87161/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                      
AGCATGATAAACCTGGGGTATCTGGCCTGGATCTCTTATTAACCGGGAGGGTTCTCCCTCCCGACACGGCTCAACCTTCTTCGCGGATCTGTGCTTGCAGATAATTTTGCAGGCCCATCTTCTGAATCAGATCAAGTTCCGTTTCCAGCCAGTCGATATGGCCTTCTTCATCACGCAAAATTTCTATCATCATATCGCGGCTGACGTAATCATGAACGCTATCGGCATAACCAATTGCCTCACGCAAATTCTTCGCGCCATCCAGCTCA  >  NC_000913/3466179‑3466447

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