Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,501,175 T→C I34V (ATC→GTC)  fecE ← Fe(3+) dicitrate ABC transporter ATP‑binding protein FecE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,501,1750TC96.3% 89.2 / ‑5.8 27I34V (ATC→GTC) fecEFe(3+) dicitrate ABC transporter ATP‑binding protein FecE
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (14/12);  total (15/12)
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.00e+00

TGAGAGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACA  >  NZ_CP009273/4501065‑4501302
                                                                                                              |                                                                                                                               
tGAGAGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGACCAGGGCGGTGATCTTCCCCGTTGGCAGt                                                                                                     >  2:268982/1‑139 (MQ=255)
tGAGAGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGACCAGGGCGGTGATCTTCCCCGTTGGCAGt                                                                                                     <  1:144689/139‑1 (MQ=255)
             aTGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGACCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAAcg                                                                                         <  1:165774/138‑1 (MQ=255)
              tGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGACCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAAcgtc                                                                                       >  1:384370/1‑139 (MQ=255)
                  aTTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGACCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTa                                                                                   <  1:484495/139‑1 (MQ=255)
                   ttATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGACCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGtt                                                                                    >  1:71147/1‑137 (MQ=255)
                    tATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGACCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTa                                                                                   <  2:71147/137‑1 (MQ=255)
                     aTCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGACCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAg                                                                                 >  2:162875/1‑138 (MQ=255)
                     aTCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGACCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAg                                                                                 >  1:507734/1‑138 (MQ=255)
                     aTCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGACCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAg                                                                                 >  1:563429/1‑138 (MQ=255)
                        gCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGACCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTAc                                                                              <  2:45864/138‑1 (MQ=255)
                            agaAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGACCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTtgt                                                                         >  2:61766/1‑139 (MQ=255)
                                      tGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGACCa                                                                                                                               <  1:53813/75‑1 (MQ=255)
                                      tGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGACCa                                                                                                                               >  2:53813/1‑75 (MQ=255)
                                        ccAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGACCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAATGTCGtt                                                                                    <  2:200489/116‑1 (MQ=255)
                                        ccAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGACCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAATGTCGtt                                                                                    >  1:200489/1‑116 (MQ=255)
                                        ccAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGACCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAAc                                                             >  2:139377/1‑139 (MQ=255)
                                                cGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGACCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAAcg                                                                                         <  1:308709/103‑1 (MQ=255)
                                                cGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGACCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAAcg                                                                                         >  2:308709/1‑103 (MQ=255)
                                                cGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGACCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGt                                                      <  1:193885/138‑1 (MQ=255)
                                                cGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGACCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGt                                                      >  2:440642/1‑138 (MQ=255)
                                                        aaaGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGACCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTc                                             <  1:268982/139‑1 (MQ=255)
                                                                                  gTCGATTTCCCGCAACCGTTAGGACCGACCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATtc                   <  1:69641/139‑1 (MQ=255)
                                                                                   tCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATtct                  >  1:254063/1‑139 (MQ=255)
                                                                                   tCGATTTCCCGCAACCGTTAGGACCGACCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTAGTAAAGTCATTTCTCGCATtct                  >  2:422083/1‑139 (MQ=255)
                                                                                                 ccGTTAGAACCGACCATGGCGGTGATCTTCCGCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTGCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGa    <  2:328768/139‑1 (MQ=255)
                                                                                                    ttAGGACCGACCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACa  >  1:165332/1‑138 (MQ=255)
                                                                                                              |                                                                                                                               
TGAGAGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGATCAGGGCGGTGATCTTCCCCGTTGGCAGTGAGAGTGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACA  >  NZ_CP009273/4501065‑4501302

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