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%	85.6 / NA	31	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 (19/12);  total (19/12)

AATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATCTCTGCTCTCA  >  NC_000913/3583326‑3583583                                                                                                                                 |                                                                                                                                  aaTGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATgg                                                                                                                         <  1:186871/139‑1 (MQ=255)                 ataatGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTc                                                                                                         >  2:456544/1‑139 (MQ=255)                      gatgatTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTg                                                                                                    <  2:735850/139‑1 (MQ=255)                           tttCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTc                                                                                               >  1:384091/1‑139 (MQ=255)                             tCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTg                                                                                             <  1:554893/139‑1 (MQ=255)                                   aCCGTGATCTACAAGAACGGTTAGTGATCGGGAATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTAAATGAATGATGGTGTTTTTGAGGTGCTCCTGTGGCTTCTGTTTCTa                                                                                       <  1:439602/139‑1 (MQ=25)                                      cTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCa                                                                                    >  1:348271/1‑139 (MQ=255)                                      cTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCa                                                                                    >  1:392978/1‑139 (MQ=255)                                         aTCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCt                                                                                 <  1:268661/139‑1 (MQ=255)                                              cAAGAACGGTTAGTGATCGGGGATTATAGCGTTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGc                                                                                  <  2:602311/133‑1 (MQ=255)                                              cAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGc                                                                                  >  1:602311/1‑133 (MQ=255)                                                aGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCcctc                                                                          <  2:721178/139‑1 (MQ=255)                                                   aCGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTg                                                                       >  1:771994/1‑139 (MQ=255)                                                      gTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTc                                                                    >  1:674880/1‑139 (MQ=255)                                                                 gggATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAg                                                                                                               <  1:693155/85‑1 (MQ=255)                                                                 gggATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAg                                                                                                               >  2:693155/1‑85 (MQ=255)                                                                 gggATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACg                                                         <  2:771994/139‑1 (MQ=255)                                                                    aTTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACgggg                                                      <  2:610430/139‑1 (MQ=255)                                                                       ataGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGggtgg                                                   >  2:345023/1‑139 (MQ=255)                                                                               tCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAAc                                           >  1:239998/1‑139 (MQ=255)                                                                                cAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACg                                          >  2:681807/1‑139 (MQ=255)                                                                                     tttACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCaaa                                     >  2:82770/1‑139 (MQ=255)                                                                                       tACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACgggg                                                      <  1:688256/120‑1 (MQ=255)                                                                                       tACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACgggg                                                      >  2:688256/1‑120 (MQ=255)                                                                                        aCTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGc                                  >  2:277654/1‑139 (MQ=255)                                                                                        aCTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTCCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTTACGGGCTGGTGCGTCACGGCAAAAGc                                  >  1:229915/1‑139 (MQ=14)                                                                                               aCATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCAccgccg                           >  2:608765/1‑139 (MQ=25)                                                                                                   tAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACa                       >  2:46386/1‑139 (MQ=17)                                                                                                   tAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACa                       >  2:384929/1‑139 (MQ=17)                                                                                                   tAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACa                       >  2:20869/1‑139 (MQ=17)                                                                                                                        acagATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTTACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATCTCTGCTCTCa  <  1:529171/136‑1 (MQ=0)                                                                                                                                 |                                                                                                                                  AATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATCTCTGCTCTCA  >  NC_000913/3583326‑3583583

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