Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,540,301 G→A D47N (GAC→AAC)  feoA → ferrous iron transport protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,540,3010GA98.4% 207.2 / NA 62D47N (GAC→AAC) feoAferrous iron transport protein A
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (28/33);  total (29/33)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGCAAAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGAAAAAATTAACCATTGGCTTA  >  NC_000913/3540171‑3540430
                                                                                                                                  |                                                                                                                                 
cacTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCAAc                                                                                                                                 <  1:600687/133‑1 (MQ=255)
   tCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTTATGCGGTACGCGTCGCTCCACCCGGCAAcccc                                                                                                                              >  2:204950/1‑133 (MQ=255)
   tCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCAAcccc                                                                                                                              >  1:537077/1‑133 (MQ=255)
   tCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCAAcccc                                                                                                                              >  1:335986/1‑133 (MQ=255)
   tCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGCGGTTCGCGTCGCTCCACACGGCAAcccc                                                                                                                              >  2:219339/1‑133 (MQ=255)
     cAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCAACCCCAt                                                                                                                            >  1:487379/1‑133 (MQ=255)
     cAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCAACCCCAt                                                                                                                            >  1:133105/1‑133 (MQ=255)
       gATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCAACCCCATTc                                                                                                                          <  2:42731/133‑1 (MQ=255)
        aTACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCAACCCCATTCa                                                                                                                         >  1:438613/1‑133 (MQ=255)
            tGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCAACCCCATTCATATc                                                                                                                     <  2:108207/133‑1 (MQ=255)
             gCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCGCTCAGCGACCCCAATCATATCg                                                                                                                    >  1:167699/1‑133 (MQ=255)
             gCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCAACCCCATTCATATCg                                                                                                                    >  1:351401/1‑133 (MQ=255)
             gCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCAACCCCATTCATATCg                                                                                                                    >  2:91418/1‑133 (MQ=255)
              cGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCAACCCCATTCATATCGa                                                                                                                   >  1:438872/1‑133 (MQ=255)
              cGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCAACCCCATTCATATCGa                                                                                                                   >  1:519850/1‑133 (MQ=255)
               gTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCAACCCCATTCATATCGaa                                                                                                                  >  1:69690/1‑133 (MQ=255)
                 ggAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCAACCCCATTCATATCGAAAc                                                                                                                <  1:91418/133‑1 (MQ=255)
                    aaaTCACTGGCTTTTCCCGTGAAATCAGCTCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCAACCCCATTCATATCGAAACCcg                                                                                                             <  2:329890/133‑1 (MQ=255)
                    aaaTCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCAACCCCATTCATATCGAAACCcg                                                                                                             <  1:117579/133‑1 (MQ=255)
                    aaaTCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCAACCCCATTCATATCGAAACCcg                                                                                                             <  2:40623/133‑1 (MQ=255)
                    aaaTCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCAACCCCATTCATATCGAAACCcg                                                                                                             <  2:80021/133‑1 (MQ=255)
                      aTCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCAACCCCATTGATATCGAAACCcgtc                                                                                                           >  1:229027/1‑133 (MQ=255)
                      aTCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCAACCCCATTCATATCGAAACCcgtc                                                                                                           >  2:587993/1‑133 (MQ=255)
                         aCTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCAACCCCATTCATATCGAAACCCGTCgtg                                                                                                        <  1:135439/133‑1 (MQ=255)
                          cTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCAACCCCATTCATATCGAAACCCGTCgtgt                                                                                                       <  2:351401/133‑1 (MQ=255)
                          cTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCAACCCCATTCATATCGAAACCCGTCgtgt                                                                                                       <  2:487379/133‑1 (MQ=255)
                          cTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCAACCCCATTCATATCGAAACCCGTCgtgt                                                                                                       <  1:557218/133‑1 (MQ=255)
                          cTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCAACCCCATTCATATCGAAACCCGTCgtgt                                                                                                       <  2:537077/133‑1 (MQ=255)
                           tGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCAACCCCATTCATATCGAAACCCGTCGtgtg                                                                                                      <  2:133105/133‑1 (MQ=255)
                             gCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCAACCCCATTCATATCGAAACCCGTCgtgt                                                                                                       <  2:357061/130‑1 (MQ=255)
                             gCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCAACCCCATTCATATCGAAACCCGTCgtgt                                                                                                       >  1:357061/1‑130 (MQ=255)
                                       tGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCAACCCCATTCATATCGAAACCCGACGAGTGAGCTTGGTATTa                                                                                          >  1:246860/1‑133 (MQ=255)
                                                     cATATCGCTAAAAACTGCTTTCGCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCAACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGCAAAAAAGATCt                                                                            <  1:280307/133‑1 (MQ=255)
                                                      atatCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCAACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGCAAAAAAGATCTg                                                                           <  2:438613/133‑1 (MQ=255)
                                                        atCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCAACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGCAAAAAAGATCTGGc                                                                         <  2:438872/133‑1 (MQ=255)
                                                               aaaaCTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCAACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGCAAAAAAGATCTGGc                                                                         >  2:626563/1‑126 (MQ=255)
                                                               aaaaCTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCAACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGCAAAAAAGATCTGGc                                                                         <  1:626563/126‑1 (MQ=255)
                                                                        ttCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCAACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGCAAAAAAGATCTGGCCTTATTAGAAGTGGaa                                                         <  2:519850/133‑1 (MQ=255)
                                                                        ttCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCAACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGCAAAAAAGATCTGGCCTTATTAGAAGTGGaa                                                         <  2:499819/133‑1 (MQ=255)
                                                                           tctTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCAACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGCAAAAAAGATCTGGCCTTATTAGAAGTGGAAg                                                        >  2:459027/1‑131 (MQ=255)
                                                                           tctTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCAACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGCAAAAAAGATCTGGCCTTATTAGAAGTGGAAg                                                        <  1:459027/131‑1 (MQ=255)
                                                                            ctTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCAACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGCAAAAAAGATCTGGCCGTAGTAGaagcggaagcgg                                                     >  1:290041/1‑133 (MQ=255)
                                                                              tGGCATGTCACCTGGCTCATCTGTTAATGTGGTGCGCGTCGCTCCACTCGGCAACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGCAAAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGtt                                                   <  1:222970/133‑1 (MQ=255)
                                                                                gCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCAACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGCAAAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTc                                                 >  2:419578/1‑133 (MQ=255)
                                                                                 cATGTTACCTGGCTCCTCTTTTAATGTGCTGCTCGTCGCTCAACTCGGCAACCCCACTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGCAAAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTcc                                                <  2:185338/133‑1 (MQ=255)
                                                                                    gTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCAACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGCAAAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGt                                             <  2:69690/133‑1 (MQ=255)
                                                                                    gTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCAACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGCAAAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGt                                             <  1:602494/133‑1 (MQ=255)
                                                                                           ggCTCCTCTTTTAATTTGGTGCGCATCGCTCCCCTCGGCAACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGCAAAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACg                                      <  2:246860/133‑1 (MQ=255)
                                                                                           ggCTCCTCTTGTAATGTGGTGCGCCTCGTTCCACTCGGCAAACCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGCAAAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACg                                      <  2:229027/133‑1 (MQ=255)
                                                                                               cctcTTTTAATGTGGTGCGCGTCGCTCCACTCGGCAACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGCAAAAAAGATCTGGCCTTATTCGAAATTGAAGCGGTTTGCTGTAAATACGGTGa                                  >  2:280382/1‑133 (MQ=255)
                                                                                               cctcTTTTAATGTGGTGCGCGTCGCTCCACTCGGCAACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGCAAAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATa                                        <  2:413945/127‑1 (MQ=255)
                                                                                               cctcTTTTAATGTGGTGCGCGTCGCTCCACTCGGCAACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGCAAAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATa                                        >  1:413945/1‑127 (MQ=255)
                                                                                                ctctTTTAATGTGGTGCGCGTCGCTCCACTCGGCAACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGCAAAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGAt                                 >  1:10578/1‑133 (MQ=255)
                                                                                                ctctTTTAATGTGGTGCGCGTCGCTCCACTCGGCAACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGCAAAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGAt                                 >  2:98209/1‑133 (MQ=255)
                                                                                                ctctTTTAATGTGGTGCGCGTCGCTCCACTCGGCAACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGCAAAAAAGATCTGGCCTTATTAAAAGGGGAGGCGGTTTCGTGTTGATACGGTGAt                                 >  2:164958/1‑133 (MQ=255)
                                                                                                ccctTTTAAAGTGGTGCGCGTCGCCCCACTCGGCAACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGCAAAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGAt                                 <  2:290041/131‑1 (MQ=255)
                                                                                                  ctTTTAATGTGGTGCGCGTCGCTCCACTCGGCAACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGCAAAAAAGATCTGGCCTTATTAGAAGGGGAAGCGATTTTCTTTTAATAGGGTGATaa                               >  1:297055/1‑133 (MQ=255)
                                                                                                           tggtgCGCGTCGCTCGACGCGGTAACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGCAAAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGa                      <  1:204950/133‑1 (MQ=255)
                                                                                                           tggtgCGCGTCGCTCCACTCGGCAACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGCAAAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGa                      <  1:587993/133‑1 (MQ=255)
                                                                                                                        tCCACTCGGCAACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGCAAAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGAAAAAATTAACCAt         >  2:409739/1‑133 (MQ=255)
                                                                                                                              cGGCAACCTCATTCACATCGAACCCCGTCGCGTGAGCCTGGTATTACGCAAAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGAAAAAATTAACCATTGGCtt   <  1:164958/133‑1 (MQ=255)
                                                                                                                               ggCAACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGCAAAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGAAAAAATTAACCATTGGCTTa  >  1:450751/1‑133 (MQ=255)
                                                                                                                                  |                                                                                                                                 
CACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCTGGTATTACGCAAAAAAGATCTGGCCTTATTAGAAGTGGAAGCGGTTTCCTGTTAATACGGTGATAACAACAATGAAAAAATTAACCATTGGCTTA  >  NC_000913/3540171‑3540430

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