Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,297,451 T→A 9.8% 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 9.8% 162.9 / 9.1 51Q65L (CAG→CTG) ccmAheme export ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (33/13);  new base A (3/2);  total (36/15)
Fisher's exact test for biased strand distribution p-value = 6.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.41e-01

GGCTTCCAGACATTGTGCGGTATCGCCATCGCGATGATAAAAGTGCAGATTTTCTAACGCCGTCAGCCGGGTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCACTCTCCTGCGTTCAGCGTAAATGACAAGCCACTAAATAAGGTTCG  >  NC_000913/2297309‑2297599
                                                                                                                                              |                                                                                                                                                    
ggCTTCCAGACATTGTGCGGTATCGCCATCGCGATGATAAAAGTGCAGATTTTCTAACGCCGTCAGCCGGGTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCtt                                                                                                                                                >  4:88516/1‑149 (MQ=255)
      cAGACATTGTGCGGTATCGCCATCGCGATGATAAAAGTGCAGATTTTCTAACGCCGTCAGCCGGGTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCAGCCCTTGcat                                                                                                                                            <  4:340327/147‑3 (MQ=255)
      cAGACATTGTGCGGTATCGCCATCGCGATGATAAAAGTGCAGATTTTCTAACGCCGTCAGCCGGGTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCAGCCCTTGcat                                                                                                                                            >  3:340327/1‑145 (MQ=255)
           aTTGTGCGGTATCGCCATCGCGATGATAAAAGTGCAGATTTTCTAACGCCGTCAGCCGGGTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAAcc                                                                                                                                     >  1:250970/1‑149 (MQ=255)
                          cATCGCGATGATAAAAGTGCAGATTTTCTAACGCCGTCAGCCGGGTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAAc                                                                                                                                      <  1:450843/133‑1 (MQ=255)
                          cATCGCGATGATAAAAGTGCAGATTTTCTAACGCCGTCAGCCGGGTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAAc                                                                                                                                      >  2:450843/1‑133 (MQ=255)
                            tCGCGATGATAAAAGTGCAGATTTTCTAACGCCGTCAGCCGGGTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCg                                                                                                                    >  1:261927/1‑149 (MQ=255)
                                                  tttCTAACGCCGTCAGCCGGGTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACgg                                                                                              <  1:577940/149‑1 (MQ=255)
                                                    tCTAACGCCGTCAGCCGGGTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGgag                                                                                            <  1:28454/149‑1 (MQ=255)
                                                      tAACGCCGTCAGCCGGGTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGaa                                                                                          >  3:371023/1‑149 (MQ=255)
                                                      tAACGCCGTCAGCCGGGTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGaa                                                                                          <  3:213545/149‑1 (MQ=255)
                                                        aCGCCGTCAGCCGGGTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGTGAGACAACCCCGTCAGCAAACGGAGAAGc                                                                                        <  2:261927/149‑1 (MQ=255)
                                                           ccGTCAGCCGGGTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGtt                                                                                     <  4:347385/149‑1 (MQ=255)
                                                             gTCAGCCGGGTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGt                                                                                   >  3:256240/1‑149 (MQ=255)
                                                                    gggTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTcccc                                                                            >  4:359727/1‑149 (MQ=255)
                                                                           gATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGt                                                                     >  3:115984/1‑149 (MQ=255)
                                                                            aTCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAAc                                                                                                <  1:396271/121‑1 (MQ=255)
                                                                            aTCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAAc                                                                                                >  2:396271/1‑121 (MQ=255)
                                                                            aTCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGtt                                                                    >  4:501214/1‑149 (MQ=255)
                                                                            aTCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGtt                                                                    >  3:577393/1‑149 (MQ=255)
                                                                               cccGGCTGATGGCCTATCCATAACGGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCt                                                                 >  1:530170/1‑149 (MQ=255)
                                                                               cccGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCt                                                                 >  3:403639/1‑149 (MQ=255)
                                                                                ccGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTa                                                                >  3:186986/1‑149 (MQ=255)
                                                                                ccGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTa                                                                >  1:449778/1‑149 (MQ=255)
                                                                                 cGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTAc                                                               >  4:322376/1‑149 (MQ=255)
                                                                                     tGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGt                                                           >  3:301271/1‑149 (MQ=255)
                                                                                          gCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGAttt                                                      >  4:464052/1‑149 (MQ=255)
                                                                                            cTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGt                                                    >  4:60008/1‑149 (MQ=255)
                                                                                            cTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGt                                                    >  2:370780/1‑149 (MQ=255)
                                                                                                 cATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCa                                               >  3:58063/1‑149 (MQ=255)
                                                                                                   tAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCAct                                             >  1:544892/1‑149 (MQ=255)
                                                                                                       aGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCACTCTcc                                         <  1:276512/149‑1 (MQ=255)
                                                                                                          ttttGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCACTCTCCTGc                                      >  2:350165/1‑149 (MQ=255)
                                                                                                                 ggTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGa                                                         >  3:470248/1‑123 (MQ=255)
                                                                                                                 ggTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGa                                                         <  4:470248/123‑1 (MQ=255)
                                                                                                                 ggTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCACTCTCCTGCGTTCAGc                               >  2:433526/1‑149 (MQ=255)
                                                                                                                  gTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCACTCTCCTGCGTTCAGCg                              >  2:86467/1‑149 (MQ=255)
                                                                                                                     gCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCACTCTCCTGCGTTCAGCGTaa                           <  4:261447/149‑1 (MQ=255)
                                                                                                                      cTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCACTCTCCTGCGTTCAGCGTaaa                          <  4:58063/149‑1 (MQ=255)
                                                                                                                           gcgTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCACTCTCCTGCGTTCAGCGTAAATGACa                     >  3:518523/1‑149 (MQ=255)
                                                                                                                              tACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCACTCTCCTGCGTTCAGCGTAAATGACAAGc                  <  4:115984/149‑1 (MQ=255)
                                                                                                                                 cTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCACTCTCCTGCGTTCAGCGTAAATGACAAGCCAc               >  3:28307/1‑149 (MQ=255)
                                                                                                                                 cTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCACTCTCCTGCGTTCAGCGTAAATGACAAGCCAc               >  4:549202/1‑149 (MQ=255)
                                                                                                                                  tGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCACTCTCCTGCGTTCAGCGTAAATGACAAGCCACt              >  4:95507/1‑149 (MQ=255)
                                                                                                                                    tggCAAGGGCAGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGAtt                                                       >  1:196901/3‑106 (MQ=255)
                                                                                                                                    tggCAAGGGCAGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGAtt                                                       <  2:196901/104‑1 (MQ=255)
                                                                                                                                      gCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCACTCTCCTGCGTTCAGCGTAAATGACAAGCCACTAAAt          >  1:396559/1‑149 (MQ=255)
                                                                                                                                        aaGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCACTCTCCTGCGTTCAGCGTAAATGACAAGCCACTAAATaa        <  4:186986/149‑1 (MQ=255)
                                                                                                                                        aaGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCACTCTCCTGCGTTCAGCGTAAATGACAAGCCACTAAATaa        >  4:140378/1‑149 (MQ=255)
                                                                                                                                        aaGGGCAGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCACTCTCCTGCGTTCAGCGTAAATGACAAGCCACTAAATaa        >  1:558968/1‑149 (MQ=255)
                                                                                                                                            gCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCACTCTCCTGCGTTCAGCGTAAATGACAAGCCACTAAATAAGgtt    >  2:104814/1‑149 (MQ=255)
                                                                                                                                              tGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCACTCTCCTGCGTTCAGCGTAAATGACAAGCCACTAAATAAGgttcg  >  2:299496/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
GGCTTCCAGACATTGTGCGGTATCGCCATCGCGATGATAAAAGTGCAGATTTTCTAACGCCGTCAGCCGGGTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCACTCTCCTGCGTTCAGCGTAAATGACAAGCCACTAAATAAGGTTCG  >  NC_000913/2297309‑2297599

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