Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 515,808 C→T intergenic (‑35/‑111) qmcA ← / → fetA PHB domain membrane‑anchored putative protease/iron export ABC transporter ATPase, peroxide resistance protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913515,8080CT100.0% 126.0 / NA 39intergenic (‑35/‑111)qmcA/fetAPHB domain membrane‑anchored putative protease/iron export ABC transporter ATPase, peroxide resistance protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/18);  total (21/18)

TGATAGCCCTGCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGCAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCTT  >  NC_000913/515691‑515937
                                                                                                                     |                                                                                                                                 
tGATAGCCCTGCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTaaa                                                                                                               <  2:331291/138‑1 (MQ=255)
tGATAGCCCTGCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTaaa                                                                                                               >  1:331291/1‑138 (MQ=255)
  aTAGCCCTGCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACaa                                                                                                            >  2:421441/1‑139 (MQ=255)
     gCCCTGCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAAtt                                                                                                         >  1:53476/1‑139 (MQ=255)
              tACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGAt                                                                                                <  2:53476/139‑1 (MQ=255)
                cGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTa                                                                                              <  2:967607/139‑1 (MQ=255)
                       gACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAag                                                                                       <  1:421441/139‑1 (MQ=255)
                               cgcCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTcc                                                                               >  2:710208/1‑139 (MQ=255)
                                      aaTGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTAt                                                                        <  1:710208/139‑1 (MQ=255)
                                          aCCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATtgctg                                                                                                                         >  1:610044/1‑86 (MQ=255)
                                          aCCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATtgctg                                                                                                                         <  2:610044/86‑1 (MQ=255)
                                            cAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCCTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATTCCTGc                                                                  <  2:720077/139‑1 (MQ=255)
                                              gcgcGACAAAAATGAGAATCGGGATACAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTaaa                                                                                                               <  2:773834/92‑1 (MQ=255)
                                              gcgcGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTaaa                                                                                                               >  1:773834/1‑92 (MQ=255)
                                               cgcgACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTaaa                                                                                                               >  1:879481/1‑91 (MQ=255)
                                               cgcgACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTaaa                                                                                                               <  2:879481/91‑1 (MQ=255)
                                               cgcgACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATTCCTGCGaa                                                               >  1:879908/1‑139 (MQ=255)
                                                                              aGCATCGGAAAAACCTCCTGTTGTACCGTTCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGaaaa                                >  2:1021915/1‑139 (MQ=255)
                                                                              aGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTgctgct                                                                                                                       >  1:315116/1‑52 (MQ=255)
                                                                              aGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTgctgct                                                                                                                       <  2:315116/52‑1 (MQ=255)
                                                                                    ggAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATTCCTGCGAAGGATaaa                                                        >  1:681670/1‑109 (MQ=255)
                                                                                    ggAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATTCCTGCGAAGGATaaa                                                        <  2:681670/109‑1 (MQ=255)
                                                                                    ggAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGaa                          >  2:24586/1‑139 (MQ=255)
                                                                                             tcctGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTc                                                                                >  1:1120692/1‑76 (MQ=255)
                                                                                             tcctGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTc                                                                                <  2:1120692/76‑1 (MQ=255)
                                                                                              cctGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATTCCTGCGAAgg                                                             <  2:465310/94‑1 (MQ=255)
                                                                                              cctGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATTCCTGCGAAgg                                                             >  1:465310/1‑94 (MQ=255)
                                                                                                    gtACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATTCCTGCGAAGGATaaa                                                        <  2:309989/93‑1 (MQ=255)
                                                                                                    gtACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATTCCTGCGAAGGATaaa                                                        >  1:309989/1‑93 (MQ=255)
                                                                                                         gTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTccccc                                                                            >  2:1046134/1‑68 (MQ=255)
                                                                                                         gTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTccccc                                                                            <  1:1046134/68‑1 (MQ=255)
                                                                                                         gTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTc     >  2:532178/1‑139 (MQ=255)
                                                                                                          tCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATg                  >  2:544565/1‑125 (MQ=255)
                                                                                                          tCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATg                  <  1:544565/125‑1 (MQ=255)
                                                                                                            cATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATAc                                                                                                     <  1:747778/40‑1 (MQ=255)
                                                                                                            cATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATAc                                                                                                     >  2:747778/1‑40 (MQ=255)
                                                                                                            cATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAag                                                                                       <  2:562175/54‑1 (MQ=255)
                                                                                                            cATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAag                                                                                       >  1:562175/1‑54 (MQ=255)
                                                                                                            cATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCtt  >  1:870598/1‑139 (MQ=255)
                                                                                                                     |                                                                                                                                 
TGATAGCCCTGCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGCAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCTT  >  NC_000913/515691‑515937

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