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%	168.7 / NA	52	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 (27/25);  total (27/25)

CTACTGTCCACTGATAGCCCTGCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGCAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCTTTGCT  >  NC_000913/515680‑515941                                                                                                                                 |                                                                                                                                      cTACTGTCCACTGATAGCCCTGCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATtgctg                                                                                                                             <  2:277403/139‑1 (MQ=255)  tACTGTCCACTGATAGCCCTGCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATtgctgc                                                                                                                            >  1:588403/1‑139 (MQ=255)       tCCACTGATAGCCCTGCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGAt                                                                                                                       >  1:1715026/1‑139 (MQ=255)                      gCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACa                                                                                                                 >  2:744562/1‑130 (MQ=255)                      gCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACa                                                                                                                 <  1:744562/130‑1 (MQ=255)                         gTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGa                                                                                                     >  2:1422851/1‑139 (MQ=255)                             gATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTAc                                                                                                 <  2:586139/139‑1 (MQ=255)                                         cgcgCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTgctgct                                                                                                                           <  1:206609/101‑1 (MQ=255)                                         cgcgCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTgctgct                                                                                                                           >  2:206609/1‑101 (MQ=255)                                         cgcgCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGtt                                                                                     <  1:1422851/139‑1 (MQ=255)                                                        cAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATTCCTGc                                                                      <  1:1405930/139‑1 (MQ=255)                                                          gcgcGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGtt                                                                                     <  2:1713280/122‑1 (MQ=255)                                                          gcgcGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGtt                                                                                     >  1:1713280/1‑122 (MQ=255)                                                           cgcgACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTc                                                                                    >  1:145044/1‑122 (MQ=255)                                                           cgcgACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTc                                                                                    <  2:145044/122‑1 (MQ=255)                                                           cgcgACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATTCCTGCGaa                                                                   >  2:1663221/1‑139 (MQ=255)                                                            gcgACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTaa                                                                                                                    >  1:641752/1‑89 (MQ=255)                                                            gcgACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTaa                                                                                                                    <  2:641752/89‑1 (MQ=255)                                                                cAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATtgctg                                                                                                                             >  1:2524/1‑76 (MQ=255)                                                                cAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATtgctg                                                                                                                             <  2:2524/76‑1 (MQ=255)                                                                cAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTc                                                                                    <  2:203261/117‑1 (MQ=255)                                                                cAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTc                                                                                    >  1:203261/1‑117 (MQ=255)                                                                        agaATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATTCCTGCGAAGGATAAACTGttt                                                      >  1:492387/1‑139 (MQ=255)                                                                          aaTCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAAt                                                                                                              >  2:210777/1‑81 (MQ=255)                                                                          aaTCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAAt                                                                                                              <  1:210777/81‑1 (MQ=255)                                                                          aaTCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATta                                                                                                            <  1:483948/83‑1 (MQ=255)                                                                          aaTCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATta                                                                                                            >  2:483948/1‑83 (MQ=255)                                                                                        tAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTaa                                                                                                                    <  1:472575/61‑1 (MQ=255)                                                                                        tAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTaa                                                                                                                    >  2:472575/1‑61 (MQ=255)                                                                                               cGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTgctgct                                                                                                                           <  2:232031/47‑1 (MQ=255)                                                                                               cGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTgctgct                                                                                                                           >  1:232031/1‑47 (MQ=255)                                                                                               cGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTccccc                                                                                >  2:1145421/1‑90 (MQ=255)                                                                                               cGGAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTccccc                                                                                <  1:1145421/90‑1 (MQ=255)                                                                                                ggAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATTCCTGc                                                                      <  2:1512972/99‑1 (MQ=255)                                                                                                ggAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATTCCTGc                                                                      >  1:1512972/1‑99 (MQ=255)                                                                                                 gAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAagagag                                                                                       <  2:296377/81‑1 (MQ=255)                                                                                                 gAAAAACCTCCTGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAagagag                                                                                       >  1:296377/1‑81 (MQ=255)                                                                                                        ctcctGTTGTACCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATg                      <  2:1715026/139‑1 (MQ=255)                                                                                                                  aCCGTCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGaaatag            >  2:340715/1‑136 (MQ=255)                                                                                                                      tCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTAc                                                                                                 >  2:215355/1‑50 (MQ=255)                                                                                                                      tCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTAc                                                                                                 >  1:413305/1‑50 (MQ=255)                                                                                                                      tCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTAc                                                                                                 <  2:413305/50‑1 (MQ=255)                                                                                                                      tCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTAc                                                                                                 <  1:215355/50‑1 (MQ=255)                                                                                                                      tCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATTCTTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAAt             >  1:1280026/1‑134 (MQ=255)                                                                                                                      tCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTcc        >  1:679641/1‑139 (MQ=255)                                                                                                                      tCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTcc        >  1:1479752/1‑139 (MQ=255)                                                                                                                      tCCATAATCAGTAAAATTGCTGCTTGATTAAACAAATGATACCTGATTACTGAAAGAGAGTTCCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAAt             <  2:1280026/134‑1 (MQ=255)                                                                                                                        cATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATTCCTGCg                                                                     >  1:39936/1‑76 (MQ=255)                                                                                                                        cATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATTCCTGCg                                                                     <  2:39936/76‑1 (MQ=255)                                                                                                                        cATAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCtt      >  2:1546770/1‑139 (MQ=255)                                                                                                                         aTAATCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCttt     <  1:1663221/139‑1 (MQ=255)                                                                                                                            aTCAGTAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCTTTGCt  <  2:492387/139‑1 (MQ=255)                                                                                                                                 |                                                                                                                                      CTACTGTCCACTGATAGCCCTGCGGTACGATTTTGACACCCGCGCCGACAATGACCAGCGCGACAAAAATGAGAATCGGGATAAAGATAAGCATCGGAAAAACCTCCTGTTGTACCGTCCATAATCAGCAAAATTGCTGCTTGATTAAACAAATTATACCTGATTACTGAAAGAGAGTTCCCCCTTATTCCTGCGAAGGATAAACTGTTTTTAGTAAAAATCAGAAAAAGGGAACAGCGATGCAGGAAAATAGTCCTTTGCT  >  NC_000913/515680‑515941

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