Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,635,023	T→C	A275A (GCA→GCG)	ydfN ←	Qin prophage; putative side tail fiber assembly protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,635,023	0	T	C	97.9%	146.3 / ‑4.3	47	A275A (GCA→GCG)	ydfN	Qin prophage; putative side tail fiber assembly protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (21/25);  total (21/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

ATGTTTTTGACGGTGTTTTCCGCGTTACCAGCAGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAATCG  >  NC_000913/1634916‑1635159                                                                                                            |                                                                                                                                         aTGTTTTTGACGGTGTTTTCCGCGTTACCAGCAGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTGCAGCAGTGCCGGa                                                                                                          <  1:532047/140‑1 (MQ=255)                             cagcagCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCtgat                                                                              <  1:31724/140‑1 (MQ=255)                                 agcGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGatgatg                                                                          >  1:817897/1‑140 (MQ=255)                                  gcGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGt                                                                         >  2:724519/1‑140 (MQ=255)                                  gcGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGt                                                                         <  2:305265/140‑1 (MQ=255)                                     ttAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCtt                                                                      >  1:94840/1‑140 (MQ=255)                                      tAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTc                                                                     <  1:437960/140‑1 (MQ=255)                                       aaCGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTcc                                                                    <  2:461590/140‑1 (MQ=255)                                        aCGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTccc                                                                   >  2:642476/1‑140 (MQ=255)                                        aCGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTccc                                                                   <  1:31133/140‑1 (MQ=255)                                             gATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCAt                                                              >  1:581638/1‑140 (MQ=255)                                                    tgtCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGCGCGTGTGGACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGt                                                       <  2:425495/140‑1 (MQ=17)                                                    tgtCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATTAGCACCATCTGATGATGTCTTCCCTGCATTACGAGt                                                       <  2:595971/140‑1 (MQ=17)                                                    tgtCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGt                                                       >  1:813228/1‑140 (MQ=255)                                                    tgtCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGt                                                       <  2:669379/140‑1 (MQ=255)                                                    tgtCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGt                                                       >  1:669379/1‑140 (MQ=255)                                                     gtCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTc                                                      >  2:203285/1‑140 (MQ=255)                                                      tCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCt                                                     >  1:562429/1‑140 (MQ=255)                                                       ccATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTg                                                    <  2:9226/140‑1 (MQ=255)                                                         aTGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGt                                                       <  2:333475/135‑1 (MQ=255)                                                         aTGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGt                                                       >  1:333475/1‑135 (MQ=255)                                                         aTGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGc                                                  <  1:724519/140‑1 (MQ=255)                                                              aaCCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTa                                             >  2:688668/1‑140 (MQ=255)                                                                ccAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTAcc                                           >  2:28068/1‑140 (MQ=255)                                                                   aTCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCt                                        >  1:131003/1‑140 (MQ=255)                                                                       cAACGGAGTGCGTATGAGCACCAATACCGACAGTACGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCttgtt                                    <  1:642476/140‑1 (MQ=18)                                                                        aaCGGAGTGCGTATGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCttgttg                                   <  1:663002/140‑1 (MQ=255)                                                                          cGGAGTGCGTATGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTtgtg                                 <  1:28068/140‑1 (MQ=255)                                                                          cGGAGTGCGTATGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTtgtg                                 <  2:832886/140‑1 (MQ=255)                                                                               tGCGTATGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCAt                            >  1:372786/1‑140 (MQ=21)                                                                                    aTGAGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATcc                       >  1:114748/1‑140 (MQ=21)                                                                                       aGCACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCAGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTACATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCcg                    >  2:974258/1‑140 (MQ=21)                                                                                         cACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGGGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCcgcg                  >  2:354576/1‑140 (MQ=21)                                                                                         cACCAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCcgcg                  >  2:163636/1‑140 (MQ=25)                                                                                           ccAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCt                <  1:427888/140‑1 (MQ=255)                                                                                           ccAATACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCt                <  1:163636/140‑1 (MQ=255)                                                                                             aaTACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTa              <  2:94840/140‑1 (MQ=255)                                                                                              aTACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCtt               >  1:537717/1‑138 (MQ=25)                                                                                              aTACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCtt               <  2:537717/138‑1 (MQ=255)                                                                                              aTACCGACAGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAg             >  2:997367/1‑140 (MQ=25)                                                                                                      aGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGaaa     <  2:817897/140‑1 (MQ=255)                                                                                                      aGTATGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGaaa     >  2:88384/1‑140 (MQ=37)                                                                                                         aTGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAATcg  <  1:203285/140‑1 (MQ=255)                                                                                                         aTGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAATcg  <  2:581638/140‑1 (MQ=255)                                                                                                         aTGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAATcg  <  2:96922/140‑1 (MQ=255)                                                                                                         aTGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAATcg  <  1:732008/140‑1 (MQ=255)                                                                                                         aTGCGCGTGTGCACCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTNAGATTTGAAATcg  <  1:839612/140‑1 (MQ=255)                                                                                                            |                                                                                                                                         ATGTTTTTGACGGTGTTTTCCGCGTTACCAGCAGCGTTAACGGTGATGGTGTGTCCATGTGAACCAATCGCAACGGAGTGCGTATGAGCACCAATACCGACAGTATGTGCATGCGCGCCTGCGCTTGCAGCAGTGCCGGACAGCGAGTGGGTATGAGCACCATCTGATGATGTCTTCCCTGCATTACGAGTCTGGCCACTACCGCTTGTTGTGCTCATAATCCCCGCGCTTAGATTTGAAATCG  >  NC_000913/1634916‑1635159

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: ‑