Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	85,481	T→C	R11G (AGA→GGA)	yabR ←	protein YabR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	85,481	0	T	C	96.7%	196.6 / ‑3.8	61	R11G (AGA→GGA)	yabR	protein YabR
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (33/26);  total (34/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

TGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCTGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGCAAGATCGGACGGTTAATGTGTTTTACACATTTTTTCCGTCAAACAGTGAGGCAGGCCATGGAGATGTTGTCTGGAGCCGAGATGGTCGTCCGATCGCTTATCGATCAGGGCGTTAAACAAGTATTCGGTTATCCCGGAGGCGCAGTCCTTGATATTTA  >  NC_000913/85238‑85730                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                          tGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCATTTTAATCCTCCGTAATgg                                                                                                                                                                                                                                                    >  1:406394/1‑251 (MQ=255) tGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATgg                                                                                                                                                                                                                                                    >  2:508207/1‑251 (MQ=255)        aaGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGAttt                                                                                                                                                                                                                                             <  1:711597/251‑1 (MQ=255)        aaGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGAttt                                                                                                                                                                                                                                             >  1:517025/1‑251 (MQ=255)                  gcgcGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCt                                                                                                                                                                                                                                   <  1:937851/251‑1 (MQ=255)                    gcgATAAAGGCCATCAGTGGATGGAAGAGCAATTAGACTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCttt                                                                                                                                                                                                                                 <  2:89001/251‑1 (MQ=255)                    gcgATAAAGGACATCAGTGGATGGAAGCGCAATTATGCACAATTTGCAAACTCTAACTTATTGCAGAATAGGTCAGACATGAATGTCTGTGTTATTCTGCATTTTTTATTGAATGTAGAAGTTGATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAAGAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTGAATCCTCCGTAATGGATGATTTTATCGTTTcatt                                                                                                                                                                                                                                 <  2:657735/251‑4 (MQ=255)                      gATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCACTTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTc                                                                                                                                                                                                                               >  1:918384/1‑251 (MQ=255)                             gCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTtc                                                                                                                                                                                                                        <  1:943031/251‑1 (MQ=255)                              ccATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTtcc                                                                                                                                                                                                                       <  1:666831/251‑1 (MQ=255)                               cATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCATCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTtcct                                                                                                                                                                                                                      <  2:435348/251‑1 (MQ=255)                               cATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATGTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTtcct                                                                                                                                                                                                                      <  2:124624/251‑1 (MQ=255)                                aTCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTtcctc                                                                                                                                                                                                                     <  2:115159/251‑1 (MQ=255)                                 tCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTtcctcc                                                                                                                                                                                                                    >  1:243183/1‑251 (MQ=255)                                        atggaAGAGCAATTAGTCTCAATTTGCAACCGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTAt                                                                                                                                                                                                             >  1:808599/1‑251 (MQ=255)                                          ggaAGAGCAAATAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTc                                                                                                                                                                                                           <  2:411362/251‑1 (MQ=255)                                                  aaTTAGTATCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTTTGGTTTACTCTGTACTTCTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGTTCTCTATTTGAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTCGTGGGATTCAGCTGATTTATTATCAACTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCCCCTTTCCTCCTGTTTATTCttatttcc                                                                                                                                                                                                   <  1:182033/251‑4 (MQ=255)                                                   aTTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTAccc                                                                                                                                                                                                  >  1:984188/1‑251 (MQ=255)                                                                        cTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGc                                                                                                                                                                             >  1:518733/1‑251 (MQ=255)                                                                          taCTGATTGGAGAATAGGTCAGAGATTAATGTCTTGTTTATTCGGCATTTTTTATTGAATGTAGAATTTGATTCTGAATGTGTGTGCTCTATAGGTTAGGATTAATTAAAAAAATAGAGAAATTGCTGGAAGTTGTGGGAGTCAGCCGATTGATTATCAATTTAATCCTCCGTAATGTAGGATTTTATCGTTGCTTTTCACTTTTCCGCCTGATTATACTTAGTTCCCCGTGTTTATGTCTCTGGCTg                                                                                                                                                                              <  2:105096/247‑1 (MQ=255)                                                                           aCTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCa                                                                                                                                                                           >  1:341245/1‑250 (MQ=255)                                                                               aTTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCGTTTCACCTTTCCTCCTGTTTATTCTTATTCCCCCGTGTTTATGTCTCTGGCTGCCAATTGc                                                                                                                                                                      <  1:724468/251‑1 (MQ=255)                                                                                   cAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTaa                                                                                                                                                                  <  1:189329/251‑1 (MQ=255)                                                                                   cAGAATAGGTCAGACATGAATGTATGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGCTTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCAGATTTATTATAAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTACTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGACAATTGCTTa                                                                                                                                                                   >  1:578553/1‑250 (MQ=255)                                                                                     gAATAGGTCAGTCATTAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTAGTTCTCTGGCTGCCAATTGCTTAAGc                                                                                                                                                                <  2:448814/251‑1 (MQ=255)                                                                                        tAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATAGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCTGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGCa                                                                                                                                                               >  2:460695/1‑249 (MQ=255)                                                                                             cAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGCAAGATCg                                                                                                                                                         <  2:867381/250‑1 (MQ=255)                                                                                                  aTGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGTATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGAATAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCATCCGATTTATTATTAATTTAATCCTCCTTAATGGATGATTTTATTGTTTCTGTTCACCTTTCCTCCTGTTTATTCTTTTTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGCAAGATCTGACGGt                                                                                                                                                   <  1:124802/251‑1 (MQ=255)                                                                                                            gggttATTCTGCATTTTTTATTGAAAGTAGAAGTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTACGTTGTGGGATTCAGCCGATTTATTATCCATTTAATCCTCCGTGATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGGGTTTATGTCTCTGTCTTCCCATTGCTTAAGCCAGCTCGGACGGTTAATGTGttt                                                                                                                                         >  2:689048/4‑251 (MQ=255)                                                                                                            ggTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGCAAGATCGGACGGTTAATGTGtt                                                                                                                                          <  1:117304/250‑1 (MQ=255)                                                                                                                             ttATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGCAAGATCGGACGGTTAATGTGTTTTACACATTTTTTCCGt                                                                                                                         <  1:653662/250‑1 (MQ=255)                                                                                                                               aTTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCAAATTGCTTAAGCAAGATCGGCCGGTTAATGTGTTTTACAACTTTTTTCCGTCa                                                                                                                       >  2:302714/1‑250 (MQ=255)                                                                                                                                   aaTGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGCAAGATCGGACGGTTAATGTGTTTTACACATTTTTTCCGTCAAAc                                                                                                                    >  1:34958/1‑249 (MQ=255)                                                                                                                                   aaTGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGCAAGATCGGACGGTTAATGTGTTTTACACATTTTTTCCGTCAAAc                                                                                                                    >  1:742949/1‑249 (MQ=255)                                                                                                                                              ttATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGCAAGATCGGACGGTTAATGTGTTTTACACATTTTTTCCGTCAAACAGTGAGGCAGGcc                                                                                                       >  1:484129/1‑251 (MQ=255)                                                                                                                                                aTTCTGAATGTGTAGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGCAAGATCGGACGGTTAATGTGTTTTACACATTTTTTCCGTCAAACAGTGAGGCAGGCCAt                                                                                                     >  1:699192/1‑251 (MQ=255)                                                                                                                                                   cTGAATGTGTGGTCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATAAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGTAAGATCGGACGGTTAATTTGTTTTACACCTTTTTTCCGTCATACAGTGAGGCAGGCCATggt                                                                                                  >  2:728283/1‑250 (MQ=255)                                                                                                                                                      aaTGTGTGGGCTCTCTATTTTAGGTTTAATTAAAAAAATAGAGAAATTGCTGTAAATTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGCAAGATCGGACGGTTAATGTGTTTTACACATTTTTTCCGTCAAACAGTGAGGCAGGCCATgg                                                                                                   >  1:860914/1‑247 (MQ=255)                                                                                                                                                           gtgGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGCAAGATCGGACGGTTAATGTGTTTTACACATTTTTTCCGTCAAACAGTGAGGCAGGCCATGGAGAtgttg                                                                                           <  1:99963/250‑1 (MQ=255)                                                                                                                                                                ctctctATTTTAGGATTAATTAAAAAAATAGCGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATACTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGAAAGATCGGCCGGTTAATGTGTTTTACACATTTTTTCCGTAAAACAGTGAGGCAGGCCATGGAGATGTTTTCTgg                                                                                      >  2:696088/1‑250 (MQ=255)                                                                                                                                                                     tATTTTAGGATTACTTAAAAAAATAGAGATATTGCTGAAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGATGTTTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCTCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGCAAGATCGGACGGTTAATGTGTTTTACACATTTTTTCCGTCAAACAGTGAGGCAGGCCATGGAGATGTTGTCTGTAGCCgg                                                                                >  1:297608/1‑250 (MQ=255)                                                                                                                                                                     tATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGCAAGATCGGACGGTTAATGTGTTTTACACATTTTTTCCGTCAAACAGTGAGGCAGGCCATGGAGATGTTGTCTGGAGCCg                                                                                 >  1:727578/1‑250 (MQ=255)                                                                                                                                                                          tAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGCTTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGCAAGATCGGACGGTTAATGTGTTTTACACATTTTTTCCGTCAAACAGTGAGGCAGGCCATGGCGATGTTGTCTGGAGCCGAGATgg                                                                           >  1:84705/1‑251 (MQ=255)                                                                                                                                                                                taattaaAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGCAAGATCGGACGGTTAATGTGTTTTACACATTTTTTCCGTCAAACAGTGAGGCAGGCCATGGAGATGTTGTCTGGAGCCGAGATGGTCGTcc                                                                     >  1:353310/1‑251 (MQ=255)                                                                                                                                                                                           aTAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTCTTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCCATTGCTTAAGCAAGATCGGACGGTTAATGTGTTTTACACATTTTTTCCGTCAAACAGTGAGGCAGGCCATGGAGATGTTGTCTGGAGCCGAGATGGTCGTCCGATCGCTTATc                                                          >  1:418163/1‑251 (MQ=255)                                                                                                                                                                                                   aTTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGCAAGATCGGACGGTTAATGTGTTTTACACAtttt                                                                                                                               <  2:45776/174‑1 (MQ=255)                                                                                                                                                                                                   aTTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGCAAGATCGGACGGTTAATGTGTTTTACACAtttt                                                                                                                               >  1:45776/1‑174 (MQ=255)                                                                                                                                                                                                     tGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCTGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGCAAGATCGGACGGTTAATGTGTTTTACACATTTTTTCCGTCAAACAGTGAGGCAGGCCATGGAGATGTTGTCTGGAGCCGAGATGGTCGTCCGATCGCTTATCGATCAGGGc                                                 <  1:727200/250‑1 (MQ=255)                                                                                                                                                                                                         gTAAGTTTTGGGATTCAGCCGATTTATTATCAATTTAATCCCCCGTAATGGAGGATGTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGCAAGATCGGACGGTTAATGTGTTTTACACATTTTTTCCGTCTAACAGTGAGGCAGGCCATGGAGATGTTGTCTGGAGCCGAGATGGTCGTCCGATCGCTTATCGATCAGGGCGTTa                                             <  2:740566/250‑1 (MQ=255)                                                                                                                                                                                                         gTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGCAAGATCGGACGGTTAATGTGTTTTACACATTTTTTCCGTCAAACAGTGAGGCAGGCCATGGAGATGTTGTCTGGAGCCGCGATGGTCGTCCGATCGCTTATCGATCATGGCGTTa                                             >  1:243946/1‑250 (MQ=255)                                                                                                                                                                                                              ttGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGCAAGATCGGACGGTTAATGTGTTTTACACATTTTTTCCGTCAAACAGTGAGGCAGGCCATGGAGATGTTGTCTGGAGCCGAGATGGTCGTCCGATCGCTTATCGATCAGGGCGTTAAACa                                         >  2:667873/1‑249 (MQ=255)                                                                                                                                                                                                                     aTTCAGCCGATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGCAAGATCGGACGGTTAATGTGTTTTACACATTTTTTCCGTCAAACAGTGAGGCAGGCCATGGAGATGTTGTCTGGAGCCGAGATGGTCGTCCGATCGCTTATCGATCAGGGCGTTAAACAAGTATTCgg                                <  2:98735/251‑1 (MQ=255)                                                                                                                                                                                                                         aGCCTATTTTTTACCAAGTTAATCCTCCTTAATGGAGGATTTTATCTTTTCTTTTCCCCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGCAAGATCGGACGGTTAATGTGTTTTACACATTTTTTCCGTCAAACAGTGAGGCAGGCCATGGAGATGTTGTCTGGAGCCGAGATGGTCGTCCGATCGCTTATCGATCAGGGCGTTAAACAATTATTCGGTTa                             <  2:266124/250‑1 (MQ=255)                                                                                                                                                                                                                             gATTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGCAAGATCGGACGGTTAATGTGTTTTACACATTTTTTCCGTCAAACAGTGAGGCAGGCCATGGAGATGTTGTCTGGAGCCGAGATGGTCGTCCGATCGCTTATCGATCAGGGCGTTAAACAAGTATTCGGTTATCCCg                        >  1:789715/1‑251 (MQ=255)                                                                                                                                                                                                                              aTTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCTCCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGCAAGATCGGACGGTTAATGTGTTTTACACATTTTTTCCGTCAAACAGTGAGGCAGGCCATGGAGATGTTGTCTGGAGCCGAGATGGTCGTCCGATCGCTTATCGATCAGGGCGTTAAACAAGTATTCGGTTATCCCgg                       <  1:473553/251‑1 (MQ=255)                                                                                                                                                                                                                              aTTTATTATCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGCAAGATCGGACGGTTAATGTGTTTTACACATTTTTTCCGTCAAACAGTGAGGCAGGCCATGGAGATGTTGTCTGGAGCCGAGATGGTCGTCCGATCGCTTATCGATCAGGGCGTTAAACAAGTATTCGGTTATCCCgg                       >  1:537071/1‑251 (MQ=255)                                                                                                                                                                                                                                 tattatCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCAGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGCAAGATCGGACGGTTAATGTGTTTTACACATTTTTTACGTAAAACAGTGAGGCATGCCGTGGATCTGTTGTCTGGCGACGAGATGGTCGTCCGATCACTTATCGATCAGGGCGTTACACACGTATTAGGTTACCACGGAgg                    >  2:214055/1‑251 (MQ=255)                                                                                                                                                                                                                                     atCAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGCAAGATCGGACGGTTAATGTGTTTTACACATTTTTTCCGTCAAACAGTGAGGCAGGCCATGGAGATGTTGTCTGGAGCCGAGATGGTCGTCCGATCGCTTATCGATCAGGGCGTTAAAAAAGTATTCGGTTATACCGGAgtagc                 >  1:829746/1‑246 (MQ=255)                                                                                                                                                                                                                                       cAATTTAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCGGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGCAAGATCGGACGGTTAATGTGTTTTACACATTTTTTCCGTCAAACAGTGAGGCAGGCCATGGAGATGTTGTCTGGAGCCGAGATGGTCGTCCGATCGCTTATCGATCAGGGCGTTAAACAAGTATTCGGTTATCCCGGAGGCGCAg               <  2:860914/250‑1 (MQ=255)                                                                                                                                                                                                                                          tttAATCCTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGCAAGATCGGACGGTTAATGTGTTTTACACATTTTTTCCGTCAAACAGTGAGGCAGGCCATGGAGATGTTGTCTGGAGCCGAGATGGTCGTCCGATCGCTTATCGATCAGGGCGTTAAACAAGTATTCGGTTATCCCGGAGGCGCAGTCCt           >  1:304819/1‑251 (MQ=255)                                                                                                                                                                                                                                                ccTCCGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGCAAGATCGGACGGTTAATGTGTTTTACACATTTTTTCCGTCAAACAGTGAGGCAGGCCATGGAGATGTTGTCTGGAGCCGAGATGGTCGTCCGATCGCTTATCGATCAGGGCGTTAAACAAGTATTCGGTTATCACGGAGGCGCAGTCCTTGatat     >  2:853766/1‑251 (MQ=255)                                                                                                                                                                                                                                                    tGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGCAAGATCGGACGGTTAATGTGTTTTACACATTTTTTCCGTCAAACAGTGAGGCAGGCCATGGAGATGTTGTCTGGAGCCGAGATGGCCGTCCGATCGCTTATCGATCAGGGCGTTAAACAAGTATTCGGTTATCCCGGAGGCGCAGTCCTTGATATTTa  <  1:428814/250‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                          TGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGATTTATTATCAATTTAATCCTCTGTAATGGAGGATTTTATCGTTTCTTTTCACCTTTCCTCCTGTTTATTCTTATTACCCCGTGTTTATGTCTCTGGCTGCCAATTGCTTAAGCAAGATCGGACGGTTAATGTGTTTTACACATTTTTTCCGTCAAACAGTGAGGCAGGCCATGGAGATGTTGTCTGGAGCCGAGATGGTCGTCCGATCGCTTATCGATCAGGGCGTTAAACAAGTATTCGGTTATCCCGGAGGCGCAGTCCTTGATATTTA  >  NC_000913/85238‑85730

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