Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,913,058	T→C	N587D (AAC→GAC)	prc ←	tail‑specific protease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,913,058	0	T	C	97.6%	136.9 / ‑4.9	41	N587D (AAC→GAC)	prc	tail‑specific protease
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (24/16);  total (24/17)
Fisher's exact test for biased strand distribution p-value = 4.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TAAAGCGTTCGTTCAAACGCGCCAGACGCGTCGCATCATCTTCATTATTCTCTTTCTCACGCACAGCGTAATTCAGAGAAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTTCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTAAATCTCCTGATTTCACATAAGTCGCGGCATCAATGCTATCCCACGGCAGCGCG  >  NC_000913/1912923‑1913179                                                                                                                                        |                                                                                                                          tAAAGCGTTCGTTCAAACGCGCCAGACGCGTCGCATCATCTTCATTATTCTCTTTCTCACGCACAGCGTAATTCAGAGAAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTgg                                                                                                                       >  2:1016624/1‑140 (MQ=255)     gCGTTCGTTCAAACGCGCCAGACGCGTCGCATCATCTTCATTATTCTCTTTCTCACGCACAGCGTAATTCAGAGAAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACt                                                                                                                   >  2:1042268/1‑140 (MQ=255)        ttcgttcAAACGCGCCAGACGCGTCGCATCATCTTCATTATTCTCTTTCTCACGCACAGCGTAATTCAGAGAAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACt                                                                                                                   <  2:596361/137‑1 (MQ=255)        ttcgttcAAACGCGCCAGACGCGTCGCATCATCTTCATTATTCTCTTTCTCACGCACAGCGTAATTCAGAGAAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACt                                                                                                                   >  1:596361/1‑137 (MQ=255)         tcgttcAAACGCGCCAGACGCGTCGCATCATCTTCATTATTCTCTTTCTCACGCACAGCGTAATTCAGAGAAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACTCAgg                                                                                                               <  1:1004648/140‑1 (MQ=255)                               tCGCATCATCTTCATTATTCTCTTTCTCACGCACAGCGTAATTCAGAGAAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTa                                                                                         >  1:236628/1‑140 (MQ=255)                                cGCATCATCTTCATTATTCTCTTTCTCACGCACAGCGTAATTCAGAGAAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTAt                                                                                        <  2:999072/140‑1 (MQ=255)                                     catcTTCATTATTCTCTTTCTCACGCACAGCGTAATTCAGAGAAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTcc                                                                                   <  1:1016624/140‑1 (MQ=255)                                       tcTTCATTATTCTCTTTCTCACGCACAGCGTAATTCAGAGAAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCtt                                                                                 >  1:409990/1‑140 (MQ=255)                                           cATTATTCTCTTTCTCACGCACAGCGTAATTCAGAGAAACGATATTGCGTTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTTCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTcagc                                                                             <  2:389049/140‑1 (MQ=255)                                                ttCTCTTTCTCACGCACAGCGTAATTCAGAGAAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTc                                                                        >  1:803615/1‑140 (MQ=255)                                                ttCTCTTTCTCACGCACAGCGTAATTCAGAGAAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTc                                                                        <  2:236628/140‑1 (MQ=255)                                                 tctcTTTCTCACGCACAGCGTAATTCAGAGAAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTcc                                                                       <  1:1072155/140‑1 (MQ=255)                                                     tttCTCACGCACAGCGTAATTCAGAGAAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGtt                                                                   >  1:1025162/1‑140 (MQ=255)                                                     tttCTCACGCACAGCGTAATTCAGAGAAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGtt                                                                   >  1:180689/1‑140 (MQ=255)                                                      ttCTCACGCACAGCGTAATTCAGAGAAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTc                                                                  >  1:627870/1‑140 (MQ=255)                                                                   cGTAATTCAGAGAAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAGGGCCGTTaaa                                                     >  2:904418/1‑140 (MQ=255)                                                                   cGTAATTCAGAGAAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTaaa                                                     >  2:863104/1‑140 (MQ=255)                                                                   cGTAATTCAGAGAAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTaaa                                                     >  2:379519/1‑140 (MQ=255)                                                                   cGTAATTCAGAGAAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTaaa                                                     >  1:836806/1‑140 (MQ=255)                                                                         tCAGAGAAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTAAATCTCCt                                               <  1:863104/140‑1 (MQ=255)                                                                           agagAAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTAAATCTCCTGa                                             >  2:124352/1‑140 (MQ=255)                                                                            gagaAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTAAATCTCCTGAt                                            <  1:379519/140‑1 (MQ=255)                                                                             agaAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTAAAtc                                                   >  2:107279/1‑132 (MQ=255)                                                                             agaAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTAAAtc                                                   <  1:107279/132‑1 (MQ=255)                                                                             agaAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTAAATCTCCTGAtt                                           <  1:124352/140‑1 (MQ=255)                                                                                  cGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTAAATCTCCTGATTTcaca                                      >  1:674722/1‑140 (MQ=255)                                                                                         gcgcTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTAAATCTCCTGATTTCACATAAGTcg                               <  2:409990/140‑1 (MQ=255)                                                                                               gTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTAAATCTCCTGATTTCACATAAGTCGCGGCAt                         <  2:1025162/140‑1 (MQ=255)                                                                                               gTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTAAATCTCCTGATTTCACATAAGTCGCGGCAt                         <  2:438385/140‑1 (MQ=255)                                                                                                tCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTAAATCTCCTGATTTCACATAAGTCGCGGCATc                        <  1:1042268/140‑1 (MQ=255)                                                                                                     cATAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTAAATCTCCTGATTTCACATAAGTCGCGGCATCAATGc                   <  1:1124712/140‑1 (MQ=255)                                                                                                       tAGCGTTGAAGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTAAATCTCCTGATTTCACATAAGTCGCGGCATCAATGCTa                 >  1:564694/1‑140 (MQ=255)                                                                                                            ttGAAGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTAAATCTCCTGATTTCACATAAGTCGCGGCATCAATGCTATCCCa            >  1:337884/1‑140 (MQ=255)                                                                                                               aaGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTAAATCTCCTGATTTCACATAAGTCGCGGCATCAATGCTATCCCACgg         >  2:1122511/1‑140 (MQ=255)                                                                                                               aaGCGCGCGATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTAAATCTCCTGATTTCACATAAGTCGCGGCATCAATGCTATCCCACgg         >  2:669208/1‑140 (MQ=255)                                                                                                                  cgcgcgATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTAAATCTCCTGATTTCACATAAGTCGCGGCATCAATGCTATCCCACGGCAg      <  2:586528/140‑1 (MQ=255)                                                                                                                    cgcgATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTAAATCTCCTGATTTCACATAAGTCGCGGCATCAATGCTATCCCACGGCAgcg    >  1:867328/1‑140 (MQ=255)                                                                                                                    cgcgATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTAAATCTCCTGATTTCACATAAGTCGCGGCATCAATGCTATCCCACGGCAgcg    >  2:707024/1‑140 (MQ=255)                                                                                                                    cgcgATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTAAATCTCCTGATTTCACATAAGTCGCGGCATCAATGCTATCCCACGGCAgcg    >  2:1098517/1‑140 (MQ=255)                                                                                                                      cgATATCCTTCATGATGTCCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTAAATCTCCTGATTTCACATAAGTCGCGGCATCAATGCTATCCCACGGCAGagcg  >  2:181189/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                          TAAAGCGTTCGTTCAAACGCGCCAGACGCGTCGCATCATCTTCATTATTCTCTTTCTCACGCACAGCGTAATTCAGAGAAACGATATTGCGCTTGTCCTTCATAGCGTTGAAGCGCGCGATATCCTTCATGATGTTCTGGAACTCAGGATCTTTCGCGATACGCGCATTATGTTCCTTCAGCAGCTCCGGTTCAAAGGCCGTTAAATCTCCTGATTTCACATAAGTCGCGGCATCAATGCTATCCCACGGCAGCGCG  >  NC_000913/1912923‑1913179

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