Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,636,290	T→C	G26G (GGA→GGG)	nohA ←	Qin prophage; putative prophage DNA‑packaging protein NohA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,636,290	0	T	C	100.0%	97.6 / NA	33	G26G (GGA→GGG)	nohA	Qin prophage; putative prophage DNA‑packaging protein NohA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (14/19);  total (14/19)

GTTCTTCAACTTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGA  >  NC_000913/1636163‑1636404                                                                                                                                |                                                                                                                   gTTCTTCAACTTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCCCCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATCCCCTGTTCCTGc                                                                                                        <  2:692360/140‑1 (MQ=255) gTTCTTCAACTTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATCCCCTGTTCCTGc                                                                                                        <  2:690902/140‑1 (MQ=37) gTTCTTCAACTTCCCGGAGCAGTTTTTCGTTATCAGTTTCAGCATCCCCTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATCCCCTGTTCCTGc                                                                                                        <  2:313283/140‑1 (MQ=11)  ttcttcAACTTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATCCCCTGTTCCTGcc                                                                                                       <  2:590359/140‑1 (MQ=37)   tcttcAACTTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGAAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATCCCCTGTTCCTGCCa                                                                                                      <  2:779577/140‑1 (MQ=255)      tcAACTTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATCCCCTGTTCCTGCCAGtt                                                                                                   <  2:63534/140‑1 (MQ=37)          cTTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATCCCCTGTTCCTGCCAGTTCTGa                                                                                               <  1:52213/140‑1 (MQ=37)            tCCCGGCGCAGCTTTTCGTTCTGAATTTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATCCCCTGTTCCTGCCAGTTCTGAAt                                                                                             <  1:466598/140‑1 (MQ=19)            tCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTAGCCACCGCATCGGAAAAAGGGCATCACCTGGTCCTGCCAGTTCTGAAt                                                                                             >  2:496265/1‑140 (MQ=11)               cGGCGCAGCTTTTCGTTCTCAATTTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCCCAGCAAGCGCATCCCCAGTTCCTGCCAGTTCTGAATGGt                                                                                          >  2:439420/1‑140 (MQ=11)                       cTTTTCGTTCTCAATTTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTGGCAGAACGGGCAGCCCCTGTTCCTGGCAGTACTGAATGGTACGGATAc                                                                                  >  2:314844/1‑140 (MQ=11)                                ctcAATTTCAGCATCCATTTCGTCATACCATATTACAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATCCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGa                                                                         <  2:457860/140‑1 (MQ=11)                                 tcAATTTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATCCCCTGTTCCTGCCATTTCGGACTGGTACGGATACTCGCACCGaa                                                                        >  1:319993/1‑140 (MQ=255)                                   aaTTTCAGCATCCCTGTCGGCATACGATCTTATAACGTCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATCCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGaaaa                                                                      <  1:392520/140‑1 (MQ=255)                                           cATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATCCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTcagcc                                                              >  2:529933/1‑140 (MQ=25)                                                   tCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATCCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCtt                                                      >  2:105048/1‑140 (MQ=37)                                                       cATACCATCTTATAACGCCGGCAGAGTCATAGAGCCCCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATCCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTttg                                                  <  2:319993/140‑1 (MQ=255)                                                              tCTTATAACGGCGGCAGTGTCACAAAGCACCCCATGACCCTTGCCACCTCCTCGCAGAACGGGCATCCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACtt                                           <  2:311032/140‑1 (MQ=11)                                                                        gcggcAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATCCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCAGCTTTTTGTAGACTTCCATTGTTCa                                 >  1:426564/1‑140 (MQ=16)                                                                         cggcAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATCCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCAt                                >  1:627463/1‑140 (MQ=18)                                                                             agagTCATAAAGCACCTCATTAGCCTTGCGACCGCCTGGCAGAACGGGCATCCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCa                            <  1:496265/140‑1 (MQ=255)                                                                               agTCATAAAGCACCTAATTACCCTTCCCACCGCCTCGCAGAACGGGCATCCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACg                          <  1:314844/140‑1 (MQ=255)                                                                                       aGCACCTCATTACCCCTGCCACCGCCTCGTAGAACGGCCATCCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCaaaaa                  <  1:293653/140‑1 (MQ=255)                                                                                           ccTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATCCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACaga              <  2:948110/140‑1 (MQ=255)                                                                                           ccTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATCCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACaga              >  1:789672/1‑140 (MQ=37)                                                                                           ccTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATCCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACaga              >  1:125200/1‑140 (MQ=37)                                                                                           ccTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATCCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACa                >  2:109122/1‑138 (MQ=37)                                                                                           ccTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATCCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACa                <  1:109122/138‑1 (MQ=37)                                                                                           ccTCATTACCCTTGCCACCGCCTCACAGAACGGGCATCCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACaga              >  2:843080/1‑140 (MQ=255)                                                                                             tCATTACCCTTGCCACCGCCTCGCAGAACGGGCATCCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAgaga            <  2:789672/140‑1 (MQ=255)                                                                                             tCATTACCCTTGCCACCGCCTCGCAGAACGGGCATCCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAgaga            <  1:105048/140‑1 (MQ=255)                                                                                                      ttGCCACCGCCTCGCAGAACGGGCATCCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCTTTGTTCAGTCCGCGGCCAAAAGCAGAGAAAGGAATCg   >  1:270236/1‑140 (MQ=11)                                                                                                      ttGCCACCGCCTCGCAGAACGGGCATCCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCC‑AAAACAGAGAAAGGAAACGa  >  2:867462/1‑140 (MQ=255)                                                                                                                                |                                                                                                                   GTTCTTCAACTTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGA  >  NC_000913/1636163‑1636404

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