Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,816,007	A→T	S181S (TCT→TCA)	rph ←	truncated RNase PH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,816,007	0	A	T	98.2%	189.7 / NA	56	S181S (TCT→TCA)	rph	truncated RNase PH
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (38/17);  minor base G (1/0);  total (39/17)
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

GCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TT‑CCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTACCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCC  >  NC_000913/3815769‑3816250                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    gccctCGCTCCTCATGTTACTTTTCTACAAACAAAAAAAAGGCGACACATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCGTCTGCGTCGCTTCAATGGATTCGC‑TT‑CCCCCTCGGGCCAGCGTCAACAAGGTGCTTAGTGCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCaat                                                                                                                                                                                                                                          <  2:54483/247‑2 (MQ=255)       gCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCTACTCATCAGTCGCCTTATCAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACACTGGATTCGA‑TT‑CCCCCTCGGGCCAGCGCCCACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACGTCAATGATTCGCCCGTCTTCGGTCCTCACTACGTTCATGTCGGTCTCTGCGGCTTGGTCTTCAACGTCTTc                                                                                                                                                                                                                                    >  2:156377/1‑251 (MQ=255)          cctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TT‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTcc                                                                                                                                                                                                                                   >  1:366602/1‑249 (MQ=255)             cATCTGACTTCTCTAACGACAAAAAAAAGGCGACACAGCAGTCGCCTTAAAAATCACTTTGCCGGCGCCGCCTTCTGCGTCGCTACAATGGATTCGACTC‑CCCC‑GCTGGCCAGCGCCAACCAGTTGTGTAGCTCTTCATTGTTGAACGGCTCGCCCTCTGCCGTCCCCTGCACTTCAATTATGCGCCCGTCTTCGGTCATCACTACGCTCATGTCGGTCTCTGCGTCTGCGTCTTCAACGTATTCCAGagc                                                                                                                                                                                                                              <  2:20196/251‑3 (MQ=255)                  tACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TT‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCCGATCGCaa                                                                                                                                                                                                                          >  2:252826/1‑250 (MQ=255)                   aCTTTTCTACAGACAAAAAAAAGGCGACTCATCAGGCGCCTTAAAAAGCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TT‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCaaa                                                                                                                                                                                                                         <  1:137367/250‑1 (MQ=255)                      tttCTACAGACAAAAAAAAGGCGACTCATCAGACGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTAAAATGGATTCGA‑TT‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGCCTCTGCGGCTGAGTCTTCAACGTATTCCCGATCGCaaattgc                                                                                                                                                                                                                     >  2:161382/1‑247 (MQ=255)                                 aaataaaaGGCGACTAATCAGTCGCCTTTCACATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TTTCACA‑TCAGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGACCGGCTCGCCTTCTCCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTTCGTTCATGTCGGTCTCGGCGGCGGCGTCTGCCACGTATCCCAGATCGCaatttgcttcaacgttt                                                                                                                                                                                                           >  2:28341/5‑235 (MQ=255)                                    aaaaaGGCGACTCATCAGTCGCCTTACAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TT‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCTGTAATCACTACGTTCATGTCGGTCTCTGCGCCTGAGTCTTCAACGTATTCCAGCTCGctaatttcttcaccgttttct                                                                                                                                                                                                        >  2:256339/1‑230 (MQ=255)                                       atggCGATTCATCAGTAACCTTAAAGATCAGTTTGACAGCCACGAATTCTGCGTCGCTACAATGGATTCGA‑TT‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGCCGGTCTCTGCGGCTGAGTCTTCAACGTATTTCAGCTCGCAAATTGCTTCCCCGTTTACAAtac                                                                                                                                                                                                    >  2:273046/3‑249 (MQ=255)                                          gCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TT‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCg                                                                                                                                                                                                  >  1:350878/1‑250 (MQ=255)                                          gCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTAAAATGGATTCGA‑TT‑CCCCCTCGGGCCAGCGCCAACAAGGTAAGCAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCACTGATGCGCCCGTCTTCGGTCCTCACTACGTTCATTTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCg                                                                                                                                                                                                  >  1:245359/1‑250 (MQ=255)                                                catcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TT‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTGCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGa                                                                                                                                                                                             <  2:377798/249‑1 (MQ=255)                                                   caGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TT‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTg                                                                                                                                                                                        >  1:355915/1‑251 (MQ=255)                                                          cTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TT‑CCCCCTCGGGCCAGCGCCAACAAGGGGAGAAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTTCGGCTAc                                                                                                                                                                                 >  1:413564/1‑251 (MQ=255)                                                              aaaaTCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TT‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATc                                                                                                                                                                             >  1:284774/1‑251 (MQ=255)                                                                   cAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TT‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCACATTGCTTCCCCGTTTACAATACCCACAGAAACTGCGGCTACCATCCCttt                                                                                                                                                                        >  2:41603/1‑251 (MQ=255)                                                                         gccacccccgccTTCTGCGTCTCTACAATGGATTCGAGTC‑CCCC‑TCGGGCCAGCGCCAACAAGGTGGGGAGCTCTTCATGTGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGGCCGGCTTCGGGCATCACTACGTTCATGTCGGTCTCTGCTGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTGCACCGTTTACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCg                                                                                                                                                                   <  1:273046/243‑1 (MQ=255)                                                                              cgccgcCTTCTGCGTCGCTACAATGGATTCGA‑TT‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCACGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTgg                                                                                                                                                             >  1:102560/1‑251 (MQ=255)                                                                               gccgccTTCTGCGTCGCTACAATGGATTCGA‑TT‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGt                                                                                                                                                            >  1:451433/1‑251 (MQ=255)                                                                               gccgccTTCTGCGTCGCTACAATGGATTCGA‑TT‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGt                                                                                                                                                            <  1:3097/251‑1 (MQ=255)                                                                               accgccTTCTGCGTCGATACAATGGATTCGA‑TT‑CCCCCTCGGTCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCCATTATGCGCCCATCTTCGGTCATCACTACGTTCATGTCGATCTCTTCGGCTGCGCCTTCAACGTATTCCCGCTCGCACATTGCTTCACCTTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCCTCggttttgt                                                                                                                                                            >  2:63/2‑245 (MQ=255)                                                                                  gccTTCTGCGTCGCTACAATGGATTCGA‑TT‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGtttt                                                                                                                                                         >  1:58825/1‑251 (MQ=255)                                                                                      tcggCGTCGCTACAATGGATTCGA‑TT‑CCCCCGCGGGCCAGCGTCAACAAGGTGAGTAGTTCGTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAg                                                                                                                                                      <  2:333823/247‑1 (MQ=255)                                                                                       cTGCGTCGCTACAATGGATTCGA‑TT‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCACTTTCATCGGGTTGGTTTTAAGCt                                                                                                                                                    >  2:484868/1‑251 (MQ=255)                                                                                         gCGTCGCTACAATGGATTCGA‑TT‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATAGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTg                                                                                                                                                  >  2:54780/1‑251 (MQ=255)                                                                                          cGTCGCTACAATGGATTCGA‑TT‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGGCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCCACGTATTCCAGATCGCAAATTGCTTCACCGTGTACAATACCGACAGAAACTGCGGCTACCATCCCCTTCATCGGATTGGTGTTAAGCTTGc                                                                                                                                                 >  2:441607/1‑251 (MQ=255)                                                                                           gTCGCTACAATGGATTCGA‑TT‑CCCCCTCGGGCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTAAGCtttcc                                                                                                                                                >  2:336014/1‑248 (MQ=255)                                                                                                                          ggCCAGCGCCAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGCGTCTTCAACGTATTCCAGATCGCCAATTGCTTCACCGTTTACAATACCGGCCGAAACTGCGGCTACCATCCCTTTCCTCGGGTTGGTTTTTAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCtg                                                                                                                   >  2:3875/1‑250 (MQ=255)                                                                                                                                cgCCAACAAGGTGAGTAGCTCTTCATGGGTGACCGGCTCGCCTTCTGCCGTCCCCTGCCCTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGCCAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAgcg                                                                                                             <  2:384951/250‑1 (MQ=255)                                                                                                                                  ccAACAAGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGcc                                                                                                           >  2:343680/1‑251 (MQ=255)                                                                                                                                    aacaaGGTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCAc                                                                                                         >  1:359826/1‑251 (MQ=255)                                                                                                                                                    ctTCATGGGTGCACGGCTCGCCTTCTGCCGTCCCCTGCAATTCATTGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGACTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTAACCGTTTACAATACCGACAGAAACTGCGGCTACCATCACTTTCATCGGGTTGGTTTTCAGCTTGCCGTGTTCCACCCGATTCTGGAGCGNATCTTACAGCGCCACGCACGCCCCCGAAATc                                                                                         >  2:233229/1‑251 (MQ=255)                                                                                                                                                      tCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGAGGCGCCCGGCTTCGGGCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGTCAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCATGCACCCGTAAt                                                                                          <  2:463088/248‑1 (MQ=255)                                                                                                                                                       cATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCCCCCGTAATCGAc                                                                                      >  1:120078/1‑251 (MQ=255)                                                                                                                                                                  ggCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGt                                                                           >  1:135351/1‑251 (MQ=255)                                                                                                                                                                       gCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGccac                                                                      >  1:340805/1‑251 (MQ=255)                                                                                                                                                                              gCCGTGCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGCCGCGGTACGCGTGCCACCATCAGc                                                               >  1:318820/1‑251 (MQ=255)                                                                                                                                                                                 gTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCACCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTg                                                            >  2:53385/1‑251 (MQ=255)                                                                                                                                                                                        gctcTTCAATGATGCGCCCGGCTTGGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTCCGCGTGCCACCATCAGCCTGAAGCACg                                                     <  2:451433/248‑1 (MQ=255)                                                                                                                                                                                         cACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGt                                                    >  1:335636/1‑251 (MQ=255)                                                                                                                                                                                          aCTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTc                                                   >  2:263490/1‑251 (MQ=255)                                                                                                                                                                                                  gatgCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGGACGCGTGCCACCATCAGCCTGACGCACGTCGCAGTccc                                           >  2:258966/1‑250 (MQ=255)                                                                                                                                                                                                      cgcCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGt                                       <  1:144854/251‑1 (MQ=255)                                                                                                                                                                                                          cGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATg                                   <  1:12517/251‑1 (MQ=255)                                                                                                                                                                                                           gTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATgg                                  <  1:204313/251‑1 (MQ=255)                                                                                                                                                                                                             cTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTg                                >  1:124303/1‑251 (MQ=255)                                                                                                                                                                                                               tCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGa                               <  1:343680/250‑1 (MQ=255)                                                                                                                                                                                                                 ggTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACt                            >  2:15107/1‑251 (MQ=255)                                                                                                                                                                                                                    catcaCTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCAc                         <  1:198140/251‑1 (MQ=255)                                                                                                                                                                                                                           aCGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAgcgc                  >  2:483049/1‑251 (MQ=255)                                                                                                                                                                                                                               tCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCTACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCttt               <  1:72766/250‑1 (MQ=255)                                                                                                                                                                                                                               tCATGTCGGTCTCTGCGGCTGAGCCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTc              >  1:278440/1‑251 (MQ=255)                                                                                                                                                                                                                                      ggTCTCTGCGGCTGAGACTTCAACGTATTCCAGCTCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCCGCTTGCCGTTTTCCACCAGCTTCTGCAGCGCATCTGCCAGCGCCACGCCGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCCCGTCGCCGTCCAGCGCAATGGCTAACTCACCCAGCGCTTTCAAATCt        >  2:32553/1‑250 (MQ=255)                                                                                                                                                                                                                                           ctGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGAGAACTCACCCAGCGCTTTCAACTCTACTGcc  >  2:74056/1‑251 (MQ=255)                                                                                                                                                                                                                                             gCGGCTGATTCTTCAACGGATTCCAGAACGCAAATTGCTTCACCGTTTACAATACCGACAGAAACTGCGGCTACCATCCCTTTCATCGGGTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACGCATGCACCCGTAATCGACGCGGTACGCGTGCCACCAGCAGCCTGAAGCACGTCGCAGTCCAGCTTAATGGTTAACTCACCCAGCGCTTTCAAATCTACTGc   <  1:414811/248‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    GCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGA‑TT‑CCCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTACCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCC  >  NC_000913/3815769‑3816250

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