Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,832,051	C→T	P532L (CCT→CTT)	yicH →	AsmA family protein YicH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,832,051	0	C	T	97.9%	152.3 / NA	48	P532L (CCT→CTT)	yicH	AsmA family protein YicH
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (27/19);  minor base G (0/1);  total (28/20)
Fisher's exact test for biased strand distribution p-value = 4.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACC‑GCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGC  >  NC_000913/3831804‑3832297                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                       gcgcTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACC‑GCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCAGTTGAAACTTAc                                                                                                                                                                                                                                                     >  1:59949/1‑251 (MQ=255)      gACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACC‑GCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACTTACGGttc                                                                                                                                                                                                                                                >  2:180339/1‑250 (MQ=255)              aCAGCAGCACGGTGAATATAAGCGAACTGAGCGCATTTACAGCAAAAGGCATTCTGGAAGCCACC‑GCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACTTACGGTTAGCGGGCa                                                                                                                                                                                                                                         >  1:330397/1‑250 (MQ=255)                  cagcaCGGTGAATATCAGCGAACTGAGTGCATTTACTGCAAAAGGCATTCTGGAAGCCACC‑GCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGTCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACTTACGGTTAGCGGGCAATTGc                                                                                                                                                                                                                                    >  2:327454/1‑251 (MQ=255)                          aaaaTATCGGCGAACTGATTGAATTTACTAAAAAAGGATTTCTGGAAGCCCCC‑GCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACGGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACTTACGGTTAGCGGGCAATTGCATGCCGTg                                                                                                                                                                                                                            <  1:246499/249‑1 (MQ=255)                                 cggaaaagtagggcatttgcataaaaaaGGCTGTCTG‑AAGCCACCCGCCAGTGTTACACAACCGCCACACCGTCAGACACATATCAGCCTGAATGGACGCGTTTTGCGGGTGAATATTTTGCAACAGTGGGGATGCCCTAAATTTCCGTTGACTGGCGAGGGCATTATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACTTACGGTTAGCGGGCAATTGCATGCCGTGAAcgccg                                                                                                                                                                                                                     <  2:254623/228‑1 (MQ=255)                                         tGAGTGCCTATACTGATAAAGGCATTCTGGAAGCCACC‑GCCAGCGTCTCACAAACGCAACATCGCCAGACACATATCAGCCTGAATGGACGCGGAGTGCCGGTGAATATTTTGCAACAGTGGGGATTGCCTGAATTACCGTTGACTGGCGACGGCACTATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAACCTTAAGTTTAGCCGGCAATTGCATGCCGTGAACGTCGCAAagc                                                                                                                                                                                                               >  2:130885/1‑249 (MQ=255)                                            gTGCATTTACTGAAAAAGGCATTCTGGAAGCCACC‑GCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGTACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGTCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAgcagca                                                                                                                                                                                                           >  2:88678/1‑250 (MQ=255)                                                   tACTGAAAAAGGCATTCTGGAAGCCACC‑GCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACTTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTc                                                                                                                                                                                                   >  2:234639/1‑251 (MQ=255)                                                                   cTGGAAGCCACC‑GCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACTTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAACGCAGCAAGTGACTCAAACCATGAATGCTgg                                                                                                                                                                                   >  1:42135/1‑251 (MQ=255)                                                                           cccc‑GCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGTCCCGGTGAATATTTTGAAACAGTGGGGTTGGCCTGAATTCCCGTTGACTGCGGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACTTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGttt                                                                                                                                                                           <  1:398110/249‑1 (MQ=255)                                                                            ccc‑GCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGGGCCGGTGAATATTTGGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACTTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTc                                                                                                                                                                          <  1:141481/250‑1 (MQ=255)                                                                                gccagtgtgtccaaacccccACATCGTCAGCCACATATCAGCCTGAAGGGACGCGGTGTGCCGGTGAATATTTTGCAACAGGGGGGATGGCCTGATTTCCCGTTGTCTGGCGCCGGCATTATTCAGCTTACCGCCAGTGGCGATAGTCAGGCCAAGGTCCCGTTGAAACTTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAg                                                                                                                                                                       <  2:54088/233‑1 (MQ=255)                                                                                gCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACTTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTcc                                                                                                                                                                         >  1:134929/1‑249 (MQ=255)                                                                                cccAGTGTTTCACAACCCCCACAACGGCAGACCCATATCAGCCTGAATGGACGCGGGTTCCCGGTGATTATTTTGCATCAGTGGGGATGGCTTGAATTTCCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACTTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAg                                                                                                                                                                       <  1:125347/250‑1 (MQ=255)                                                                                     tgtTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACTTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTg                                                                                                                                                                  >  2:28524/1‑251 (MQ=255)                                                                                        ttCACAAACGCACCAACGTAAACCACATAACAGCCTGAATGGACGCGGGGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTAACGGGGGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACTTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGaa                                                                                                                                                                <  2:238589/250‑1 (MQ=255)                                                                                                         gTCAGACACATCTCAGCCTGAATGGACGCGGTGAGACGGTGAATATTTTGCAACAGTGAGGATGGCCTGACTTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACTTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGGTACATCGACGGAGCCg                                                                                                                                              >  2:249198/1‑251 (MQ=255)                                                                                                               caaatatCAGCCTGAATGGACGCGGTGTGCCGGTAAATATTTTGCAACAGTGGGGATGGCCTGAATTCCCGTTGACGGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACTTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCgg                                                                                                                                        <  1:344380/248‑1 (MQ=255)                                                                                                                 caTATCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACTTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGt                                                                                                                                       >  1:209378/1‑250 (MQ=255)                                                                                                                    atCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGCCCAATGTCCCGTTGAAACTTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCggaaacc                                                                                                                                   >  2:36896/1‑246 (MQ=255)                                                                                                                       aGCCTGGATGGACGGGGTGTGCCGATGAATATTTTGCAACAGTGGGGAGGGCCTGAATACCCGTTGGCTGGGTGCGGCAATATTCAGCTTCCCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACTTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAACCCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGGCGGAGCCGGTGCGGTAAACATc                                                                                                                                <  2:22235/251‑1 (MQ=255)                                                                                                                       aGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACTTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATc                                                                                                                                >  2:328329/1‑251 (MQ=255)                                                                                                                                ggACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACTTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCa                                                                                                                                                                        >  2:249353/1‑202 (MQ=255)                                                                                                                                ggACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACTTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGcc                                                                                                                       >  1:20219/1‑251 (MQ=255)                                                                                                                                             cGGTGAATTTTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGGGACGGCAATTTTCAGCTTACCGCCAGTGGCGATATTCAGGGCAATGTCCCGTTGAAACTTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCGCGCCGGGTACGACTCAg                                                                                                          <  2:389710/251‑1 (MQ=255)                                                                                                                                               gTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACTTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCacac                                                                                                        >  2:48716/1‑248 (MQ=255)                                                                                                                                                    tatcttGCAACAGTGGGGCAGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCACTGGCGATATTCAGGCCAATGTCCCGTTGAAACTTACGGTTAGCGGGCAATTGAATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGCAAGCTGGCATCGTTTCTAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTAACACTCAGACGAAcc                                                                                                   >  2:144597/5‑251 (MQ=255)                                                                                                                                                      ttttGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACTTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCg                                                                                                 >  1:127713/1‑251 (MQ=255)                                                                                                                                                             acattGGGGATGGCCTAAATTCCCGTTGACGGGCGCCGCAAAGATTCCGCTTCCCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACTTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAAGCTCGCCGGGTACGACTCAGACGTACCCGGCATTcc                                                                                          <  2:20219/247‑1 (MQ=255)                                                                                                                                                                  ggggATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACTTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCa                                                                                      >  1:209946/1‑250 (MQ=255)                                                                                                                                                                            tGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACTTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACa                                                                            <  1:197590/250‑1 (MQ=255)                                                                                                                                                                               aTTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACTTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGAAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGt                                                                        >  2:179537/1‑251 (MQ=255)                                                                                                                                                                                 tACCGTTAACGGGCGTCGGCAATATTCAGTTAACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACTTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATTGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTa                                                                       <  1:231466/250‑1 (MQ=255)                                                                                                                                                                                  ccctttGACTGGCGACGGAAATTTTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACTTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAAt                                                                     <  2:93299/247‑1 (MQ=255)                                                                                                                                                                                         aCTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACTTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAgc                                                              >  2:307610/1‑251 (MQ=255)                                                                                                                                                                                                aCGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACTTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTcc                                                       >  1:96955/1‑251 (MQ=255)                                                                                                                                                                                                    cAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACTTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTg                                                   >  1:351590/1‑251 (MQ=255)                                                                                                                                                                                                      atatTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACTTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTg                                                   >  1:76503/1‑249 (MQ=255)                                                                                                                                                                                                                           gCGATATTCAGGCCAATGTCCCGTTGAAACTTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCt                            <  1:180339/251‑1 (MQ=255)                                                                                                                                                                                                                                      gCCAATGTCCCGTGCAAACGTACGGTTAGCGGTCAATTGCAGGCCGTGAACCCCGCAAAGCAGCAAGTGACTCTAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCTGTAAACATTAATCCCGCCGGGTACGCCTCAGACGTACCCGGCATTCCATCAATAGATTGCAACGTAATTGGCAGCGCATTTCTTTGAGGCTTTACCCCCAGTCTCTGCTCACTTTCAGcc                 <  2:181128/251‑1 (MQ=255)                                                                                                                                                                                                                                       ccAATGTCCCGTTGAAACTTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCt                <  2:66102/251‑1 (MQ=255)                                                                                                                                                                                                                                        cAATGTCCCGTTGAAACTTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTg               >  1:176225/1‑251 (MQ=255)                                                                                                                                                                                                                                          atttCCCGTTGAAACTTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAATCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACATAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCg             <  1:215464/248‑1 (MQ=255)                                                                                                                                                                                                                                                cGTTGAAACTTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACa                                                                                                                                  <  1:187231/128‑1 (MQ=255)                                                                                                                                                                                                                                                cGTTGAAACTTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACa                                                                                                                                  >  2:187231/1‑128 (MQ=255)                                                                                                                                                                                                                                                cGTTGAAACTTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGt       <  1:107855/251‑1 (MQ=255)                                                                                                                                                                                                                                                     aaaCTTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGc  >  1:366332/1‑251 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                       GCGCTGACCGCCAACAGCAGCACGGTGAATATCAGCGAACTGAGTGCATTTACTGAAAAAGGCATTCTGGAAGCCACC‑GCCAGTGTTTCACAAACGCCACAACGTCAGACACATATCAGCCTGAATGGACGCGGTGTGCCGGTGAATATTTTGCAACAGTGGGGATGGCCTGAATTACCGTTGACTGGCGACGGCAATATTCAGCTTACCGCCAGTGGCGATATTCAGGCCAATGTCCCGTTGAAACCTACGGTTAGCGGGCAATTGCATGCCGTGAACGCCGCAAAGCAGCAAGTGACTCAAACCATGAATGCTGGCATCGTTTCCAGCGGTGAAGTTACATCGACGGAGCCGGTGCGGTAAACATCAATCTCGCCGGGTACGACTCAGACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGC  >  NC_000913/3831804‑3832297

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