Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,050,688	(C)5→4	intergenic (+158/+275)	insB‑4 → / ← cspH	IS1 protein InsB/CspA family protein CspH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,050,684	0	C	.	100.0%	188.3 / NA	44	intergenic (+154/+279)	insB‑4/cspH	IS1 protein InsB/CspA family protein CspH
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (28/16);  total (28/16)

CAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTTGTGAAATCTCCCTTTCCTCATACCGCTTGCACTATATACTGCGGGTAGGAAAGGCGCGCAGAGGGAAATATAAGATTGTTAACTATCCCCCTCTGGAAATACGACTTGTGAGCATAAGATAAATCTTAGCAAAATACAAATTAAAGATAAGGAAATTGTTACGAAAGCTATTAATTATTGTTGGTAATATAGTTTCAAGTGGTACTATAACTGGCTGTT  >  NC_000913/1050548‑1050816                                                                                                                                         |                                                                                                                                     cAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTTGTGAAATCTCCCTTTCCTCATACCGCTTGCACTATATACTGCGGGTAGGAAAGGCGCGCAGAGGGAAATATAAGATTGTTAACTATcccc                                                                                                                                   <  1:23497/140‑1 (MQ=255)   aTACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTTGTGAAATCTCCCTTTCCTCATACCGCTTGCACTATATACTGCGGGTAGGAAAGGCGCGCAGAGGGAAATATAAGATTGTTAACTATcccctc                                                                                                                                 >  2:431326/1‑138 (MQ=255)    tacgtacgCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTTGTGAAATCTCCCTTTCCTCATACCGCTTGCACTATATACTGCGGGTAGGAAAGGCGCGCAGAGGGAAATATAAGATTGTTAACTATcccctct                                                                                                                                >  2:101249/1‑137 (MQ=255)    tacgtacgCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTTGTGAAATCTCCCTTTCCTCATACCGCTTGCACTATATACTGCGGGTAGGAAAGGCGCGCAGAGGGAAATATAAGATTGTTAACTATcccctct                                                                                                                                <  1:462613/140‑4 (MQ=255)            cTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTTGTGAAATCTCCCTTTCCTCATACCGCTTGCACTATATACTGCGGGTAGGAAAGGCGCGCAGAGGGAAATATAAGATTGTTAACTAT‑CCCCTCTGGAAATAc                                                                                                                       >  2:346380/1‑140 (MQ=255)            cTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTTGTGAAATCTCCCTTTCCTCATACCGCTTGCACTATATACTGCGGGTAGGAAAGGCGCGCAGAGGGAAATATAAGATTGTTAACTAT‑CCCCTCTGGAAATAc                                                                                                                       >  1:504760/1‑140 (MQ=255)              tGAGATCAGGATTGTCTACTTTGAGTATAACATCCTTGTGAAATCTCCCTTTCCTCATACCGCTTGCACTATATACTGCGGGTAGGAAAGGCGCGCAGAGGGAAATATAAGATTGTTAACTAT‑CCCCTCTGGAAATACGa                                                                                                                     >  2:226497/1‑140 (MQ=255)               gagaTCAGGATTGTCTACTTTGAGTATAACATCCTTGTGAAATCTCCCTTTCCTCATACCGCTTGCACTATATACTGCGGGTAGGAAAGGCGCGCAGAGGGAAATATAAGATTGTTAACTAT‑CCCCTCTGGAAATACGAc                                                                                                                    >  2:220661/1‑140 (MQ=255)                        aTTGTCTACTTTGAGTATAACATCCTTGTGAAATCTCCCTTTCCTCATACCGCTTGCACTATATACTGCGGGTAGGAAAGGCGCGCAGAGGGAAATATAAGATTGTTAACTAT‑CCCCTCTGGAAATACGACTTGTGAGCa                                                                                                           <  1:226497/140‑1 (MQ=255)                        aTTGTCTACTTTGAGTATAACATCCTTGTGAAATCTCCCTTTCCTCATACCGCTTGCACTATATACTGCGGGTAGGAAAGGCGCGCAGAGGGAAATATAAGATTGTTAACTAT‑CCCCTCTGGAAATACGACTTGTGAGCa                                                                                                           <  2:356913/140‑1 (MQ=255)                            tCTACTTTGAGTATAACATCCTTGTGAAATCTCCCTTTCCTCATACCGCTTGCACTATATACTGCGGGTAGGAAAGGCGCGCAGAGGGAAATATAAGATTGTTAACTAT‑CCCCTCTGGAAATACGACTTGTGAGCa                                                                                                           <  1:338452/136‑1 (MQ=255)                            tCTACTTTGAGTATAACATCCTTGTGAAATCTCCCTTTCCTCATACCGCTTGCACTATATACTGCGGGTAGGAAAGGCGCGCAGAGGGAAATATAAGATTGTTAACTAT‑CCCCTCTGGAAATACGACTTGTGAGCa                                                                                                           >  2:338452/1‑136 (MQ=255)                             cTACTTTGAGTATAACATCCTTGTGAAATCTCCCTTTCCTCATACCGCTTGCACTATATACTGCGGGTAGGAAAGGCGCGCAGAGGGAAATATAAGATTGTTAACTAT‑CCCCTCTGGAAATACGACTTGTGAGCATAAGa                                                                                                      <  2:449182/140‑1 (MQ=255)                              tACTTTGAGTATAACATCCTTGTGAAATCTCCCTTTCCTCATACCGCTTGCACTATATACTGCGGGTAGGAAAGGCGCGCAGAGGGAAATATAAGATTGTTAACTAT‑CCCCTCTGGAAATACGACTTGTGAGCATAAGAt                                                                                                     <  2:316697/140‑1 (MQ=255)                                      gTATAACATCCTTGTGAAATCTCCCTTTCCTCATACCGCTTGCACTATATACTGCGGGTAGGAAAGGCGCGCAGAGGGAAATATAAGATTGTTAACTAT‑CCCCTCTGGAAATACGACTTGTGAGCATAAGATAAATCTTa                                                                                             >  2:57047/1‑140 (MQ=255)                                      gTATAACATCCTTGTGAAATCTCCCTTTCCTCATACCGCTTGCACTATATACTGCGGGTAGGAAAGGCGCGCAGAGGGAAATATAAGATTGTTAACTAT‑CCCCTCTGGAAATACGACTTGTGAGCATAAGATAAATCTTa                                                                                             >  2:364009/1‑140 (MQ=255)                                          aaCATCCTTGTGAAATCTCCCTTTCCTCATACCGCTTGCACTATATACTGCGGGTAGGAAAGGCGCGCAGAGGGAAATATAAGATTGTTAACTAT‑CCCCTCTGGAAATACGACTTGTGAGCATAAGATAAATCTTAGCaa                                                                                         >  2:874/1‑140 (MQ=255)                                               ccTTGTGAAATCTCCCTTTCCTCATACCGCTTGCACTATATACTGCGGGTAGGAAAGGTGCGCAGAGGGAAATATAAGATTGTTAACTAT‑CCCCTCTGGAAATACGACTTGTGAGCATAAGATAAATCTTAGCAAAATAc                                                                                    >  1:187882/1‑140 (MQ=255)                                                                 tCCTCATACCGCTTGCACTATATACTGCGGGTAGGAAAGGCGCGCAGAGGGAAATATAAGATTGTTAACTAT‑CCCCTCTGGAAATACGACTTGTGAGCATAAGATAAATCTTAGCAAAATACAAATTAAAGATAAGGaaa                                                                  <  1:369384/140‑1 (MQ=255)                                                                         ccGCTTGCACTATATACTGCGGGTAGGAAAGGCGCGCAGAGGGAAATATAAGATTGTTAACTAT‑CCCCTCTGGAAATACGACTTGTGAGCATAAGATAAATCTTAGCAAAATACAAATTAAAGATAAGGAAATTGTTACg                                                          <  1:516015/140‑1 (MQ=255)                                                                            cTTGCACTATATACTGCGGGTAGGAAAGGCGCGCAGAGGGAAATATAAGATTGTTAACTAT‑CCCCTCTGGAAATACGACTTGTGAGCATAAGATAAATCTTAGCAAAATACAAATTAAAGATAAGGAAATTGTTACGaaa                                                       >  2:20938/1‑140 (MQ=255)                                                                            cTTGCACTATATACTGCGGGTAGGAAAGGCGCGCAGAGGGAAATATAAGATTGTTAACTAT‑CCCCTCTGGAAATACGACTTGTGAGCATAAGATAAATCTTAGCAAAATACAAATTAAAGATAAGGAAATTGTTACGaaa                                                       >  2:297436/1‑140 (MQ=255)                                                                            cTTGCACTATATACTGCGGGTAGGAAAGGCGCGCAGAGGGAAATATAAGATTGTTAACTAT‑CCCCTCTGGAAATACGACTTGTGAGCATAAGATAAATCTTAGCAAAATACAAATTAAAGATAAGGAAATTGTTACGaaa                                                       >  1:12042/1‑140 (MQ=255)                                                                            cTTGCACTATATACTGCGGGTAGGAAAGGCGCGCAGAGGGAAATATAAGATTGTTAACTAT‑CCCCTCTGGAAATACGACTTGTGAGCATAAGATAAATCTTAGCAAAATACAAATTAAAGATAAGGAAATTGTTACGaaa                                                       >  1:210043/1‑140 (MQ=255)                                                                            cTTGCACTATATACTGCGGGTAAGAAAGGCGCGCAGAGGGAAATATAAGATTGTTAACTAT‑CCCCTCTGGAAATACGACTTGTGAGCATAAGATAAATCTTAGCAAAATACAAATTAAAGATAAGGAAATTGTTACGaaa                                                       >  1:247792/1‑140 (MQ=255)                                                                             ttGCACTATATACTGCGGGTAGGAAAGGCGCGCAGAGGGAAATATAAGATTGTTAACTAT‑CCCCTCTGGAAATACGACTTGTGAGCATAAGATAAATCTTAGCAAAATACAAATTAAAGATAAGGAAATTGTTACGAAAg                                                      >  1:292996/1‑140 (MQ=255)                                                                             ttGCACTATATACTGCGGGTAGGAAAGGCGCGCAGAGGGAAATATAAGATTGTTAACTAT‑CCCCTCTGGAAATACGACTTGTGAGCATAAGATAAATCTTAGCAAAATACAAATTAAAGATAAGGAAATTGTTACGAAAg                                                      <  1:196670/140‑1 (MQ=255)                                                                               gCACTATATACTGCGGGTAGGAAAGGCGCGCAGAGGGAAATATAAGATTGTTAACTAT‑CCCCTCTGGAAATACGACTTGTGAGCATAAGATAAATCTTAGCAAAATACAAATTAAAGATAAGGAAATTGTTACGAAAGCt                                                    >  1:324001/1‑140 (MQ=255)                                                                               gCACTATATACTGCGGGTAGGAAAGGCGCGCAGAGGGAAATATAAGATTGTTAACTAT‑CCCCTCTGGAAATACGACTTGTGAGCATAAGATAAATCTTAGCAAAATACAAATTAAAGATAAGGAAATTGTTACGAAAGCt                                                    >  1:289617/1‑140 (MQ=255)                                                                                cACTATATACTGCGGGTAGGAAAGGCGCGCAGAGGGAAATATAAGATTGTTAACTAT‑CCCCTCTGGAAATACGACTTGTGAGCATAAGATAAATCTTAGCAAAATACAAATTAAAGATAAGGAAATTGTTACGAAAGCTa                                                   <  2:151279/140‑1 (MQ=255)                                                                                cACTATATACTGCGGGTAGGAAAGGCGCGCAGAGGGAAATATAAGATTGTTAACTAT‑CCCCTCTGGAAATACGACTTGTGAGCATAAGATAAATCTTAGCAAAATACAAATTAAAGATAAGGAAATTGTTACGAAAGCTa                                                   <  2:504760/140‑1 (MQ=255)                                                                                         cTGCGGGTAGGAAAGGCGCGCAGAGGGAAATATAAGATTGTTAACTAT‑CCCCTCTGGAAATACGACTTGTGAGCATAAGATAAATCTTAGCAAAATACAAATTAAAGATAAGGAAATTGTTACGAAAGCTATTAattatt                                          <  1:410692/140‑1 (MQ=255)                                                                                           gCGGGTAGGAAAGGCGCGCAGAGGGAAATATAAGATTGTTAACTAT‑CCCCTCTGGAAATACGACTTGTGAGCATAAGATAAATCTTAGCAAAATACAAATTAAAGATAAGGAAATTGTTACGAAAGCTa                                                   <  1:511553/129‑1 (MQ=255)                                                                                           gCGGGTAGGAAAGGCGCGCAGAGGGAAATATAAGATTGTTAACTAT‑CCCCTCTGGAAATACGACTTGTGAGCATAAGATAAATCTTAGCAAAATACAAATTAAAGATAAGGAAATTGTTACGAAAGCTa                                                   >  2:511553/1‑129 (MQ=255)                                                                                           gCGGGTAGGAAAGGCGCGCAGAGGGAAATATAAGATTGTTAACTAT‑CCCCTCTGGAAATACGACTTGTGAGCATAAGATAAATCTTAGCAAAATACAAATTAAAGATAAGGAAATTGTTACGAAAGCTATTAATTAttgt                                        >  1:184807/1‑140 (MQ=255)                                                                                             gggTAGGAAAGGCGCGCAGAGGGAAATATAAGATTGTTAACTAT‑CCCCTCTGGAAATACGACTTGTGAGCATAAGATAAATCTTAGCAAAATACAAATTAAAGATAAGGAAATTGTTACGAAAGCTATTAATTAttgttg                                      >  2:104857/1‑140 (MQ=255)                                                                                             gggTAGGAAAGGCGCGCAGAGGGAAATATAAGATTGTTAACTAT‑CCCCTCTGGAAATACGACTTGTGAGCATAAGATAAATCTTAGCAAAATACAAATTAAAGATAAGGAAATTGTTACGAAAGCTATTAATTAttgttg                                      >  1:279227/1‑140 (MQ=255)                                                                                             gggTAGGAAAGGCGCGCAGAGGGAAATATAAGATTGTTAACTAT‑CCCCTCTGGAAATACGACTTGTGAGCATAAGATAAATCTTAGCAAAATACAAATTAAAGATAAGGAAATTGTTACGAAAGCTATTAATTAttgttg                                      >  1:231639/1‑140 (MQ=255)                                                                                               gTAGGAAAGGCGCGCAGAGGGAAATATAAGATTGTTAACTAT‑CCCCTCTGGAAATACGACTTGTGAGCATAAGATAAATCTTAGCAAAATACAAATTAAAGATAAGGAAATTGTTACGAAAGCTATTAATTATTGTTgg                                     <  2:430559/139‑1 (MQ=255)                                                                                               gTAGGAAAGGCGCGCAGAGGGAAATATAAGATTGTTAACTAT‑CCCCTCTGGAAATACGACTTGTGAGCATAAGATAAATCTTAGCAAAATACAAATTAAAGATAAGGAAATTGTTACGAAAGCTATTAATTATTGTTgg                                     >  1:430559/1‑139 (MQ=255)                                                                                               gTAGGAAAGGCGCGCAGAGGGAAATATAAGATTGTTAACTAT‑CCCCTCTGGAAATACGACTTGTGAGCATAAGATAAATCTTAGCAAAATACAAATTAAAGATAAGGAAATTGTTACGAAAGCTATTAATTATTGTTGGt                                    >  2:356513/1‑140 (MQ=255)                                                                                               gTAGGAAAGGCGCGCAGAGGGAAATATAAGAATGTTAACTAT‑CCCCTCTGGAAATACGACTTGTGAGCATAAGATAAATCTTAGCAAAATACAAATTAAAGATAAGGAAATTGTTACGAAAGCTATTAATTATTGTTGGt                                    >  1:213899/1‑140 (MQ=255)                                                                                                             cAGAGGGAAATATAAGATTGTTAACTAT‑CCCCTCTGGAAATACGACTTGTGAGCATAAGATAAATCTTAGCAAAATACGAATTAAAGATAAGGAAATTGTTACGAAAGCTATTAATTATTGTTGGTAATATAGTTTCAAg                      <  2:521674/140‑1 (MQ=255)                                                                                                             cAGAGGGAAATATAAGATTGTTAACTAT‑CCCCTCTGGAAATACGACTTGTGAGCATAAGATAAATCTTAGCAAAATACAAATTAAAGATAAGGAAATTGTTACGAAAGCTATTAATTATTGTTGGTAATATAGTTTCAAg                      <  1:96503/140‑1 (MQ=255)                                                                                                                agGGAAATATAAGATTGTTAACTAT‑CCCCTCTGGAAATACGACTTGTGAGCATAAGATAAATCTTAGCAAAATACAAATTAAAGATAAGGAAATTGTTACGAAAGCTATTAATTATTGTTGGTAATATAGTTTCAAGTgg                   <  1:63303/140‑1 (MQ=255)                                                                                                                  ggAAATATAAGATTGTTAACTAT‑CCCCTCTGGAAATACGACTTGTGAGCATAAGATAAATCTTAGCAAAATACAAATTAAAGATAAGGAAATTGTTACGAAAGCTATTAATTATTGTTGGTAATATAGTTTCAAGTGGTa                 >  1:192613/1‑140 (MQ=255)                                                                                                                              ttGTTAACTAT‑CCCCTCTGGAAATACGACTTGTGAGCATAAGATAAATCTTAGCAAAATACAAATTAAAGATAAGGAAATTGTTACGAAAGCTATTAATTATTGTTGGTAATATAGTTTCAAGTGGTACTATAACTGGCt     >  1:430818/1‑140 (MQ=255)                                                                                                                                 ttAACTAT‑CCCCTCTGGAAATACGACTTGTGAGCATAAGATAAATCTTAGCAAAATACAAATTAAAGATAAGGAAATTGTTACGAAAGCTATTAATTATTGTTGGTAATATAGTTTCAAGTGGTACTATAACTGGCTGtt  >  2:366341/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                     CAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTTGTGAAATCTCCCTTTCCTCATACCGCTTGCACTATATACTGCGGGTAGGAAAGGCGCGCAGAGGGAAATATAAGATTGTTAACTATCCCCCTCTGGAAATACGACTTGTGAGCATAAGATAAATCTTAGCAAAATACAAATTAAAGATAAGGAAATTGTTACGAAAGCTATTAATTATTGTTGGTAATATAGTTTCAAGTGGTACTATAACTGGCTGTT  >  NC_000913/1050548‑1050816

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