Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,126,894:1	+G	5.6%	coding (89/924 nt)	yceM →	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,126,894	1	.	G	5.6%	280.7 / 10.8	71	coding (89/924 nt)	yceM	putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base . (35/32);  new base G (2/2);  total (37/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.80e-01

TCAGTGGTGAGGTTGAAGATCAGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCCGCGTCGAAGCCCTGGAAATCGAAGTGGCAGAACTGAAACAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTG‑CGCAAAAAGCGTGGTTAC‑CGGTGCTGGCGGCAGCGTCTGACTGGAC‑G‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAGCTGGCGCATTCCTTATGCCGATTCGTTATCCAGCCTTGCCGCCAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGACGTGGTCAGTACGTTACTCAATGCGGGGGTACATGTCTGTGTCGATAAACCGCTGGCAGAAAATCTGCGCGATGCTGAACGGCTGGTGGA  >  NC_000913/1126656‑1127128                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                             tCAGTGGTGAGGTTGAAGATCAGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCCGCGTCGAAGCCCTGGAAATCGAAGTGGCAGAACTGAAACAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTG‑CGCAAAAAGCGTGGTTAC‑CGGTGCTGGCGGCAGCGTCTGACTGGAc                                                                                                                                                                                                                                               <  2:1391608/239‑1 (MQ=255)      ggggaagTTGAAGATCAGCTGGCGGTCACTGATATGTCGAACGCGGTTGACGGTGATTTACAGGACCGCGTCGAGGCCCTGGAAATCGAAGTGGCAGAACTGAAACAGCGTCTTGATTCGTTGCTGGCTCATCTGGGAGATTAAAGTGAAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTG‑CGCAAAAAGCGTGGTTAC‑CGGTGCTGGCGGCAGCGTCTGACTGGAC‑G‑t                                                                                                                                                                                                                                           <  2:1415190/230‑1 (MQ=255)       gtgAGGTTGAAGATCAGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCCGCGTCGAAGCCCTGGAAATCGAAGTGGCAGAACTGAAACAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAAATTACGTATCGGTGTAGTGGGATTAGGTGGCATTG‑CGCAAAAAGCGTGGTTAC‑CGGTGCTGGCGGCAGCGTCTGACTGGAC‑G‑TTACa                                                                                                                                                                                                                                       >  1:1073373/1‑239 (MQ=255)          aGGTTGAAGATCAGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCCGCGTCGAAGCCCTGGAAATCGAAGTGGCAGAACTGAAACAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTG‑CGCAAAAAGCGTGGTTAC‑CGGTGCTGGCGGCAGCGTCTGACTGGAC‑G‑TTACAAgg                                                                                                                                                                                                                                    >  1:979274/1‑239 (MQ=255)          aGGTTGAAGATCAGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCCGAGTCGAAGCCCTGGAAATCGAAGTGGCAGAACTGAAACAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGGGAAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTG‑CGCAAAAAGCGTGGTTAC‑CGGTGCTGGCGGCAGCGTCTGACTGTAC‑G‑TTACAAgg                                                                                                                                                                                                                                    >  1:807411/1‑239 (MQ=255)             ttGAAGATCAGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCCGCGTCGAAGCCCTGGAAATCGAAGTGGCAGAACTGAAACAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTG‑CGCAAAAAGCGTGGTTAC‑CGGTGCTGGCGGCAGCGTCTGACTGGAC‑G‑TTACAAGGa                                                                                                                                                                                                                                   <  2:174540/237‑1 (MQ=255)                      aGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCCGCGTCGAAGCCCTGGAAATCGAAGTGGCAGAACTGAAACAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTG‑CGCAAAAAGCGTGGTTAC‑CGGTGCTGGCGGCAGCGTCTGACTGGAC‑G‑TTACAAGGAGCCTGGTCg                                                                                                                                                                                                                          >  1:154708/1‑237 (MQ=255)                                                   gTTGACGGTGATTTACAGGCCCGCGTCGAAGCCCTGGAAATCGAAGTGGCAGAACTGAAACAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTG‑CGCAAAAAGCGTGGTTAC‑CGGAGCTGGCGGCAGCGTCTGACTGGAC‑G‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCACTTtgtg                                                                                                                                                                                          >  1:141859/1‑239 (MQ=255)                                                    ttGACGGTGATTTACAGGCCCGCGTCGAAGCCCTGGAAATCGAAGTGGCAGAACTGAAACAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTG‑CGCAAAAAGCGTGGTTAC‑CGGTGCTGGCGGCAGCGTCTGACTGGAC‑G‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGa                                                                                                                                                                                         >  1:377645/1‑239 (MQ=255)                                                        cGGTGATTTACAGGCCCGCTTCGAAGCCCTGGTAATCGAAGTGGCAGAACTGAAACAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAAATTACGTATCGGCGTAGTGGGGTTAGGTGGCATTG‑CGCAAAAAGCGTGGTTAC‑CGGTGCTGGCGGCAGCGTCTGACTGGAC‑G‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAg                                                                                                                                                                                      <  2:803913/238‑1 (MQ=255)                                                                tACAGGCCCGCGTCGAAGCCCTGGAAATCGAAGTGGCAGAACTGAAACAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAAATTACGTATCGGCGTAGTGGGATTAGGGGGCATTG‑CGCAAAAAGCGTGGTTAC‑CGGTGCTGGCGGCAGCGTCTGACTGGAC‑G‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAGCTGGCGCAt                                                                                                                                                                             <  2:1132645/239‑1 (MQ=255)                                                                     gCCCGCGTCGAAGCCCTGGAAATCGAAGAGGCAGAACTGAAACAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTG‑CGCAAAAAGCGTGGTTAC‑CGGTGCTGGCGGCAGCGTCTGACTGGAG‑G‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAGCTCGCGCATTCCtt                                                                                                                                                                        >  1:368542/1‑239 (MQ=255)                                                                       ccGCGTCGAAGCCCTGGAAATCGAAGTGGCAGAACTGAAACAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTG‑CGCAAAAAGCGTGGTTAC‑CGGTGCTGGCGGCAGCGTCTGACTGGAC‑G‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAGCTGGCGCATTCCTTAt                                                                                                                                                                      >  1:97573/1‑239 (MQ=255)                                                                          cgTCGAAGCCCTGGAAATCGAAGTGGCAGAACTGAAACAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTG‑CGCAAAAAGCGTGGTTAC‑CGGTGCTGGCGGCAGCGTCTGACTGGAC‑G‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAGCTGGCGCATTCCTTATGcc                                                                                                                                                                   >  1:651361/1‑239 (MQ=255)                                                                          cgTCGAAGCCCTGGAAATCGAAGTGGCAGAACTGAAACAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTG‑CGCAAAAAGCGTGGTTAC‑CGGTGCTGGCGGCAGCGTCTGACTGGAC‑G‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAGCTGGCGCATTCCTTATGcc                                                                                                                                                                   >  1:740511/1‑239 (MQ=255)                                                                           gTCGAAGCCCTGGAAATCGAAGTGGCAGAACTGAAACAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTG‑CGCAAAAAGCGTGGTTAC‑CGGTGCTGGCGGCAGCGTCTGACTGGAC‑G‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAGCTGGCGCATTCCTTATGCCg                                                                                                                                                                  <  2:1236241/239‑1 (MQ=255)                                                                                aGCCCTGGAAATCGAAGTGGCAGAACTGAAACAGCGTCTTGATTCGTTGCTGGCCCAGCTGGGAGATTAAAGTGAAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTG‑CGCAAAAAGCGTGGTTAC‑CGGTGCTGGCGGCAGCGTCTGACTGGAC‑G‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAGCTGGCGCATTCCTTATGCCGATTCg                                                                                                                                                             <  2:451537/239‑1 (MQ=255)                                                                                            cGAAGTGGCAGAACTGAAACAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTG‑CGCAAAAAGCGTGGTTAC‑CGGTGCTGGCGGCAGCGTGTGACTGGAC‑G‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAGCTGGCGCATTCCTTATGCCGATTCGTTATCCAGCCtt                                                                                                                                                 >  2:491671/1‑239 (MQ=255)                                                                                                           gAAACAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTG‑CGCAAAAAGCGTGGTTAC‑CGGTGCTGGCGGCAGCGTCTGACTGGAC‑G‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAGCTGGCGCATTCCTTATGCCGATTCGTTATCCAGCCTTGCCGCCAGTTGCGAt                                                                                                                                  >  2:1417950/1‑239 (MQ=255)                                                                                                                 gCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTG‑CGCAAAAAGCGTGGTTAC‑CGGTGCTGGCGGCAGCGTCTGACTGGAC‑G‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAGCTGGCGCATTCCTTATGCCGATTCGTTATCCAGCCTTGCCGCCCGTTGCGATGCGGtt                                                                                                                            >  1:87436/1‑239 (MQ=255)                                                                                                                       tGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTG‑CGCAAAAAGCGTGGTTAC‑CGGTGCTGGCGGCAGCGTCTGACTGGAC‑G‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAGCTGGCGCATTCCTTATGCCGATTCGTTATCCAGCCTTGCCGCCAGTTGCGATGCGGTTTTtgtg                                                                                                                      <  2:1061157/239‑1 (MQ=255)                                                                                                                       tGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTG‑CGCAAAAAGCGTGGTTAC‑CGGTGCTGGCGGCAGCGTCTGACTGGAC‑G‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAGCTGGCGCATTCCTTATGCCGATTCGTTATCCAGCCTTGCCGCCAGTTGCGATGCGGTTTTtgtg                                                                                                                      <  2:468290/239‑1 (MQ=255)                                                                                                                       tGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTG‑CGCAAAAAGCGTGGTTAC‑CGGTGCTGGCGGCAGCGTCTGACTGGAC‑G‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAGCTGGCGCATTCCTTATGCCGATTCGTTATCCAGCCTTGCCGCCAGTTGCGATGCGGTTTTtgtg                                                                                                                      <  2:979274/239‑1 (MQ=255)                                                                                                                          ttCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTG‑CGCAAAAAGCGTGGTTAC‑CGGTGCTGGCGGCAGCGTCTGACTGGAC‑G‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAGCTGGCGCATTCCTTATGCCGATTCGTTATCCAGCCTTGCCGCCAGTTGCGATGCGGTTTTTGTGCAt                                                                                                                   <  2:97573/239‑1 (MQ=255)                                                                                                                                         tCTGGGAGATTAAAGTGAAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTG‑CGCAAAAAGCGTGGTTAC‑CGGTGCTGGCGGCAGCGTCTGACTGGAC‑G‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAGCTGGCGCATTCCTTATGCCGATTCGTTATCCAGCCTTGCCGCCAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCgccag                                                                                                     >  1:422428/1‑238 (MQ=255)                                                                                                                                         tCTGGGAGATTAAAGTGAAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTG‑CGCAAAAAGCGTGGTTAC‑CGGTGCTGGCGGCAGCGTCTGACTGGAC‑G‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAGCTGGCGCATTCCTTATGCCGATTCGTTATCCAGCCTTGCCGCCAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCgccag                                                                                                     >  1:1009959/1‑238 (MQ=255)                                                                                                                                              gagaTTAAAGTGAAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTG‑CGCAAAAAGCGTGGTTAC‑CGGTGCTGGCGGCAGCGTCTGACTGGAC‑G‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAGCTGGCGCATTCCTTATGCCGATTCGTTATCCAGCCTTGCCGCCAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACtt                                                                                               >  1:1278024/1‑239 (MQ=255)                                                                                                                                              gagaTTAAAGTGAAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTG‑CGCAAAAAGCGTGGTTAC‑CGGTGCTGGCGGCAGCGTCTGACTGGAC‑G‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAGCTGGCGCATTCCTTATGCCGATTCGTTATCCAGCCTTGCCGCCAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACt                                                                                                <  1:828835/238‑1 (MQ=255)                                                                                                                                              gagaTTAAAGTGAAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTG‑CGCAAAAAGCGTGGTTAC‑CGGTGCTGGCGGCAGCGTCTGACTGGAC‑G‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAGCTGGCGCATTCCTTATGCCGATTCGTTATCCAGCCTTGCCGCCAGTTGCGATGCGGTTTTTGTGCACTCCAGCACCGCCAGCCACtt                                                                                               >  1:112143/1‑239 (MQ=255)                                                                                                                                                  ttAAAGTGAAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTG‑CGCAAAAAGCGTGGTTAC‑CGGTGCTGGCGGCAGCGTCTGACTGGAC‑G‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAGCTGGCGCATTCCTTATGCCGATTCGTTATCCAGCCTTGCCGCCAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGAc                                                                                           <  2:377645/239‑1 (MQ=255)                                                                                                                                                  ttAAAGTGAAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTG‑CGCAAAAAGCGTGGTTAC‑CGGTGCTGGCGGCAGCGTCTGACTGGAC‑G‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAGCTGGCGCATTCCTTATGCCGATTCGTTATCCAGCCTTGCCGCCAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGAc                                                                                           >  1:1343807/1‑239 (MQ=255)                                                                                                                                                  ttAAAGTGAAAAAATTACGTAGCGGCGGAGTGGGATTAGGGGGCATGG‑CGCAAAAAGCGTGGTTAC‑CGGTGCTGGCGGATGCGTCTGACTGGAC‑G‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAGCTGGCGCATTCCTTATGCCGATTCGTTATCCAGCCTTGCCGCCAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGAc                                                                                           <  2:1031104/239‑1 (MQ=255)                                                                                                                                                   tAAGGTGAAAAAATTGCGTTATGGCGCAGTGGGTTTA‑GTGGCAATGTAGGAAAAAG‑ATGGGTACACGGTCCTGGAGGCGGCTGTTGAC‑GGAG‑GTTTACAAGACGCCTGGTCGCCTTCAGCG‑GTTAAACCCCGCCCGATTTGTGAGAGTTATCGCTTTCCATACGCCGTTACTTAATCCACCCTGGCCGCCTGGTGCTATGTGGTTTTTGTGCATTCCAGCACCGCCAGCCACCTTGACg                                                                                          <  2:807411/239‑1 (MQ=255)                                                                                                                                                   tAAAGTGAAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTG‑CGCAAAAAGCGTGGTTAC‑CGGTGCTGGCGGCAGCGTCTGACTGGAC‑G‑TTACAAGGAGCTTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAGCAGGCGGATTCCTTCTGCCGATACGTTATGCAGCCTTGCCGCAAGTTCCGATGCGGTTTTAGTGCATTCCCGCACCGCCAGCCCACTTGACg                                                                                          >  2:788914/1‑239 (MQ=255)                                                                                                                                                   tAAAGTGAAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTG‑CGCAAAAAGCGTGGTTAC‑CGGTGCTGGCGGCAGCGTCTGACTGGAC‑G‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAGCTGGCGCATTCCTTATGCCGATTCGTTATCCAGCCTTGCCGCCAGTTGCGATGCGGTTTTtgtg                                                                                                                      >  1:1277552/1‑211 (MQ=255)                                                                                                                                                   tAAAGTGAAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTG‑CGCAAAAAGCGTGGTTAC‑CGGTGCTGGCGGCAGCGTCTGACTGGAC‑G‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAGCTGGCGCATTCCTTATGCCGATTCGTTATCCAGCCTTGCCGCCAGTTGCGATGCGGTTTTtgtg                                                                                                                      <  2:1277552/211‑1 (MQ=255)                                                                                                                                                   tAAAGTGAAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTG‑CGCAAAAAGCGTGGTTAC‑CGGTGCTGGCGGCAGCGTCTGACTGGAC‑G‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAGCTGGCGCATTCCTTATGCCGATTCGTTATCCAGCCTTGCCGCCAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGACg                                                                                          >  1:138227/1‑239 (MQ=255)                                                                                                                                                   tAAAGTGAAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTG‑CGCAAAAAGCGTGGTTAC‑CGGTGCTGGCGGCAGCGTCTGACTGGAC‑G‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAGCTGGCGAATTCCTTATGCCGATTCGTTATCCAGCCTTGCCGCCAGTTGCGAGGCGGTTTTTGTGCATTCCCACACCCCCAGCCACttgcacg                                                                                          >  2:788883/1‑234 (MQ=255)                                                                                                                                                   tAAAGGGAAAAAATTACGTATCGGCGTAGTGGGATTAGGAGGCATGG‑CGCAAAAAGCGTGGTTGC‑CGGTGCTGGCGGCAGCGTCTGACTGGAC‑G‑TTACAAGGAGCCTGGTCGCCTCC‑GCGCGCGAAAGCCCTGCCATTTTGTGAAAGCTGGCGCATTCCTTATGCCGATTCGTTATCCAGCCTTGCCGCCAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGACg                                                                                          <  2:141859/239‑1 (MQ=255)                                                                                                                                                        tGAAAAATTTACTTTGCGGGGTAGTGGGGTGAGGTGGCAGTC‑CGAAAAAAGGGTGTTTAC‑CGGGGCTGGGGGAGGCGGCTGCCTGGAC‑G‑TTACAAGGAGCCTGGTCGCTTAC‑GCGCGGGAAAGCCCTGCCAATTTGTGAAAGCTGGCGCATTCCTTATGCCGATTCGTTATCCAGCCTTGCCGCCAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGACGTGGt                                                                                      <  2:154708/238‑1 (MQ=255)                                                                                                                                                         gAAAAAATTTCGGATGGGCGTAGTGGGATTGGGGGGCTTTG‑CGCAAAAGCCGTGGTTAC‑TGGTGCTGGCGGCAGCGCCTGACTGGCC‑G‑TTGCAGGTAGCCTGGTCGCCTCG‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAGCTGGCGCATTCCTTATGCCGATTCGTTAGCAAGCCTTGCCGCCAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGACGTGGt                                                                                      <  2:793941/237‑1 (MQ=255)                                                                                                                                                         gAAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTG‑CGCAAAAAGCGTGGTTAC‑CGGTGCTGGCGGCAGCGTCTGACTGGAC‑G‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAGCTGGCGCATTCCTTATGCCGATTCGTTATCCAGCCTTGCCGCCAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGACGTGGt                                                                                      <  1:1417950/237‑1 (MQ=255)                                                                                                                                                         gAAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTG‑CGCAAAAAGCGTGGTTAC‑CGGTGCTGGCGGCAGCGTCTGACTGGAC‑G‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAGCTGGCGCATTCCTTATGCCGATTCGTTATCCAGCCTTGCCGCCAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGACGTGGTCa                                                                                    <  1:961991/239‑1 (MQ=255)                                                                                                                                                         gAAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTG‑CGCAAAAAGCGTGGTTAC‑CGGTGCTGGCGGCAGCGTCTGACTGGAC‑G‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAGCTGGCGCATTCCTTATGCCGATTCGTTATCCAGCCTTGCCGCCAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGACGTGGTCa                                                                                    <  2:302978/239‑1 (MQ=255)                                                                                                                                                         gAAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTG‑CGCAAAAAGCGTGGTTAC‑CGGTGCTGGCGGCAGCGTCTGACTGGAC‑G‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAGCTGGCGCATTCCTTATGCCGATTCGTTATCCAGCCTTGCCGCCAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGACGTGGTCa                                                                                    <  1:813246/239‑1 (MQ=255)                                                                                                                                                                              aGTGGGATTAGGTGGCATTG‑CGCAAAAAGCGTGGTTAC‑CGGTGCTGGCGGCAGCGTCTGACTGGAC‑G‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAGCTGGCGCATTCCTTATGCCGATTCGTTATCCAGCCTTGCCGCCAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGACGTGGTCAGTACGTTACTCAATGCgggg                                                                <  2:1427165/238‑1 (MQ=255)                                                                                                                                                                              aGTGGGATTAGGTGGCATTG‑CGCAAAAAGCGTGGTTAC‑CGGTGCTGGCGGCAGCGTCTGACTGGAC‑G‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAGCTGGCGCATTCCTTATGCCGATTCGTTATCCAGCCTTGCCGCCAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGACGTGGTCAGTACGTTACTCAATGCgggg                                                                <  1:287706/238‑1 (MQ=255)                                                                                                                                                                                 gggATTGGGTGGCTTGG‑CGCAAAAAGGGGGGTGAC‑CGGGGGGGGCGGCAGCGTCTGACTGGAC‑G‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAGCTGGCGCATTCCTTATGCCGATTCGTTATCCAGCCTTGCCGCCAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGACGTGGTCAGTACGTTACTCAATGCGGGGGTa                                                             <  2:740511/238‑1 (MQ=255)                                                                                                                                                                                    aTTAGGTGGCATTG‑CGCAAAAAGCGTGGTTAC‑CGGTGCTGGCGGCAGCGTCTGACTGGACGG‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAGCTGGCGCATTCCTTATGCCGATTCGTTATCCAGCCTTGCCGCCAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGACGTGGTCa                                                                                    <  2:749583/213‑1 (MQ=255)                                                                                                                                                                                    aTTAGGTGGCATTG‑CGCAAAAAGCGTGGTTAC‑CGGTGCTGGCGGCAGCGTCTGACTGGACGG‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAGCTGGCGCATTCCTTATGCCGATTCGTTATCCAGCCTTGCCGCCAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGACGTGGTCa                                                                                    >  1:749583/1‑213 (MQ=255)                                                                                                                                                                                          tGGCATTG‑CGAAAAAAGCTGGGTTCC‑CGGGGCGGGCGGCAGCGTCTGACTGGAC‑G‑TTACAAGGAGCCCGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAGCTGGCGCATTCCTTATGCCGATTCGTTATCCAGCCTTGCCGCCAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGACGTGGTCAGTACGTTACTCAATGCGGGGGTACATGTCtgt                                                    <  2:1173823/238‑1 (MQ=255)                                                                                                                                                                                             cATTG‑CGCAAAAAGCGTGGTTAC‑CGGTGCTGGCGGCAGCGTCTGACTGGAC‑G‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAGCTGGCGCATTCCTTATGCCGATTCGTTATCCAGCCTTGCCGCCAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGACGTGGTCAGTACGTTACTCAATGCGGGGGTACATGTCTGTGTCg                                                >  1:35313/1‑239 (MQ=255)                                                                                                                                                                                             cATTG‑CGCAAAAAGCGTGGTTAC‑CGGTGCTGGCGGCAGCGTCTGACTGGAC‑G‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAGCTGGCGCATTCCTTATGCCGATTCGTTATCCAGCCTTGCCGCCAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGACGTGGTCAGTACGTTACTCAATGCGGGGGTACATGTCTGTGTCg                                                >  1:418905/1‑239 (MQ=255)                                                                                                                                                                                                   cgcAAAAAGCGTGGTTAC‑CGGTGCTGGCGGCAGCGTCTGACTGGAC‑G‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAGCTGGCGCATTCCTTATGCCGATTCg                                                                                                                                                             <  2:990086/125‑1 (MQ=255)                                                                                                                                                                                                   cgcAAAAAGCGTGGTTAC‑CGGTGCTGGCGGCAGCGTCTGACTGGAC‑G‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAGCTGGCGCATTCCTTATGCCGATTCg                                                                                                                                                             >  1:990086/1‑125 (MQ=255)                                                                                                                                                                                                   cgcAAAAAGCGTGGTTAC‑CGGTGCTGGCGGCAGCGTCTGACTGGAC‑G‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAGCTGGCGCATTCCTTATGCCGATTCGTTATCCAGCCTTGCCGCAAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGACGTGGTCAGTACGTTACTCAATGCGGGGGTACATGTCTGTGTCGATaa                                            >  1:1126966/1‑238 (MQ=255)                                                                                                                                                                                                     cAAAAAGCGTGGTTAC‑CGGTGCTGGCGGCAGCGTCTGACTGGAC‑G‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAGCTGGCGCATTCCTTATGCCGATTCGTTATCCAGCCTTGCCGCCAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGACGTGGTCAGTACGTTACTCAATGCGGGGGTACATGTCTGTGTCGATAAAcc                                         >  1:694980/1‑239 (MQ=255)                                                                                                                                                                                                      aaaaaGCGGGGTTAC‑CGGTGCTGGCGGCAGCGTCGGACTGGAC‑G‑TTACAAGGAGCCGGTCCGCCTAC‑ACGCGCGAAAGCCCTGCCAATTTGGGAAAGCTGGCGCATTCCTTATGCCGATACGTTATCCAGCCTTGCCGCCAGTTGCGAGGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGACGTGGTCAGTACGTTACTCAATGCGGGGGTACATGTCTGTGTCGATAAACCg                                        <  2:1126051/239‑1 (MQ=255)                                                                                                                                                                                                        aaaGCGTGGTTAC‑CGGGGCTGGGGGCAGCGTCGGACGGGAC‑G‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAGCTGGCGCATTCCTTATGCCGATTCGTTATCCAGCCTTGCCGCCAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGACGTGGTCAGTACGTTACTCAATGCGGGGGTACATGTCTGTGTCGATAAACCGCt                                      <  2:651361/239‑1 (MQ=255)                                                                                                                                                                                                              tGGTTAC‑CGGTGCTGGCGGCAGCGTCTGACTGGACGG‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAGCTGGCGCATTCCTTATGCCGATTCGTTATCCAGCCTTGCCGCCAGTTg                                                                                                                                      >  1:302997/1‑138 (MQ=255)                                                                                                                                                                                                              tGGTTAC‑CGGTGCTGGCGGCAGCGTCTGACTGGACGG‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAGCTGGCGCATTCCTTATGCCGATTCGTTATCCAGCCTTGCCGCCAGTTg                                                                                                                                      <  2:302997/138‑1 (MQ=255)                                                                                                                                                                                                                       ggTGCTGGCGGCAGCGTCTGACTGGAC‑G‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAGCTGGCGCATTCCTTATGCCGATTCGTTATCCAGCCTTGCCGCCAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGACGTGGTCAGTACGTTACTCAATGCGGGGGTACATGTCTGTGTCGATAAACCGCTGGCAGAAAATCTgc                        >  1:18074/1‑239 (MQ=255)                                                                                                                                                                                                                         tGCTGGCGGCAGCGTCTGACTGGAC‑G‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAGCTGGCGCATTCCTTATGCCGATTCGTTATCCAGCCTTGCCGCCAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGACGTGGTCAGTACGTTACTCAATGCGGGGGTACATGTCTGTGTCGATAAACCGCTGGCAGAAAATCTgcgc                      >  1:1413780/1‑239 (MQ=255)                                                                                                                                                                                                                         tGCTGGCGGCAGCGTCTGACTGGAC‑G‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAGCTGGCGCATTCCTTATGCCGATTCGTTATCCAGCCTTGCCGCCAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGACGTGGTCAGTACGTTACTCAATGCGGGGGTACATGTCTGTGTCGATAAACCGCTGGCAGAAAATCTgcg                       <  2:1009959/238‑1 (MQ=255)                                                                                                                                                                                                                         tGCTGGCGGCAGCGTCTGACTGGAC‑G‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAGCTGGCGCATTCCTTATGCCGATTCGTTATCCAGCCTTGCCGCCAGTTGCGATGCGGTTTTTGTGCACTCCAGCACCGCCAGCCACTTTGACGTGGTCAGTACGTTACTCAATGCGGGGGTACATGTCTGTGTCGATAAACCGCTGGCAGAAAATCTgcg                       <  2:112143/238‑1 (MQ=255)                                                                                                                                                                                                                           cTGGCGGCAGCGTCTGACTGGAC‑G‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAGCTGGCGCATTCCTTATGCCGATTCGTTATCCAGCCTTGCCGCCAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGACGTGGTCAGTACGTTACTCAATGCGGGGGTACATGTCTGTGTCGATAAACCGCTGGCAGAAAATCTGCGCGa                    >  1:1224907/1‑239 (MQ=255)                                                                                                                                                                                                                           cTGGCGGCAGCGCCTGACTGGAC‑G‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAGCTGGCGCATTCCTTATGCCGATTCGTTATCCATCCTTGCCGCCAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTCTGACGTGTTCAGCACGTTACTCAACGCAGGGGTACCTGTCTTGGTCACTAAACCGCTGGCAGAACACCTGCGCGa                    >  2:1049552/1‑239 (MQ=255)                                                                                                                                                                                                                             ggcggcAGCGTCTGACTGGAC‑G‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGCAAGCTGGCGCATTCCTTATGCCGATTCGTTATCCAGCCTTGCCGCCAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGACGTGGTCAGTACGTTACTCAATGCGGGGGTACATGTCTGGGTCGATAAACCGCTGGCAGAAAATATGCGCGAt                   >  1:974113/1‑238 (MQ=255)                                                                                                                                                                                                                             ggcggcAGCGTCTGACTGGAC‑G‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAGCTGGCGCATTCCTTATGCCGATTCGTTATCCAGCCTTGCCGCCAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGACGTGGTCAGTACGTTACTCAATGCGGGGGTACATGTCTGTGTCGATAAACCGCTGGCACAAAATCTGCGGGAt                   >  2:820964/1‑238 (MQ=255)                                                                                                                                                                                                                              ggggcggtgtttGACTGGAC‑G‑TTGAGAGAGGCCTGGTCGCCTGC‑GGGCGCGAAAGCCTTGCCAAGTTGGGAAAGCTGGCGCATTCCTTATGCCGTTCGGTTAGCGCGCCTGGCCCCCAGTTGCGATGGGGTTTTTGAGAATTCCAGCACCGCCAGCCACTTTGACGTGGTCAGTACGTGACTCAATGCGGGGGTACATGTCTGTGTCGATAAACCGCTGTCAGAAAATCTGCGCGATg                  <  2:562097/227‑1 (MQ=255)                                                                                                                                                                                                                                               gAC‑G‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAGCTGGCGCATTCCTTATGCCGATTCGTTATCCAGCCTTGCCGCCAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGACGTGGTCAGTACGTTACTCAATGCGGGGGTACATGTCTGTGTCGATAAACCGCTGGCAGAAAATCTGCGCGATGCTGACCGGCTGGTGGa  >  2:500511/1‑237 (MQ=255)                                                                                                                                                                                                                                               gAC‑G‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAGCTGGCGCATTCCTTATGCCGATTCGTTATCCAGCCTTGCCGCCAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGACGTGGTCAGTACGTTACTCAATGCGGGGGTACATGTCTGTGTCGATAAACCGCTGGCAGAAAATCTGCGCGATGCTGAACGGCTGGTGGa  >  1:325350/1‑237 (MQ=255)                                                                                                                                                                                                                                                 c‑G‑TTACAAGGAGCCAGGTCGCATAC‑GCGCGCGAACGCCCTGCCAATTTGTGAAAGCTGGCGCATTCCTTATGCCGATTCGTTATCCAGCCTTGCCGCCAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGACGTGGTCAGTACGTTACTCAATGCGGGGGTACATGTCTGTGTCGATAAACCGCt                                      <  2:93603/199‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                             TCAGTGGTGAGGTTGAAGATCAGCCGGCGGTGACGGATATGTCGAACGCGGTTGACGGTGATTTACAGGCCCGCGTCGAAGCCCTGGAAATCGAAGTGGCAGAACTGAAACAGCGTCTTGATTCGTTGCTGGCCCATCTGGGAGATTAAAGTGAAAAAATTACGTATCGGCGTAGTGGGATTAGGTGGCATTG‑CGCAAAAAGCGTGGTTAC‑CGGTGCTGGCGGCAGCGTCTGACTGGAC‑G‑TTACAAGGAGCCTGGTCGCCTAC‑GCGCGCGAAAGCCCTGCCAATTTGTGAAAGCTGGCGCATTCCTTATGCCGATTCGTTATCCAGCCTTGCCGCCAGTTGCGATGCGGTTTTTGTGCATTCCAGCACCGCCAGCCACTTTGACGTGGTCAGTACGTTACTCAATGCGGGGGTACATGTCTGTGTCGATAAACCGCTGGCAGAAAATCTGCGCGATGCTGAACGGCTGGTGGA  >  NC_000913/1126656‑1127128

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑