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% 85.0 / NA 27intergenic (‑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 (13/14);  total (13/14)

CTACTGTCCACTGATAGCCCTGCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGCAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCTT  >  NC_000913/515680‑515937
                                                                                                                                |                                                                                                                                 
cTACTGTCCACTGATAGCCCTGCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATtgctg                                                                                                                         <  1:737877/139‑1 (MQ=255)
        cACTGATAGCCCTGCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTa                                                                                                                 >  1:999428/1‑139 (MQ=255)
          cTGATAGCCCTGCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTaaa                                                                                                               <  2:496277/139‑1 (MQ=255)
          cTGATAGCCCTGCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTaaa                                                                                                               <  2:468780/139‑1 (MQ=255)
                    tGCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATAc                                                                                                     <  1:462416/139‑1 (MQ=255)
                       ggTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTg                                                                                                  >  1:270354/1‑139 (MQ=255)
                                          cgcCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTc                                                                                >  1:1102067/1‑138 (MQ=255)
                                          cgcCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTc                                                                                <  2:1102067/138‑1 (MQ=255)
                                                cAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTa                                                                         <  2:383186/139‑1 (MQ=255)
                                                          cgcgACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTg                                                                                                                     >  2:575287/1‑85 (MQ=255)
                                                          cgcgACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTg                                                                                                                     <  1:575287/85‑1 (MQ=255)
                                                                aaaaaTGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATTCCTGCGAAGGATaa                                                         <  2:270354/139‑1 (MQ=255)
                                                                                         aGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATTc                                                                      >  2:1066734/1‑101 (MQ=255)
                                                                                         aGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATTc                                                                      <  1:1066734/101‑1 (MQ=255)
                                                                                         aGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGaaaa                                >  2:1044718/1‑139 (MQ=255)
                                                                                         aGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGaaaa                                >  2:318511/1‑139 (MQ=255)
                                                                                         aGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGaaaa                                >  1:658452/1‑139 (MQ=255)
                                                                                                gAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATTCCTGCGAAgg                                                             <  1:538928/103‑1 (MQ=255)
                                                                                                gAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATTCCTGCGAAgg                                                             >  2:538928/1‑103 (MQ=255)
                                                                                                       ctcctGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATg                  <  2:673084/139‑1 (MQ=255)
                                                                                                                    gTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATTCCTGCg                                                                 >  2:357392/1‑79 (MQ=255)
                                                                                                                    gTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATTCCTGCg                                                                 <  1:357392/79‑1 (MQ=255)
                                                                                                                     tCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATTCCTGCg                                                                 <  2:260785/78‑1 (MQ=255)
                                                                                                                     tCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATTCCTGCg                                                                 >  1:687745/1‑78 (MQ=255)
                                                                                                                     tCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATTCCTGCg                                                                 >  1:260785/1‑78 (MQ=255)
                                                                                                                     tCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATTCCTGCg                                                                 <  2:687745/78‑1 (MQ=255)
                                                                                                                       cATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCtt  >  2:768982/1‑139 (MQ=255)
                                                                                                                                |                                                                                                                                 
CTACTGTCCACTGATAGCCCTGCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGCAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCTT  >  NC_000913/515680‑515937

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