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_000913	515,808	0	C	T	100.0%	163.5 / NA	51	intergenic (‑35/‑111)	qmcA/fetA	PHB 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 (24/27);  total (24/27)

CGTTCTACTGTCCACTGATAGCCCTGCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGCAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATT‑CCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCTT  >  NC_000913/515676‑515937                                                                                                                                     |                                                                                                                                   cGTTCTACTGTCCACTGATAGCCCTGCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAAtt                                                                                                                              <  1:918666/139‑1 (MQ=255)   ttCTACTGTCCACTGATAGCCCTGCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATtgc                                                                                                                            >  2:50645/1‑139 (MQ=255)     cTACTGTCCACTGATAGCCCTGCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATtgctg                                                                                                                          >  2:963388/1‑139 (MQ=255)     cTACTGTCCACTGATAGCCCTGCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATtgctg                                                                                                                          <  1:963388/139‑1 (MQ=255)           tCCACTGATAGCCCTGCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGAt                                                                                                                    >  1:118501/1‑139 (MQ=255)               cTGATAGCCCTGCGGTACGATTTTGGCACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTaaa                                                                                                                <  2:118501/139‑1 (MQ=255)               cTGATAGCCCTGCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTaaa                                                                                                                <  2:50560/139‑1 (MQ=255)               cTGATAGCCCTGCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTaaa                                                                                                                <  2:818095/139‑1 (MQ=255)               cTGATAGCCCTGCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTaaa                                                                                                                <  2:534740/139‑1 (MQ=255)                           cGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATtgctg                                                                                                                          <  1:610525/117‑1 (MQ=255)                           cGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATtgctg                                                                                                                          >  2:610525/1‑117 (MQ=255)                              tACGATTTTGGCACCCGCGTCCAAAAAGACCAGCGCGACAAAAATGAGAATCGGGATAAAGGTAAGCATCGGAAAATCCCCCTGTTGTAACGTCCATAATCCATAAAATTGCTTCTTTATTAAACAAAATATACCTGAt                                                                                                 <  2:859902/139‑1 (MQ=255)                              tACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGAt                                                                                                 <  2:434886/139‑1 (MQ=255)                              tACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGAt                                                                                                 <  2:278972/139‑1 (MQ=255)                                cGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTaaa                                                                                                                <  1:405607/122‑1 (MQ=255)                                cGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTaaa                                                                                                                >  2:405607/1‑122 (MQ=255)                                        acacCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAaga                                                                                       >  1:983735/1‑139 (MQ=255)                                             cgcgCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGtt                                                                                  >  1:62714/1‑139 (MQ=255)                                             cgcgCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGtt                                                                                  >  2:913936/1‑139 (MQ=255)                                             cgcgCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGt                                                                                   <  1:270675/138‑1 (MQ=255)                                             cgcgCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGt                                                                                   >  2:270675/1‑138 (MQ=255)                                             cgcgCCGACAATGACCAGCGCAACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGtt                                                                                  <  1:1312358/139‑1 (MQ=255)                                              gcgcCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTc                                                                                 <  1:474516/139‑1 (MQ=255)                                                gcCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTccc                                                                               <  1:913936/139‑1 (MQ=255)                                                   gACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCt                                                                            <  2:983735/139‑1 (MQ=255)                                                        tGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATT‑c                                                                      >  1:904458/1‑139 (MQ=255)                                                              gcgcGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATT‑CCTGCGa                                                                <  2:1168971/139‑1 (MQ=255)                                                               cgcgACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATT‑CCTGCGaa                                                               >  1:1209243/1‑139 (MQ=255)                                                                   aCAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATT‑CCTGCGAAGGAt                                                           <  2:350205/139‑1 (MQ=255)                                                                       aaaTGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATT‑CCTGCGAAGGATAAAc                                                       <  2:1209243/139‑1 (MQ=255)                                                                                           aTAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAag                                                                                        >  1:1008339/1‑87 (MQ=255)                                                                                           aTAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAag                                                                                        <  2:1008339/87‑1 (MQ=255)                                                                                                 aTCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATT‑CCTGc                                                                  >  2:1094617/1‑102 (MQ=255)                                                                                                 aTCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATT‑CCTGc                                                                  <  1:1094617/102‑1 (MQ=255)                                                                                                   cGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATT‑CCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGa                           >  1:530664/1‑139 (MQ=255)                                                                                                     gAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATtata                                                                                                       <  2:397148/62‑1 (MQ=255)                                                                                                     gAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATtata                                                                                                       >  1:397148/1‑62 (MQ=255)                                                                                                     gAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGAt                                                                                                 <  1:943579/68‑1 (MQ=255)                                                                                                     gAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGAt                                                                                                 >  2:943579/1‑68 (MQ=255)                                                                                                     gAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATT‑CCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAGGGGAAc                         >  2:159998/1‑139 (MQ=255)                                                                                                       aaaaCCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATT‑CCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGAAACAg                       <  1:61122/139‑1 (MQ=255)                                                                                                                tgttgtACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATTCCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAg               <  1:448204/139‑1 (MQ=255)                                                                                                                tgttgtACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATT‑CCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAgg              >  2:336761/1‑139 (MQ=255)                                                                                                                         gTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATT‑CCTGCGAAgg                                                             <  1:573132/83‑1 (MQ=255)                                                                                                                         gTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATT‑CCTGCGAAgg                                                             >  2:573132/1‑83 (MQ=255)                                                                                                                         gTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATT‑CCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGaaaaa                               >  2:1011581/1‑113 (MQ=255)                                                                                                                         gTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATT‑CCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGaaaaa                               <  1:1011581/113‑1 (MQ=255)                                                                                                                         gTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATT‑CCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTc     >  2:873512/1‑139 (MQ=255)                                                                                                                          tCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATT‑CCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTcc    >  2:415408/1‑139 (MQ=255)                                                                                                                            cATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATT‑CCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCtt  >  1:13179/1‑139 (MQ=255)                                                                                                                            cATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCGTATT‑CCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCtt  >  2:895111/1‑139 (MQ=255)                                                                                                                                     |                                                                                                                                   CGTTCTACTGTCCACTGATAGCCCTGCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGCAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATT‑CCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCTT  >  NC_000913/515676‑515937

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: ‑