Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,541,397 T→C F331L (TTC→CTC)  feoB → Fe(2(+)) transporter FeoB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,541,3970TC100.0% 142.8 / NA 41F331L (TTC→CTC) feoBFe(2(+)) transporter FeoB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (19/22);  total (19/22)

TTTCCTCTTTGTGATGTACCTGATGTTCCTGCTGGCTATCAACATCGGCGGGGCGTTACAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACAC‑GCTCCACTTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTG  >  NC_000913/3541267‑3541532
                                                                                                                                   |                                                                                                                                       
tttCCTCTTTGTGATGTACCTGATGTTCCTGCTGGCTATCAACATCGGCGGGGCGTTACAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACAC‑GCTCCACCTCCCGGACt                                                                                                                                <  1:545797/140‑1 (MQ=255)
          gtgATGTACCTGATGTTCCTGCTGGCTATCAACATCGGCGGGGCGGTACAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGTTATTCAATGGATAGGCTACAC‑GCTCACCCTACCGGGCTGGCTGACTAt                                                                                                                      >  2:267213/1‑140 (MQ=255)
              tGTACCTGATGTTCCTGCTGGCTATCAACATCGGCGGGGCGTTACAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACAC‑GCTCCACCTCCCGGACTGGCTGACTATCTTc                                                                                                                  <  1:8391/140‑1 (MQ=255)
              tGTACCTGATGTTCCTGCTGGCTATCAACATCGGCGGGGCGTTACAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACAC‑GCTCCACCTCCCGGACTGGCTGACTATCTTc                                                                                                                  <  1:106094/140‑1 (MQ=255)
                tACCTGATGTTCCTGCTGGCTATCAACATCGGCGGGGCGTTACAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACAC‑GCTCCACCTCCCGGACTGGCTGACTATCTTCCt                                                                                                                <  1:610270/140‑1 (MQ=255)
                   cTGATGTTCCTGCTGGCTATCAACATCGGCGGGGCGTTACAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGGATTCAATGGATTTGCTACAC‑GCTCCACCTCCCGGACTGGCTGACTATCTTCCTCGc                                                                                                             >  1:445489/1‑140 (MQ=255)
                    tGATGTTCCTGCTGGCTATCAACATCGGCGGGGCGTTACAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACAC‑GCTCCACCTCCCGGACTGGCTGACTATCTTCCTCGcc                                                                                                            <  2:364739/140‑1 (MQ=255)
                     gATGTTCCTGCTGGCTATCAACATCGGCGGGGCGTTACAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACAC‑GCTCCACCTCCCGGACTGGCTGACTATCTTCCTCGccc                                                                                                           >  2:409008/1‑140 (MQ=255)
                          tCCTGCTGGCTATCAACATCGGCGGGGCGTTACAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACAC‑GCTCCACCTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGc                                                                                                      >  2:538119/1‑140 (MQ=255)
                           ccTGCTGGCTATCAACATCGGCGGGGCGTTACAGCCGCTGTTTTACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACAC‑GCTCCACCTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGcc                                                                                                     >  2:9985/1‑140 (MQ=255)
                            ctgctgGCTATCAACATCGGCGGGGCGTTACAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACAC‑GCTCCACCTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCt                                                                                                    >  1:595406/1‑140 (MQ=255)
                                 ggCTATCAACATCGGCGGGGCGTTACAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACAC‑GCTCCACCTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTg                                                                                               <  1:634263/140‑1 (MQ=255)
                                 ggCTATCAACATCGGCGGGGCGGTACAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACAC‑GTTCCCCCTCTCGGCCTGGCTGACTATATTCCTCGCCCCGGGGCTggggg                                                                                               >  2:263211/1‑138 (MQ=255)
                                   cTATCAACATCGGCGGGGCGTTACAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACAC‑GCTCCACCTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTggc                                                                                             <  1:450837/140‑1 (MQ=255)
                                            tCGGCGGGGCGTTACAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACAC‑GCTCCCCCTCCCGGACTGGCTAGCTGTCTTTCTCGCCCCGGGCGCGGGGGGCGGCATTAac                                                                                    >  2:281392/1‑140 (MQ=255)
                                             cggcggGGCGTTACAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACAC‑GCTCCACCTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAaca                                                                                   >  1:656711/1‑140 (MQ=255)
                                              ggcggGGCGTTACAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACAC‑GCTCCACCTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAacac                                                                                  <  2:350471/140‑1 (MQ=255)
                                                 ggggCGTTACAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACAC‑GCTCCACCTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGt                                                                               <  1:561181/140‑1 (MQ=255)
                                                          cAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACAC‑GCTCCACCTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACt                                                                      >  1:508726/1‑140 (MQ=255)
                                                                         gTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACAC‑GCTCCACCTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTgg                                                       <  1:9985/140‑1 (MQ=255)
                                                                         gTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACAC‑GCTCCACCTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTgg                                                       <  1:538119/140‑1 (MQ=255)
                                                                         gTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACAC‑GCTACACCTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTgg                                                       <  2:623238/140‑1 (MQ=255)
                                                                         gTCGGCTCCGTGGCGATATTTGGGCGTGGTATTCAATGGATTGGCTACAC‑GCTCCACCTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTgg                                                       <  1:306689/140‑1 (MQ=255)
                                                                           ccgttCCGTAGCGGTAACTGTGCATGGGTTTGAATCGATATGCTACACAGCCCTACTT‑CCGGACTGGCTGACTATCTTCCTCCGCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTGGCa                                                     <  2:282810/136‑1 (MQ=255)
                                                                           cGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACAC‑GCTCCACCTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTGGCa                                                     <  2:508726/140‑1 (MQ=255)
                                                                           cGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACAC‑GCTCCACCTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTGGCa                                                     <  2:161712/140‑1 (MQ=255)
                                                                           cGGCTCCGTGGCGCTATTTGTGAATGGTATTCATTCGATTGGCTACAC‑TCTCCACCTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTGGCa                                                     <  2:288436/140‑1 (MQ=255)
                                                                             gCTCCGTGGCGCTATTTGTGCATGCTATTCAATGGATTGGCTACAC‑GCTCCACCTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTGGCatg                                                   >  2:444381/1‑140 (MQ=255)
                                                                                      cgcTATTTGTGCATGGTATTCAATGGATTGGCTACAC‑GCTCCACCTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCTGCAGATTGGCATGATGTACCTg                                          <  2:595406/140‑1 (MQ=255)
                                                                                                gCATGGTATTCAATGGATTGGCTACAC‑GCTCCACCTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCtcct                                >  2:28850/1‑140 (MQ=255)
                                                                                                        ttCAATGGATTGGCTACAC‑GCTCCACCTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAg                        <  1:35823/140‑1 (MQ=255)
                                                                                                        ttCAATGGATTGGCTACAC‑GCTCCACCTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAg                        <  2:445489/140‑1 (MQ=255)
                                                                                                           aaTGGATTGGCTACAC‑GCTCCACCTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGAc                     <  2:497518/140‑1 (MQ=255)
                                                                                                             tGGATTGGCTACAC‑GCTCCACCTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTc                   <  1:409008/140‑1 (MQ=255)
                                                                                                              ggATTGGCTACAC‑GCTCCACCTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTcc                  >  2:7522/1‑140 (MQ=255)
                                                                                                                aTTGGCTACAC‑GCTCCACCTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTTGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCgg                >  2:500401/1‑140 (MQ=255)
                                                                                                                aTTGGCTACAC‑GCTCCACCTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTAGCTGGTCATCTCCTTCCTTGAGGACTCCgg                >  1:216392/1‑140 (MQ=255)
                                                                                                                aTTGGCTACAC‑GCTCCACCTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCgg                >  1:142095/1‑140 (MQ=255)
                                                                                                                aTTGGCTACAC‑GCTCCACCTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCgg                >  1:56595/1‑140 (MQ=255)
                                                                                                                aTTGGCTACAC‑GCTCCACCTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCgg                >  1:335885/1‑140 (MQ=255)
                                                                                                                aTTGGCTACAC‑GCTCCACCTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTAACTGGTCGTCTCCTTCCTTGAGGACTCCgg                >  1:235252/1‑140 (MQ=255)
                                                                                                                    gCTACAC‑GCTCCACCTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCtccttc                              >  2:407172/1‑122 (MQ=255)
                                                                                                                    gCTACAC‑GCTCCACCTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCtccttc                              <  1:407172/122‑1 (MQ=255)
                                                                                                                         ac‑GCTCCACCTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGgc       >  2:481283/1‑140 (MQ=255)
                                                                                                                         ac‑GCTCCACCTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGgc       >  1:328878/1‑140 (MQ=255)
                                                                                                                               ccACCTCCCGGACTGACTGACTCTCATCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTg  <  1:281392/140‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                       
TTTCCTCTTTGTGATGTACCTGATGTTCCTGCTGGCTATCAACATCGGCGGGGCGTTACAGCCGCTGTTTGACGTCGGCTCCGTGGCGCTATTTGTGCATGGTATTCAATGGATTGGCTACAC‑GCTCCACTTCCCGGACTGGCTGACTATCTTCCTCGCCCAGGGCCTGGGTGGCGGCATTAACACCGTGCTGCCACTGGTGCCGCAGATTGGCATGATGTACCTGTTCCTCTCCTTCCTTGAGGACTCCGGGTATATGGCGCGTG  >  NC_000913/3541267‑3541532

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