Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,204,860	A→G	T231A (ACC→GCC)  pseudogene (460/1137 nt)	ymfN → aaaE →	e14 prophage; chimeric replication protein/phage terminase YmfN phage terminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,204,860	0	A	G	100.0%	258.5 / NA	72	T231A (ACC→GCC)  pseudogene (460/1137 nt)	ymfN aaaE	e14 prophage; chimeric replication protein/phage terminase YmfN phage terminase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (36/36);  total (36/36)

TCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCGGCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGACCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTT  >  NC_000913/1204727‑1204996                                                                                                                                      |                                                                                                                                         tCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCGGCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGCTGAGCTGAGAGGATAAATCACTGGCCCTTg                                                                                                                                    >  1:804602/1‑140 (MQ=255) tCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCGGCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGGCCCTTg                                                                                                                                    >  2:1246323/1‑140 (MQ=255)   agcagcagcGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCGGCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGGCCCTTGAg                                                                                                                                  <  1:1246323/140‑1 (MQ=255)          gcGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCGGCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGGCCCTTGAGCAGTTCg                                                                                                                           >  1:287516/1‑140 (MQ=255)           cGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCGGCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGGCCCTTGAGCAGTTCGa                                                                                                                          >  2:979618/1‑140 (MQ=255)             tGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCGGCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGGCCCTTGAGCAGTTCGAgg                                                                                                                        >  2:240562/1‑140 (MQ=255)                gAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCGGCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGGCCCTTGAGCAGTTCGAGGGGc                                                                                                                     >  1:128850/1‑140 (MQ=255)                      tAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCGGCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGGCCCTTGAGCAGTTCGAGGGGCAGCCGt                                                                                                               <  2:287516/140‑1 (MQ=255)                      tAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCGGCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGGCCCTTGAGCAGTTCGAGGGGCAGCCGt                                                                                                               <  2:128850/140‑1 (MQ=255)                         cGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCGGCGTGTTCCGCGTATTTCAACCTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGGCCCTTGAGCAGTTCGAGGGGCAGCCGTGCa                                                                                                            <  1:1076191/140‑1 (MQ=255)                             cGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCGGCGCGTTCGGCGTAATTCAACCTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGGCCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCt                                                                                                        <  1:67346/140‑1 (MQ=255)                              gTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCGGCGCGTTCGGCGTATTTCGACCTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGGCCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTg                                                                                                       <  2:8855/140‑1 (MQ=255)                              gTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCGGCGCGTTCGGCGTATTCCAACCTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGGCCCTTGAGCAGTTCGAGGGGCAGCCGTg                                                                                                              <  2:37629/133‑1 (MQ=255)                              gTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCGGCGCGTTCGGCGTATTCCAACCTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGGCCCTTGAGCAGTTCGAGGGGCAGCCGTg                                                                                                              >  1:37629/1‑133 (MQ=255)                               tCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCGGCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTAGCCCTTGAGCATTTGGAGGGGAAGCCGTGCATTCTgg                                                                                                      >  1:501677/1‑140 (MQ=255)                                       aCGTCTTTAAAACAAATTACCTCAATATCTGGGCGTCGGCGCGTTCGGGGTATTTCAACCTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGGCCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGAc                                                                                              <  2:804602/140‑1 (MQ=255)                                            tttAAAACAAAACACCTCAATATCTGGGCGTCGGCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGGCCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGgc                                                                                         >  1:289169/1‑140 (MQ=255)                                             ttAAAACAAAACACCTCAATATCTGGGCGTCGGCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGGCCCTTGAGCAGTTCGAAGGGCAGCCGTGCATTAAGGCCTTTCACCTGgcg                                                                                        >  2:932346/1‑140 (MQ=255)                                               aaaacaaaacACCTCGATATCTAGGAGTCAGCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGGCCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGcgcg                                                                                      <  1:896393/140‑1 (MQ=255)                                               aaaacaaaacACCTCCATATCTAGGAGTCAGCGCGTTTGGCGTATTTCAACCTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGGCCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGcgcg                                                                                      <  1:596941/140‑1 (MQ=255)                                               aaaacaaaacACCTCAATATCTGGGCGTCGGCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGGCCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGcgcg                                                                                      <  1:385405/140‑1 (MQ=255)                                                    aaaacaCCTCAATATCTGGGCGTCGGCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGGCCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGc                                                                                 <  2:289169/140‑1 (MQ=255)                                                       acacCTCAATATCTGGGCGTCGCCGCGTTCGCCGTCTTTCAACCTGGTCAGCTGGCAGAGCTGCGAGGATAAATCACTGGCCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTgg                                                                              <  2:635176/140‑1 (MQ=255)                                                         cccTCAATATCTGGGCGACGGCACGTTCGGCTTATATCAACCTGGTGATCTGGCAGAGCTGCGAGGATAAATCACTGGCCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGat                                                                            <  2:501677/139‑1 (MQ=255)                                                           cTCAATATCTGGGCGTCGGCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGCAGAGCTGCGAGGATAAGTCACTGGCCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGc                                                                                 >  1:1752069/1‑133 (MQ=255)                                                           cTCAATATCTGGGCGTCGGCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGGCCCTTGAGCAGTTCGTGGGGCAGCTGTGGATTATGGCCTTTCACCTAGCGCGTAAGCTGGatat                                                                          >  2:600979/1‑140 (MQ=255)                                                           cTCAATATCTGGGCGTCGGCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGGCCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGc                                                                                 <  2:1752069/133‑1 (MQ=255)                                                                        cGTCGGCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGGCCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGgcg                                                             >  1:336336/1‑140 (MQ=255)                                                                        cGTCGGCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGGCCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGgcg                                                             >  1:1319791/1‑140 (MQ=255)                                                                        cGTCGGCGCGTTCGGCGAATTTCAACCTGGTGAGCAGGCCGAGCTGCGGGGATAAATCACTGGACCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGgcg                                                             <  2:735446/140‑1 (MQ=255)                                                                         gTCGGCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGGCCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGgcgc                                                            <  1:240562/140‑1 (MQ=255)                                                                         gTCGGCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGGCCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGgcgc                                                            <  2:336336/140‑1 (MQ=255)                                                                         gTCGGCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGGCCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGgcgc                                                            <  2:1319791/140‑1 (MQ=255)                                                                          tCGGCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGGCCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGcgcg                                                           >  2:374429/1‑140 (MQ=255)                                                                            ggCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGGCCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGAc                                                         >  1:996699/1‑140 (MQ=255)                                                                            ggCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGGCCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGAACTGGCGCGTAATCTGGATATGAACAGCATGGCGCGAc                                                         >  1:617173/1‑140 (MQ=255)                                                                              cgcgTTCGGCGTATTTCAACCTGGTGAGCTGGCGGAGCTGCGAGGATAAATCACTGGCCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGgcg                                                             <  2:781888/134‑1 (MQ=255)                                                                              cgcgTTCGGCGTATTTCAACCTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGGCCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGgcg                                                             >  1:781888/1‑134 (MQ=255)                                                                               gcgTTCGGCGTATTTCAACCTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGGCCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACttt                                                      <  2:996699/140‑1 (MQ=255)                                                                                 gTTCGGCGTATTTCAACCTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGGCCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTtat                                                    <  2:1406686/140‑1 (MQ=255)                                                                                       cGTATTTCAACCTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGGCCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGgcg                                                             >  2:1432292/1‑125 (MQ=255)                                                                                       cGTATTTCAACCTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGGCCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGgcg                                                             <  1:1432292/125‑1 (MQ=255)                                                                                               aaCCTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGGCCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGa                                      >  2:943183/1‑140 (MQ=255)                                                                                                aCCTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGGCCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGAt                                         <  1:471774/136‑1 (MQ=255)                                                                                                aCCTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGGCCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGAt                                         >  2:471774/1‑136 (MQ=255)                                                                                                 ccTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGGCCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACg                                    >  1:1623604/1‑140 (MQ=255)                                                                                                 ccTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGGCCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCAAGGCGCGACTTTATACCCGCGAGATTGACg                                    >  1:831919/1‑140 (MQ=255)                                                                                                    ggTGAGCTGGCAGAGCTGCGAGGATAAATCACTGGCCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTa                                 >  1:342982/1‑140 (MQ=255)                                                                                                     gTGAGCTGGCAGAGCTGCGAGGATAAATCACTGGCCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGa                                      <  1:117236/134‑1 (MQ=255)                                                                                                     gTGAGCTGGCAGAGCTGCGAGGATAAATCACTGGCCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGa                                      >  2:117236/1‑134 (MQ=255)                                                                                                      tGAGCTGGCAGAGCTGCGAGGATAAATCACTGGCCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTaaa                               >  2:1536577/1‑140 (MQ=255)                                                                                                      tGAGCTGGCAGAGCTGCGAGGATAAATCACTGGCCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTaaa                               >  2:241695/1‑140 (MQ=255)                                                                                                      tGAGCTGGCAGAGCTGCGAGGATAAATCACTGGCCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAACATGGCGCGACTTTTTACCCGAGAGATTGACGGTaaa                               >  1:833079/1‑140 (MQ=255)                                                                                                       gAGCTGGCAGAGCTGCGAGGATAAATCACTGGCCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTaaaa                              <  1:241695/140‑1 (MQ=255)                                                                                                       gAGCTGGCAGAGCTGCGAGGATAAATCACTGGCCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTaaaa                              <  2:342982/140‑1 (MQ=255)                                                                                                         gagacaggagCTGCGAGTATAAGTCACTGGCCCTTGAGCGGTTCGAGGTGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTAAAACg                            <  2:816006/133‑1 (MQ=255)                                                                                                          cTGGCAGAGCTGCGAGGATAAATCACTGGCCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTaaaa                              <  1:1150594/137‑1 (MQ=255)                                                                                                          cTGGCAGAGCTGCGAGGATAAATCACTGGCCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTaaaa                              >  2:1150594/1‑137 (MQ=255)                                                                                                           tGGCAGAGCTGCGAGGATAAATCACTGGCCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTAAAACGCa                          <  1:979618/140‑1 (MQ=255)                                                                                                             gCAGAGCTGCGAGGATAAATCACTGGCCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTAAAACGCAtt                        >  2:81142/1‑140 (MQ=255)                                                                                                               agagCTGCGAGGATAAATCACTGGCCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCCACTTTATTCCCGCGGGATTGACGGTAAAACGCATtac                      >  1:840591/1‑140 (MQ=255)                                                                                                                 agCTGCGAGGATAAATCACTGGCCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTAAAACGCATtacta                    <  1:1584550/140‑1 (MQ=255)                                                                                                                 agCTGCGAGGATAAAGCACTGGTCCCTGAGCAGTTCGTGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTAAAACGCATtacta                    <  1:600979/140‑1 (MQ=255)                                                                                                                   cTGCGAGGATAAATCACTGGCCCTTGAGCAGTTCGAGGGACAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTAAAACGCATTACTACa                  >  1:1445104/1‑140 (MQ=255)                                                                                                                       gAGGATAAATCACTGGCCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAgtgt              >  2:1615196/1‑140 (MQ=255)                                                                                                                       gAGGATAAATCACTGGCCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAgtgt              >  2:1331168/1‑140 (MQ=255)                                                                                                                           aTAAATCACTGGCCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCACGAGATTAACGCTAAAACGCATTACTACAgctgtctc          >  1:816006/1‑133 (MQ=255)                                                                                                                             aaaTCACTGGCCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTAAAACGCa                          <  2:1087359/122‑1 (MQ=255)                                                                                                                             aaaTCACTGGCCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTAAAACGCa                          >  1:1087359/1‑122 (MQ=255)                                                                                                                                tcccTGGCCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCcgcg     <  1:1660792/137‑1 (MQ=255)                                                                                                                                tCACTGGCCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCcgcg     <  1:81410/140‑1 (MQ=255)                                                                                                                                tCACTGGCCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCcgcg     <  1:1615196/140‑1 (MQ=255)                                                                                                                                   cTGGCCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGttt  <  2:1623604/140‑1 (MQ=255)                                                                                                                                      |                                                                                                                                         TCAGCAGCAGCGTGCGAAAAATAACGCCCGTCTGGCAAACGTCTTTAAAACAAAACACCTCAATATCTGGGCGTCGGCGCGTTCGGCGTATTTCAACCTGGTGAGCTGGCAGAGCTGCGAGGATAAATCACTGACCCTTGAGCAGTTCGAGGGGCAGCCGTGCATTCTGGCCTTTGACCTGGCGCGTAAGCTGGATATGAACAGCATGGCGCGACTTTATACCCGCGAGATTGACGGTAAAACGCATTACTACAGTGTGGCCCCGCGTTT  >  NC_000913/1204727‑1204996

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