Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,605,176 T→C 12.5% Q163R (CAG→CGG)  fhuF ← hydroxamate siderophore iron reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,605,1760TC12.5% 68.3 / 2.3 33Q163R (CAG→CGG) fhuFhydroxamate siderophore iron reductase
Reads supporting (aligned to +/- strand):  ref base T (13/14);  new base C (2/2);  total (16/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.22e-01

GATTCAACAGTAGCCTCGCCGAGCAGTTGCTTCATCTCAGTGAGATACCAGTTGATCAAATAACCGGTATTACTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGACCCAGAAACAGGCGACGCGTCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGACACATCTAATGC  >  NC_000913/4605038‑4605315
                                                                                                                                          |                                                                                                                                           
gATTCAACAGTAGCCTCGCCGAGCAGTTGCTTCATCTCAGTGAGATACCAGTTGATCAAATAACCGGTATTACTCCAGATAAGTTTGCCGTTGATCCCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCt                                                                                                                                         <  1:432823/143‑1 (MQ=255)
   tCAACAGTAGCCTCGCCGAGCAGTTGCTTCATCTCAGTGAGATACCAGTTGATCAAATAACCGGGATTACTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCGTCTAGTGCTTGCACAACCGGACCCAGCGCCTGGCGGAt                                                                                                                                      >  2:139888/1‑143 (MQ=255)
            gCCTCGCCGAGCAGTTGCTTCATCTCAGTGAGCTACCAGTTGATCAACTAACCGGTATTACTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTc                                                                                                                             <  1:387647/143‑1 (MQ=255)
             ccTCGCCGAGCAGTTGCTTCATCTCAGTGAGATACCAGTTGATCAAATAACCGGTATTACTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTcc                                                                                                                            <  1:524428/143‑1 (MQ=255)
             ccTCGCCGAGCAGTTGCTTCATCTCAGTGAGATACCAGTTGATCAAATAACCGGTATTACTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTcc                                                                                                                            <  1:523297/143‑1 (MQ=255)
                cggcGAGCAGTTGCTTCATCTCAGTGAGATACCAGTTGATCAAATAACCGGTATTACTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCAtt                                                                                                                         <  1:203702/140‑1 (MQ=255)
                  ccGAGCAGTTGCTTCATCTCAGTGAGATACCAGTTGATCAAATAACCGGTATTACTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCGTCTAGTGCTTGCACAACCGGAAAGTGCGACAGGCTGATTAACGTGTCCATTCg                                                                                                                       >  2:548338/1‑143 (MQ=255)
                          ttGCTTCATCTCAGTGAGATACCAGTTGATCAAATAACCGGTATTACTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTgcg                                                                                                               <  1:280720/143‑1 (MQ=255)
                             cTTCATCTCAGTGAGATACCAGTTGATCAAATAACCGGTATTACTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTACCGTTTCCATTCGATGCTgcggcg                                                                                                            >  1:227036/1‑143 (MQ=255)
                              ttCATCTCAGTGAGATACCAGTTGATCAAATAACCGGTATTACTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGa                                                                                                           >  1:402611/1‑143 (MQ=255)
                                      aGTGAGATACCAGTTGATCAAATAACCGGTATTACTCCAGATAAGCATGCCGTTGAACTCTCCGGGCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATCAGATGCTGCGGCGAATgtggtg                                                                                                   >  2:243102/1‑143 (MQ=255)
                                              aCCAGTTGATCAAATAACCGGTATTACTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGCACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGttt                                                                                           >  2:482394/1‑143 (MQ=255)
                                               ccAGTTGATCAAATAACCGGTATTACTCCAGAGACCGTTCCATTTTATTTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGTCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGtttt                                                                                          <  2:118853/143‑1 (MQ=255)
                                                 aGTTGATCAAATAACCGGTATTACTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAGTGTGGTGTTGCGTTTGTa                                                                                        >  1:195315/1‑143 (MQ=255)
                                                      aTCAAATAACCGGTATTACTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACAGTTCCATGCGACGCTACGGCGAGTGTGGTGTTGCGTTTTTATCGTc                                                                                   >  2:273490/1‑143 (MQ=255)
                                                            tAACCGGTATTACTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTGTAGTGATTGCACAACCGGAACCAGCGCCTGGCTGATGAACGTTTCCCTCCGATGCTGAGCCGCAGGGGGTGTTGCGGTTTGATCTTcccacac                                                                             >  2:185680/1‑143 (MQ=255)
                                                              aCCGGTATTACTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTCTATCTTCACTCTCAt                                                                           >  1:286845/1‑143 (MQ=255)
                                                               ccGGTATTACTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCCCACACATc                                                                          >  1:122198/1‑143 (MQ=255)
                                                               ccGGTATTACTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTGTTTATCTTCACACACCTc                                                                          >  1:378148/1‑143 (MQ=255)
                                                                 cgTATTACTCCAGCTAAGTGGGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCTGCGGCTGGCTGATTAACGTTTCCTTTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGa                                                                        <  2:463345/142‑1 (MQ=255)
                                                                    aTTACTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGAccc                                                                     <  1:420959/143‑1 (MQ=255)
                                                                       aCTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGACCCAGa                                                                  <  2:180955/143‑1 (MQ=255)
                                                                                       ccGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGACCCAGAAACAGGCGACGCGTcc                                                  >  1:165473/1‑143 (MQ=255)
                                                                                        cGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCCGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGACCCAGaaa                                                                <  1:561122/128‑1 (MQ=255)
                                                                                        cGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCCGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGACCCAGaaa                                                                <  1:410555/128‑1 (MQ=255)
                                                                                        cGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCCGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGACCCAGaaa                                                                >  2:410555/1‑128 (MQ=255)
                                                                                        cGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCCGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGACCCAGaaa                                                                >  2:561122/1‑128 (MQ=255)
                                                                                                cctcCGGTCGCGTCTCGTGCGTGCGCAACCGGAACCAGCGCTTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACAGCGACCCAGAAACAGGCGACGCGTCCGGTTTCGTg                                         <  2:214497/142‑1 (MQ=255)
                                                                                                 ctcCGGTCGCTTCTAGTGCTGGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGACCCAGAAACAGGCGACGCGTCCGGTTTCGTGa                                        <  1:326386/143‑1 (MQ=255)
                                                                                                        cgcgTCTAGTGCTTGCACAAGCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGATGGTGGTGTGGCGTTTTGATCTTCAGACACATCGACCCAGAAACAGGCGACGCGTCCGGTTTCGTGAAACTCAg                                 <  2:121624/143‑1 (MQ=255)
                                                                                                               aGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGCTTTTATCTTCACACCCATCGACCCAGCAACAGGCGACGCCTCCCGTCTCGTGCCACCCAGCACGGaa                          >  2:207394/1‑143 (MQ=255)
                                                                                                                           aacGGAACGAGAGCCTGGCTGATTAACGCTTCCATGCGATGCTGCGGCGAATGTGGTGTTGCGGTTTTATCTTCACACACATCGACCCAGAAACAGGCGACGCGTCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGa              <  1:19441/141‑1 (MQ=255)
                                                                                                                                       gcCGGGCTGATTAACGTTCCCATTCGCTGCTGCGGCGAAGCTGGTGTTGCGCTTTTATCTTCACACACATCGATCCACAAACAGGCGTCGCGTCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGACACATCTAATGc  <  1:361827/143‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                           
GATTCAACAGTAGCCTCGCCGAGCAGTTGCTTCATCTCAGTGAGATACCAGTTGATCAAATAACCGGTATTACTCCAGATAAGTTTGCCGTTGATCTCTCCGGTCGCTTCTAGTGCTTGCACAACCGGAACCAGCGCCTGGCTGATTAACGTTTCCATTCGATGCTGCGGCGAATGTGGTGTTGCGTTTTTATCTTCACACACATCGACCCAGAAACAGGCGACGCGTCCGGTTTCGTGAAACTCAGCATGGAAATGTTCCGGCGACACATCTAATGC  >  NC_000913/4605038‑4605315

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