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%	103.1 / NA	33	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 (18/15);  total (18/15)

GTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACAT  >  NC_000913/3583323‑3583563                                                                                                                                    |                                                                                                              gTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTatga                                                                                                        <  2:1241662/139‑1 (MQ=255)    aaTGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATgg                                                                                                     <  1:1432275/139‑1 (MQ=255)         aTTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGttt                                                                                                <  1:794644/139‑1 (MQ=255)                       atgatgATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAg                                                                                  >  2:1269337/1‑139 (MQ=255)                            gatTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCt                                                                             >  1:447248/1‑139 (MQ=255)                             atTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCtt                                                                            <  1:782627/139‑1 (MQ=255)                                    aaaCCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTc                                                                     >  2:431819/1‑139 (MQ=255)                                           gATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGc                                                              <  2:1033810/139‑1 (MQ=255)                                            aTCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCt                                                                    >  2:406401/1‑132 (MQ=255)                                            aTCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCt                                                                    <  1:406401/132‑1 (MQ=255)                                            aTCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCt                                                             <  2:447248/139‑1 (MQ=255)                                            aTCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCt                                                             >  2:1364639/1‑139 (MQ=255)                                                 cAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTg                                                                         >  2:1149578/1‑122 (MQ=255)                                                 cAAGAACGGTTAGTGATCGGGGATTAGAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTg                                                                         <  1:1149578/122‑1 (MQ=255)                                                  aaGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCcct                                                       >  2:400476/1‑139 (MQ=255)                                                  aaGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCcct                                                       >  2:12120/1‑139 (MQ=255)                                                      aCGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTg                                                                         >  1:1131395/1‑117 (MQ=255)                                                      aCGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTg                                                                         <  2:1131395/117‑1 (MQ=255)                                                         gTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTc                                                >  2:1340038/1‑139 (MQ=255)                                                              tGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTa                                           >  2:385069/1‑139 (MQ=255)                                                                   ggggATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACAGAc                                      >  1:1145000/1‑139 (MQ=255)                                                                    gggATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACg                                     <  1:922664/139‑1 (MQ=255)                                                                    gggATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACg                                     >  1:449621/1‑139 (MQ=255)                                                                    gggATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACg                                     <  1:385069/139‑1 (MQ=255)                                                                    gggATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTATTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACg                                     >  2:922664/1‑139 (MQ=255)                                                                                   cAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACAGACGGGGTGGTGCGTAACg                      <  2:1145000/139‑1 (MQ=255)                                                                                         ttACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCt                                                                             >  1:758704/1‑78 (MQ=255)                                                                                         ttACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCt                                                                             <  2:758704/78‑1 (MQ=255)                                                                                           aCTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGc              >  2:746043/1‑139 (MQ=255)                                                                                                tGACATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCAccgc         >  1:1370839/1‑139 (MQ=37)                                                                                                  aCATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCt                                                                    >  1:61973/1‑78 (MQ=255)                                                                                                  aCATTAAAGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCt                                                                    <  2:61973/78‑1 (MQ=255)                                                                                                       aaaGGTGATGCTGCCAACTTACTGATTTAATGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACAt  <  1:12120/139‑1 (MQ=14)                                                                                                                                    |                                                                                                              GTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACAT  >  NC_000913/3583323‑3583563

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