Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,563,778	(G)6→5	coding (148/258 nt)	ytiA ←	uncharacterized protein YtiA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,563,778	0	G	.	97.8%	191.3 / ‑4.2	45	coding (148/258 nt)	ytiA	uncharacterized protein YtiA
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base . (28/16);  total (29/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTC‑ATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCACTCGTTGACCAATTACCGCCGCGCCAATTGCGCCAGCAAATTGCGCATCAGGATGGGTATTTACCGGCATTCGCAGGTGAGATTCCAGC  >  NC_000913/4563539‑4564027                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                           gCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGcc                                                                                                                                                                                                                                                 <  2:333501/250‑1 (MQ=255) gCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGcc                                                                                                                                                                                                                                                 >  1:364773/1‑250 (MQ=255)  ccACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTGTGTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGAGCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTACAATGCAGCGCGAGACACCAACCTGCTGGCAATAAATGGACTCAGTCAGCCCAGGGATCACTCATTGTATTTCCTCGCCGGGGG‑ACAGGTCATCCaa                                                                                                                                                                                                                                               <  2:84171/251‑1 (MQ=255)   cacaATAGCAAAAAGAATCACAATCATGGAGATGCGTTTACTTACCGAAATACTTTTTGTCACGGTAATCACATTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACACAGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCACTAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCaa                                                                                                                                                                                                                                               >  2:137050/1‑250 (MQ=255)    acaATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAAt                                                                                                                                                                                                                                             >  1:63440/1‑251 (MQ=255)     caATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCACTAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATg                                                                                                                                                                                                                                            >  1:343246/1‑251 (MQ=255)         ggCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACGTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCAt                                                                                                                                                                                                                                         <  2:131189/250‑1 (MQ=255)            aaaaaGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAAcc                                                                                                                                                                                                                                     <  2:345159/251‑1 (MQ=255)                   aTCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAgcc                                                                                                                                                                                                                              >  1:98735/1‑251 (MQ=255)                           aatGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGt                                                                                                                                                                                                                       >  1:165012/1‑250 (MQ=255)                             tGGAGATGCGTTTACTTAACGAAATACTTTTTGTACCGGCAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGAGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTaa                                                                                                                                                                                                                     >  1:276589/1‑250 (MQ=255)                                      gTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCCAATAAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTCACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCAAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACgggg                                                                                                                                                                                                            >  1:52568/1‑250 (MQ=255)                                               aCGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTg                                                                                                                                                                                                  <  1:79636/251‑1 (MQ=255)                                                 gAAATACTTTTTGACACGGCAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTCTACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGCTCACCCATTGTATTTCCTAGC‑AGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAACTCACGGGTAATATCGCTGAc                                                                                                                                                                                                >  2:279780/1‑251 (MQ=255)                                                                        cccTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATGTCGCTGACGCGAAAAGTC‑ATGCCCGctgcct                                                                                                                                                                        >  2:81899/1‑251 (MQ=255)                                                                                    aaTCAAAAAGGTATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTC‑ATGCCCGCTGCCTGCAACGCTTTGc                                                                                                                                                            >  2:136715/1‑251 (MQ=255)                                                                                               tatCGCGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTC‑ATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATg                                                                                                                                                 >  2:191768/1‑251 (MQ=255)                                                                                                    agaAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTC‑ATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGcc                                                                                                                                            <  1:192048/251‑1 (MQ=255)                                                                                                          ataCGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTC‑ATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCg                                                                                                                                       <  2:132895/250‑1 (MQ=255)                                                                                                              gTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCCTCCTGCTAGCAATAAATTGACTCAGTCAGCCCAGGGATCGCCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCCGCCGCCGCGTAAATCACGGGGAATATCGCTGCCGCGAAAAGTC‑ATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCgcacgggg                                                                                                                                    >  2:249383/1‑242 (MQ=255)                                                                                                              gTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTC‑ATGCCCGCTGCCTGCAACGCTGTGCGGGTTTGTATGACGCCGACCGTGGAg                                                                                                                                  <  2:418222/251‑1 (MQ=255)                                                                                                               gagctaCCTTAAGGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAAGTGTCTCAGTCAGCCCAGGGATCACCCATTGTATTGCCTCGCCGGGGT‑ACAGTTCAGCCAAATGCATAACCCGCGTCCGCCGCGAAAATCAGCGGTAAGATCGGTAACGCGAAA‑GTCCGGGCCCGCGGCCTGCAACGCTTTGCGGGTTGGTATGAGGcccacgctggagt                                                                                                                                 <  2:384208/246‑12 (MQ=255)                                                                                                                        tAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTC‑ATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGCGGGaaataaa                                                                                                                         >  2:266730/1‑250 (MQ=255)                                                                                                                            gAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGACCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTC‑ATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGaaataaaaata                                                                                                                     >  1:392612/1‑250 (MQ=255)                                                                                                                                   tAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTC‑ATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAActc                                                                                                              <  1:201692/250‑1 (MQ=255)                                                                                                                                      cgcgGAGAAGACCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGGCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTC‑ATGCCCGCTGCCTGCAAAGCTTTGAGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTc                                                                                                          <  2:337882/251‑1 (MQ=255)                                                                                                                                       gcgGAGCAACCCCGTAATGTGCAATGCAGCGCCAGACAACATCCTGCTCGCCATAACTTGACTCAGTCAGCCAAGGGATCACCCATTCTATTTCCTCGCCGGGGG‑ACAGGTCAGCCCAATGCATAACCCGCAGCCGACGCGTAAATCCCGGGGAATATCGCTGACGCGAAAAGTC‑ATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAACTAAAAATAAAAACTCTTTc                                                                                                          >  1:41180/1‑250 (MQ=255)                                                                                                                                            gaAACCCCGTAATGTGCAATGCAGCGCCAGACTCCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGTGGAATATCGCTGACGCGAAAAGTC‑ATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATTTCGCCGTCCGTGGAGTGAAATATAAATAAAAACTCTTTCATCGGc                                                                                                    <  1:137050/251‑1 (MQ=255)                                                                                                                                             aaaCCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTC‑ATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCg                                                                                                   >  2:316855/1‑251 (MQ=255)                                                                                                                                               aCCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTC‑ATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCgtg                                                                                                 >  1:404266/1‑251 (MQ=255)                                                                                                                                                ccccGTACTGTGCAATGCAGCGCCAGAAACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACACATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAACTGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTC‑ATGCCCGCCGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCgtgt                                                                                                >  2:212784/1‑251 (MQ=255)                                                                                                                                                             aaTGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTC‑ATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGGGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCACTCGTTg                                                                                    >  2:70981/1‑250 (MQ=255)                                                                                                                                                              aTGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTC‑ATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCACTCGTTGAc                                                                                  >  1:342924/1‑251 (MQ=255)                                                                                                                                                                   gcgcCAGACCCCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTC‑ATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGCGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCCCTCTTTgccccaat                                                                             >  2:351208/1‑244 (MQ=255)                                                                                                                                                                       cAGACAACATCCTGCCGAAAATACATTGACCCAGTCAGACCAGGGTTCACCCATTGTATTTCCTCGCCGGCGG‑ACAGGTCAGCCAAATGCATAACCCGCAGAAGCCGCGTCAATCACGCGGAATATAGCTGACACGAAAAGTC‑CTGCCCGCTGCCTGCAACGCTTTGCGGGCTTGTATGACGCCGACCGTGGCGAGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCACTCGTTGACACATTAccc                                                                         >  2:272107/1‑250 (MQ=255)                                                                                                                                                                           cacCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTC‑ATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCACTCGTTGACCAATTACcgccgc                                                                     >  1:331633/1‑251 (MQ=255)                                                                                                                                                                                  cggctgGCAATAAATTGCATCAGTCAGCCCAGGGATCACCCATTGTATTGCCCCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTC‑ATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCACTCGTTGACCAATTACCGCCGCGCCAAtt                                                              <  2:125824/249‑1 (MQ=255)                                                                                                                                                                                         cAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCGCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTC‑ATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCACTCGTTGACCAATTACCGCCGCGCCAATTGCGCCAg                                                       <  1:282511/251‑1 (MQ=255)                                                                                                                                                                                           aaaaatggtccacgccagCCCAGGGATCCCCCATTGTATTTCCTCGCCGGGGG‑ACAGGACAGCCAAATGCAAACCCCGCAGCCGCCGCGTAAATCCCGGGGAATATCGCTGACGCGAAAAGTC‑ATGCCCGCGGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCGCACTCGTTGACCAATTACCGCCGCGCCAATTGCGCCAGCa                                                     <  2:8366/236‑1 (MQ=255)                                                                                                                                                                                                            gcacAGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTC‑ATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCCCTCGTTTACCAATTACCACCGCCCCAATTGCGCCAGCCAATTGCGCATCCGGATg                                    >  2:80678/4‑251 (MQ=255)                                                                                                                                                                                                                aGGGATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGAAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTC‑ATGCCCGCTGCCTGCAACGCTTTGAGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCACTCGTTGACCAATTTCCGCCGCGGCAATTACGCCAGCACATTGCGCATAAgggtgggt                                 >  2:107430/1‑250 (MQ=255)                                                                                                                                                                                                                 gggATCACCCATTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATGTCGCTGACGCGAAAAGTC‑ATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCACTCGTTGACCAATTACCGCCGCGCCAATTGCGCCAGCAAATTGCGCATCAGGATGGGTAt                               <  1:81899/251‑1 (MQ=255)                                                                                                                                                                                                                           aTTGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTC‑ATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCACTCGTTGACCAATTACCGCCGCGCCAATTTCGCCAGCAAATTGCGCATCAGGCTGGGTATTTACCGGCAt                     >  1:193501/1‑251 (MQ=255)                                                                                                                                                                                                                            ttGTATTTCCTCGCCGGGGG‑ACAGGTCAGCCAAAAGCATAACCCGCACCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAATC‑ATGACCGCTGCCTGCAACGCTTTGCGGGGTTGTATGACGCCAACCGTGGAGAGAAATCAAAATAAAAACTCTTTCATGGGCGTGTCCTCACTCGTTGACCAATTACCGCCGCGCCAATTTCGCCCACAAATTTCGCAACAGGAGGGGTCTTTACCAGCAtt                    >  1:383546/1‑251 (MQ=255)                                                                                                                                                                                                                                     cTCGCCGGGGG‑ACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTC‑ATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCACTCGGTGACCAATTACCGCCGCGCCAATTGCGCCAGCAAATTGCGCATCAGGGTGGGGATTTACCGGCCTTCGCAGGTgg           >  2:139746/1‑250 (MQ=255)                                                                                                                                                                                                                                        cgccgggggACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATAGCGCGGACGCGAAAAGTC‑ATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCAGCGGCGTGTCCTCACTCGTTGACCAATTACCGCCGCGCCAATTGCGCCAGCAAATTGCGCATCAGGATGGGTATTTACCGGCATTCGCAGGTgaga         <  2:364773/247‑1 (MQ=255)                                                                                                                                                                                                                                            gggggACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAACTCACGGGGAATATCGCTGACGCGAAAAGTC‑ATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCACTCGTTTACCAATTACCGCCGCCCCAATTGCGCCAGCAAATTGCGCATCAAGATGGGTATTTACCGGCCTTCggcggcgagagtcc     >  1:63327/1‑238 (MQ=255)                                                                                                                                                                                                                                             ggggACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTC‑ATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCACTCGTTGACCAATTACCGCCGCGCCAATTGCGCCCGCAAATTGCGCATCCGGATGGGGCTTTACCGGCATTCGCAggggggcttccc    >  2:15703/1‑241 (MQ=255)                                                                                                                                                                                                                                               ggACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGCAAAGTC‑ATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCACTCGTTGACCAATTAACGCCGCGCCAATTGCGCCAGCAAATTGCGCATCAGGATGGGTATTTACCGGCATTCGCAGGGGAGATTcccgc  >  1:215387/1‑248 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                           GCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTC‑ATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCACTCGTTGACCAATTACCGCCGCGCCAATTGCGCCAGCAAATTGCGCATCAGGATGGGTATTTACCGGCATTCGCAGGTGAGATTCCAGC  >  NC_000913/4563539‑4564027

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