Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,583,454	G→A	intergenic (+36/‑29)	yrhA → / → insA	pseudogene, interrupted by IS1E/IS1 repressor TnpA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,583,454	0	G	A	100.0%	67.8 / NA	24	intergenic (+36/‑29)	yrhA/insA	pseudogene, interrupted by IS1E/IS1 repressor TnpA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (14/10);  total (14/10)

GATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCC  >  NC_000913/3583350‑3583557                                                                                                         |                                                                                                        gatTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTCCTTATGACATTAAGGGTTATGCTGCCAACTTATTGATTTAATGTATGTTGGGTTTTTTGAGTTGTGCCAGTGGCt                                                                       >  1:474366/1‑139 (MQ=255)  atTTCATAAACCCTGTTATACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCtt                                                                      >  2:48270/1‑139 (MQ=255)  atTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCtt                                                                      <  2:312484/139‑1 (MQ=255)  atTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGATGTTTTTGAGGTTCTCCAGTGACtt                                                                      >  2:519273/1‑139 (MQ=255)   tttCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTc                                                                     <  2:240105/139‑1 (MQ=255)       aTAAACCCTGATCTACAAGAAAGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGtt                                                                 <  2:474366/139‑1 (MQ=255)                gATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGc                                                        >  2:222456/1‑139 (MQ=255)                 aTCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCt                                                       >  1:230591/1‑139 (MQ=255)                      cAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTccc                                                  >  2:304503/1‑139 (MQ=255)                            cGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCcctcc                                               >  2:158176/1‑136 (MQ=255)                            cGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCcctcc                                               <  1:158176/136‑1 (MQ=255)                                    gATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTAc                                    >  1:363421/1‑139 (MQ=255)                                    gATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTAc                                    <  2:609315/139‑1 (MQ=255)                                             ttATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGggt                           >  2:268753/1‑139 (MQ=255)                                              tataGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGggtg                          <  2:230591/139‑1 (MQ=255)                                                 aGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGgtggtg                       >  2:545720/1‑139 (MQ=255)                                                     tttCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTa                   <  1:304503/139‑1 (MQ=255)                                                     tttCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTa                   >  1:198785/1‑139 (MQ=255)                                                                aCTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAgg                                                                                    >  1:140194/1‑63 (MQ=255)                                                                aCTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAgg                                                                                    <  2:140194/63‑1 (MQ=255)                                                                aCTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGc        >  2:699247/1‑139 (MQ=255)                                                                      gACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCAccgcc  <  2:318214/139‑1 (MQ=37)                                                                           tAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACgggg                            >  2:102778/1‑108 (MQ=17)                                                                           tAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACgggg                            <  1:102778/108‑1 (MQ=17)                                                                                                         |                                                                                                        GATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCC  >  NC_000913/3583350‑3583557

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