Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,297,451 T→A 7.4% Q65L (CAG→CTG)  ccmA ← heme export ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,297,4510TA 7.4% 177.9 / 4.6 54Q65L (CAG→CTG) ccmAheme export ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (26/24);  new base A (2/2);  total (28/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

CGCCAGGGCTTCCAGACATTGTGCGGTATCGCCATCGCGATGATAAAAGTGCAGATTTTCTAACGCCGTCAGCCGGGTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCACTCTCCTGCGTTCAGCGTAAATGACAAGCCACTAAATAAGG  >  NC_000913/2297303‑2297595
                                                                                                                                                    |                                                                                                                                                
cgcCAGGGCTTCCAGACATTGTGCGGTATCGCCATCGCGATGATAAAAGTGCAGATTTTCTAACGCCGTCAGCCGGGTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCt                                                                                                                                                  <  1:166392/149‑1 (MQ=255)
  ccAGGGCTTCCAGACATTGTGCGGTATCGCCATCGCGATGATAAAAGTGCAGATTTTCTAACGCCGTCAGCCGGGTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGc                                                                                                                                                >  2:356798/1‑149 (MQ=255)
         ttCCAGACATTGTGCGGTATCGCCATCGCGATGATAAAAGTGCAGATTTTCTAACGCCGTCAGCCGGGTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGcc                                                                                                                                         >  4:490362/1‑149 (MQ=255)
         ttCCAGACATTGTGCGGTATCGCCATCGCGATGATAAAAGTGCAGATTTTCTAACGCCGTCAGCCGGGTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGcc                                                                                                                                         >  2:3397/1‑149 (MQ=255)
          tCCAGACATTGTGCGGTATCGCCATCGCGATGATAAAAGTGCAGATTTTCTAACGCCGTCAGCCGGGTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCAGCCCTTGcat                                                                                                                                        <  4:556210/149‑3 (MQ=255)
            cAGACATTGTGCGGTATCGCCATCGCGATGATAAAAGTGCAGATTTTCTAACGCCGTCAGCCGGGTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCaga                                                                                                                                      >  3:465513/1‑149 (MQ=255)
              gACATTGTGCGGTATCGCCATCGCGATGATAAAAGTGCAGATTTTCTAACGCCGTCAGCCGGGTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAgaga                                                                                                                                    >  3:308843/1‑149 (MQ=255)
              gACATTGTGCGGTATCGCCATCGCGATGATAAAAGTGCAGATTTTCTAACGCCGTCAGCCGGGTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAgaga                                                                                                                                    >  3:244392/1‑149 (MQ=255)
               aCATTGTGCGGTATCGCCATCGCGATGATAAAAGTGCAGATTTTCTAACGCCGTCAGCCGGGTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGaa                                                                                                                                   >  3:295350/1‑149 (MQ=255)
                cATTGTGCGGTATCGCCATCGCGATGATAAAAGTGCAGATTTTCTAACGCCGTCAGCCGGGTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAAc                                                                                                                                  <  4:580688/149‑1 (MQ=255)
                         gTATCGCCATCGCGATGATAAAAGTGCAGATTTTCTAACGCCGTCAGCCGGGTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGc                                                                                                                         <  4:583424/149‑1 (MQ=255)
                         gTATCGCCATCGCGATGATAAAAGTGCAGATTTTCTAACGCCGTCAGCCGGGTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGc                                                                                                                         >  4:402491/1‑149 (MQ=255)
                               ccATCGCGATGATAAAAGTGCAGATTTTCTAACGCCGTCAGCCGGGTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAgg                                                                                                                   >  1:116066/1‑149 (MQ=255)
                                        tgatAAAAGTGCAGATTTTCTAACGCCGTCAGCCGGGTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACaa                                                                                                          <  3:312603/149‑1 (MQ=255)
                                           tAAAAGTGCAGATTTTCTAACGCCGTCAGCCGGGTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAAccc                                                                                                       <  2:116066/149‑1 (MQ=255)
                                            aaaaGTGCAGATTTTCTAACGCCGTCAGCCGGGTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAAcccc                                                                                                      <  4:308843/149‑1 (MQ=255)
                                                gTGCAGATTTTCTAACGCCGTCAGCCGGGTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCa                                                                                                  <  1:191655/149‑1 (MQ=255)
                                                     gATTTTCTAACGCCGTCAGCCGGGTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCAGCCCTTGcatca                                                                                                                                      <  4:221045/108‑5 (MQ=255)
                                                     gATTTTCTAACGCCGTCAGCCGGGTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCAGCCCTTGcatca                                                                                                                                      >  3:221045/1‑104 (MQ=255)
                                                      aTTTTCTAACGCCGTCAGCCGGGTTTTGATCCCCGGCTGATGGGCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACATGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGAAAACCCCGTCAGCAAAc                                                                                            <  2:419650/149‑1 (MQ=255)
                                                      aTTTTCTAACGCCGTCAGCCGGGTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTCGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAAc                                                                                            <  4:465513/149‑1 (MQ=255)
                                                            tAACGCCGTCAGCCGGGTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGaa                                                                                      <  1:191125/149‑1 (MQ=255)
                                                            tAACGCCGTCAGCCGGGTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGaa                                                                                      <  3:490362/149‑1 (MQ=255)
                                                                gCCGTCAGCCGGGTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGt                                                                                  >  2:422769/1‑149 (MQ=255)
                                                                    tCAGCCGGGTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGTGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCCGCAAACGGAGAAGCGTTGTc                                                                              >  1:204891/1‑149 (MQ=255)
                                                                     cAGCCGGGTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCt                                                                             <  1:229396/149‑1 (MQ=255)
                                                                            gTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCcgc                                                                      <  3:28353/149‑1 (MQ=255)
                                                                             ttttGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCcgcg                                                                     >  2:248578/1‑149 (MQ=255)
                                                                                   tCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTg                                                               >  2:133341/1‑149 (MQ=255)
                                                                                   tCCCCGGCTGATGGCCTATCCATAACAGGTGTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGCCAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTg                                                               >  4:521436/1‑149 (MQ=255)
                                                                                   tCCCCGGCTGATGGCCTATCCATAACAGGTATTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCCGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTg                                                               >  2:252234/1‑149 (MQ=255)
                                                                                           tGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGt                                                       >  2:180825/1‑149 (MQ=255)
                                                                                             aTGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGa                                                     <  2:108574/149‑1 (MQ=255)
                                                                                              tGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGAt                                                    <  1:3397/149‑1 (MQ=255)
                                                                                                  cTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGt                                                >  2:519477/1‑149 (MQ=255)
                                                                                                  cTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGt                                                >  1:21344/1‑149 (MQ=255)
                                                                                                   tATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGtt                                                                <  4:310867/132‑1 (MQ=255)
                                                                                                   tATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGtt                                                                >  3:310867/1‑132 (MQ=255)
                                                                                                    aTCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTAc                                              <  2:340155/149‑1 (MQ=255)
                                                                                                       cATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCa                                           <  1:248578/149‑1 (MQ=255)
                                                                                                        aTAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCAc                                          <  1:252234/149‑1 (MQ=255)
                                                                                                         tAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCAct                                         >  4:280712/1‑149 (MQ=255)
                                                                                                           aCAGGTTTAGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCActct                                       >  3:539440/1‑149 (MQ=255)
                                                                                                                 tttGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCACTCTCCTGCg                                 <  2:521210/149‑1 (MQ=255)
                                                                                                                   tGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCACTCTCCTGCGtt                               <  1:180825/149‑1 (MQ=255)
                                                                                                                    gATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCACTCTCCTGCGTTc                              <  4:456240/149‑1 (MQ=255)
                                                                                                                                       cTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCACTCTCCTGCGTTCAGCGTAAATGACAAGCCAc           >  3:273042/1‑149 (MQ=255)
                                                                                                                                            gCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCACTCTCCTGCGTTCAGCGTAAATGACAAGCCACTAAAt      >  2:427138/1‑149 (MQ=255)
                                                                                                                                            gCAAGGGCAGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCACTCTCCTGCGTTCAGCGTAAATGACAAGCCACTAAAt      >  1:510295/1‑149 (MQ=255)
                                                                                                                                             cAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCACTCTCCTGCGTTCAGCGTAAATGACAAGCCACTAAATa     >  3:379770/1‑149 (MQ=255)
                                                                                                                                             cAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCACTCTCCTGCGTTCAGCGTAAATGACAAGCCACTAAATa     >  3:252311/1‑149 (MQ=255)
                                                                                                                                              aaGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCACTCTCCTGCGTTCAGCGTAAATGACAAGCCa            <  1:222934/141‑1 (MQ=255)
                                                                                                                                              aaGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCACTCTCCTGCGTTCAGCGTAAATGACAAGCCa            >  2:222934/1‑141 (MQ=255)
                                                                                                                                               aGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCACTCTCCTGCGTTCAGCGTAAATGACAAGCCACTAAATAAg   <  4:223140/149‑1 (MQ=255)
                                                                                                                                                gggCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCACTCTCCTGCGTTCAGCGTAAATGACAAGCCACTAAATAAgg  <  4:539440/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
CGCCAGGGCTTCCAGACATTGTGCGGTATCGCCATCGCGATGATAAAAGTGCAGATTTTCTAACGCCGTCAGCCGGGTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCACTCTCCTGCGTTCAGCGTAAATGACAAGCCACTAAATAAGG  >  NC_000913/2297303‑2297595

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