Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,297,451 T→A 6.2% 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 6.2% 236.2 / 5.7 65Q65L (CAG→CTG) ccmAheme export ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (32/29);  new base A (2/2);  total (34/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.52e-01

CCAGACATTGTGCGGTATCGCCATCGCGATGATAAAAGTGCAGATTTTCTAACGCCGTCAGCCGGGTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCACTCTCCTGCGTTCAGCGTAAATGACAAGCCACTAAATAA  >  NC_000913/2297314‑2297593
                                                                                                                                         |                                                                                                                                              
ccAGACATTGTGCGGTATCGCCATCGCGATGATAAAAGTGCAGATTTTCTAACGCCGTCAGCCGGGTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCa                                                                                                                                      >  1:605696/1‑148 (MQ=255)
     cATTGTGCGGTATCGCCATCGCGATGATAAAAGTGCAGATTTTCTAACGCCGTCAGCCGGGTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCAAGAGaa                                                                                                                                 <  2:605696/148‑1 (MQ=255)
                                    aGTGCAGATTTTCTAACGCCGTCAGCCGGGTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGaa                                                                                                                                 <  2:374269/117‑1 (MQ=255)
                                    aGTGCAGATTTTCTAACGCCGTCAGCCGGGTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGaa                                                                                                                                 >  1:374269/1‑117 (MQ=255)
                                                      ccGTCAGCCGGGTTTTGATGCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAAc                                                                                                       >  1:893662/1‑125 (MQ=255)
                                                      ccGTCAGCCGGGTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACaaa                                                                                                       <  2:893662/125‑2 (MQ=255)
                                                         tCAGCCGGGTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCAGCCCTTGcat                                                                                                                                      <  1:37789/91‑3 (MQ=255)
                                                         tCAGCCGGGTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCAGCCCTTGcat                                                                                                                                      >  2:37789/1‑89 (MQ=255)
                                                            gCCGGGTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCAGCCCTTGc                                                                                                                                        >  1:174331/1‑86 (MQ=255)
                                                            gCCGGGTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCAGCCCTTGc                                                                                                                                        <  2:174331/86‑1 (MQ=255)
                                                             ccGGGTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCt                                                                                                                                           >  1:226896/1‑82 (MQ=255)
                                                             ccGGGTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCt                                                                                                                                           <  2:226896/82‑1 (MQ=255)
                                                              cGGGTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCag                                                                                                                                     <  4:641132/87‑1 (MQ=255)
                                                              cGGGTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCag                                                                                                                                     >  3:641132/1‑87 (MQ=255)
                                                               gggTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTc                                                                                                                             <  4:454648/94‑1 (MQ=255)
                                                               gggTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTc                                                                                                                             >  3:454648/1‑94 (MQ=255)
                                                                       aTCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAAccc                                                                                                     >  3:306251/1‑110 (MQ=255)
                                                                       aTCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAAccc                                                                                                     <  4:306251/110‑1 (MQ=255)
                                                                         ccccGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGt                                                                                >  4:515588/1‑129 (MQ=255)
                                                                         ccccGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGt                                                                                <  3:515588/129‑1 (MQ=255)
                                                                            cGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAg                                                                                                                  >  4:281526/1‑92 (MQ=255)
                                                                            cGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAg                                                                                                                  <  3:281526/92‑1 (MQ=255)
                                                                                 gATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGt                                                     >  4:213560/1‑148 (MQ=255)
                                                                                     gCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGAtt                                                 <  3:213560/148‑1 (MQ=255)
                                                                                        tATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTg                                                                                                                        <  4:153449/74‑1 (MQ=255)
                                                                                        tATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTg                                                                                                                        >  3:153449/1‑74 (MQ=255)
                                                                                        tATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACg                                                                                         >  3:570532/1‑105 (MQ=255)
                                                                                        tATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACg                                                                                         <  4:570532/105‑1 (MQ=255)
                                                                                        tATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGt                                                                             >  4:704850/1‑117 (MQ=255)
                                                                                        tATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGt                                                                             <  3:704850/117‑1 (MQ=255)
                                                                                        tATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGt                                              >  3:734714/1‑148 (MQ=255)
                                                                                        tATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGt                                              <  4:734714/148‑1 (MQ=255)
                                                                                           ccATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACa                                                                                                         >  3:1087070/1‑86 (MQ=255)
                                                                                           ccATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACa                                                                                                         <  4:1087070/86‑1 (MQ=255)
                                                                                           ccATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGgaga                                                                                     <  2:61514/106‑1 (MQ=255)
                                                                                           ccATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGgaga                                                                                     >  1:61514/1‑106 (MQ=255)
                                                                                                     ttttGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACgg                                                                                        >  3:356562/1‑93 (MQ=255)
                                                                                                     ttttGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACgg                                                                                        <  4:356562/93‑1 (MQ=255)
                                                                                                          aTGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGc                                                            <  2:898531/116‑1 (MQ=255)
                                                                                                          aTGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGc                                                            >  1:898531/1‑116 (MQ=255)
                                                                                                           tGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCACTCTCCTGCGTTCa                           >  3:99071/1‑148 (MQ=255)
                                                                                                            ggTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCACTCTCCTGCGTTCAg                          >  3:722959/1‑148 (MQ=255)
                                                                                                                 cTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCACTCTCCTGCGTTCAGCGTaa                     <  4:722959/148‑1 (MQ=255)
                                                                                                                    tCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCACTCTcc                                   <  1:1066093/131‑1 (MQ=255)
                                                                                                                    tCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCACTCTcc                                   >  2:1066093/1‑131 (MQ=255)
                                                                                                                         tACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCActct                                     <  3:854061/124‑1 (MQ=255)
                                                                                                                         tACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCActct                                     >  4:854061/1‑124 (MQ=255)
                                                                                                                            cTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCaa                                                                                            >  3:451455/1‑66 (MQ=255)
                                                                                                                            cTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCaa                                                                                            <  4:451455/66‑1 (MQ=255)
                                                                                                                            cTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCACTCTCCTGCGTTCAGCGTAAATGACAAGCCa          >  2:1061781/1‑148 (MQ=255)
                                                                                                                             tGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCACTCTCCTGCGTTCAGCGTaaa                    >  1:20777/1‑137 (MQ=255)
                                                                                                                             tGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCACTCTCCTGCGTTCAGCGTaaa                    <  2:20777/137‑1 (MQ=255)
                                                                                                                             tGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCACTCTCCTGCGTTCAGCGTAAATg                  >  4:352481/1‑139 (MQ=255)
                                                                                                                             tGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCACTCTCCTGCGTTCAGCGTAAATg                  <  3:352481/139‑1 (MQ=255)
                                                                                                                             tGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCACTCTCCTGCGTTCAGCGTAAATGACAAGCCAc         <  4:99071/148‑1 (MQ=255)
                                                                                                                               aTGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCaaa                                                                                           >  1:754819/1‑64 (MQ=255)
                                                                                                                               aTGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCaaa                                                                                           <  2:754819/64‑1 (MQ=255)
                                                                                                                                tGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCACTCTCCTGCGTTCAg                          <  2:308826/128‑1 (MQ=255)
                                                                                                                                tGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCACTCTCCTGCGTTCAg                          >  1:308826/1‑128 (MQ=255)
                                                                                                                                  cAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCACTCTCCTGCGTTCAGCGTAAATGACAAGCCACTAAAt    >  1:375090/1‑148 (MQ=255)
                                                                                                                                   aaGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCACTCTCCTGCGTTCAGCGTAAATGACAAGCCACTAAATa   <  1:1061781/148‑1 (MQ=255)
                                                                                                                                   aaGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCACTCTCCTGCGTTCAGCGTAAATGACAAGCCACTAAATa   >  3:1042470/1‑148 (MQ=255)
                                                                                                                                    aGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCACTCTCCTGCGTTCAGCGTAAATGACAAGCCACTAAATaa  >  1:477139/1‑148 (MQ=255)
                                                                                                                                     gggCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCACtctc                                    >  4:335609/1‑113 (MQ=255)
                                                                                                                                     gggCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCACtctc                                    <  3:335609/113‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                              
CCAGACATTGTGCGGTATCGCCATCGCGATGATAAAAGTGCAGATTTTCTAACGCCGTCAGCCGGGTTTTGATCCCCGGCTGATGGCCTATCCATAACAGGTTTTGATGGTAGCTGTCGCGTACCTGATGCAAGGGCTGCCCTTGCCAGAGAACCTCGCCTGCGTCAGGGCGAGACAACCCCGTCAGCAAACGGAGAAGCGTTGTCTTCCCCGCGCCGTTGCTACCGGTGATTTGTACCCACTCTCCTGCGTTCAGCGTAAATGACAAGCCACTAAATAA  >  NC_000913/2297314‑2297593

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