Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,638,600	A→C	P141P (CCA→CCC)	creD →	putative inner membrane protein CreD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,638,600	0	A	C	95.6%	136.6 / NA	46	P141P (CCA→CCC)	creD	putative inner membrane protein CreD
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base C (20/23);  total (22/24)
Fisher's exact test for biased strand distribution p-value = 2.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

CATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGTTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAAC  >  NC_000913/4638354‑4638845                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       taTCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCttt                                                                                                                                                                                                                                                     <  2:161447/249‑1 (MQ=255)    cGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTAtt                                                                                                                                                                                                                                                  <  2:177790/250‑1 (MQ=255)    cGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCACTTTAtt                                                                                                                                                                                                                                                  <  2:488695/250‑1 (MQ=255)     gCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGCGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATtgt                                                                                                                                                                                                                                                >  2:370931/1‑251 (MQ=255)          tCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGAtt                                                                                                                                                                                                                                            <  2:18549/250‑1 (MQ=255)                        tttATACGGTGCAGGACGAGGATAAAACCGTGTAGCGGAACCGCTGTTTTATCCAATTTTGGTTAACTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGACATTTATTGTGATTCGCGTCGGGGATgcg                                                                                                                                                                                                                             >  2:213149/1‑251 (MQ=255)                            taCGGTGCAGGCAGAGGATAAAACCGTGGAGAGGAAACGAAGTTTTCTCCATTTTTGGTTACCTGAGTCATTGCTGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGATAAGGACTGGCAAAGTGATTTAACGCTCAAAGCCGATTTCACTGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCCAGCCCTTTCTTGAGATTAGCGTCGGGAATGCGCGTg                                                                                                                                                                                                                         >  2:572509/1‑251 (MQ=255)                                         gagtaTAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑tggtgtggt                                                                                                                                                                                                           <  1:126207/247‑1 (MQ=255)                                                     gTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGa                                                                                                                                                                                               <  1:263800/251‑1 (MQ=255)                                                        gAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCTATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCTTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGt                                                                                                                                                                                            <  2:496634/251‑1 (MQ=255)                                                                aCGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAgc                                                                                                                                                                                                     <  2:506901/234‑1 (MQ=255)                                                                aCGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAgc                                                                                                                                                                                                     >  1:506901/1‑234 (MQ=255)                                                                  gAAGTTTTATCCATTTTTCGTTACCTGAGCCATTGATGCTTGATGACAATCAGCACGTGGACGAACGCAAGCTAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTCT‑TGGTGTGGGGAAAGCGCCTGAAGTTaacggaccg                                                                                                                                                                                  >  2:93339/1‑248 (MQ=255)                                                                         tATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTg                                                                                                                                                                            <  1:468220/250‑1 (MQ=255)                                                                           tCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGAcc                                                                                                                                                                         >  1:130359/1‑251 (MQ=255)                                                                              aTTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAGTCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCAtt                                                                                                                                                                      <  2:113067/251‑1 (MQ=255)                                                                                ttttGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGAGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGAGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTg                                                                                                                                                                     >  1:362060/1‑250 (MQ=255)                                                                                   tGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAAcc                                                                                                                                                                 >  1:337832/1‑251 (MQ=255)                                                                                   tGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAAc                                                                                                                                                                  <  1:107097/250‑1 (MQ=255)                                                                                     gTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCg                                                                                                                                                               <  1:308262/251‑1 (MQ=255)                                                                                        aCCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTATGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAt‑t                                                                                                                                                                                                                   <  1:383953/196‑1 (MQ=255)                                                                                        aCCTGAGTAATTGATGGTTGATGGCAATCAGCACGTGGCAGAACGCAAGATAGGGCTTTATACCGGTCAGGTCTGGCACAGAGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCACATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGGAt‑t                                                                                                                                                                                                                   >  2:383953/1‑196 (MQ=255)                                                                                                     atggttgatggCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTTATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGGGGTGAAAGCGACTGAAGTTAACGGAACGGCGCt                                                                                                                                                                             <  1:203069/221‑1 (MQ=255)                                                                                                                  aTCAGAACATGGAAGAACGCAAGATAGGGCTTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAACGCCGCTTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGAAAGGCGGGCAGGGc                                                                                                                                  >  1:229075/1‑251 (MQ=255)                                                                                                                         cGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCata                                                                                                                           <  2:448762/251‑1 (MQ=255)                                                                                                                            ggAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTCGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGCTAACGGAACGGCGCGGGCACTTGAACCCGGCACCGGGTTAGAGCACGGCGGACCGGGAGTGCatat                                                                                                                          >  2:598920/1‑249 (MQ=255)                                                                                                                            ggAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTGTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGGCCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCatat                                                                                                                          <  1:598920/249‑1 (MQ=255)                                                                                                                               agaaCGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCtt                                                                                                                     <  1:511508/251‑1 (MQ=255)                                                                                                                                  actcAAGATAGTGACTTGTACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATGTCGATGGTTCGCGTCTTAGCTAGCGCAACGCGCCAAAGATCATCTTAGGCAAGCCCTTTAGTGTGATTGTCGTCGGGGATTCGCGTTGTAC‑TGGTGTTGTGTAAGCGCCTGATGTGAACGGAACGTCGCTGACCATTGAACCCCGCACCTGGTTAGAGCAAGGCTGGCTTGTCGTGCATAGCCCTTTAc                                                                                                                  <  2:50620/248‑1 (MQ=255)                                                                                                                                         aTAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGGACCGGGTTAGAGCAAGTCGGGCAGGGCGTGCATATCCCGTTACCCGAAgg                                                                                                           <  2:82772/251‑1 (MQ=255)                                                                                                                                            gggATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGa                                                                                                        <  1:302400/251‑1 (MQ=255)                                                                                                                                                  tataCCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGGAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCGGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGc                                                                                                   <  1:227906/250‑1 (MQ=255)                                                                                                                                                        ggTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTAGATGTTTCACGTCTTAGCGAGCTAAACGCGCCAAATATCATCTTAGGCCAGCCCTTTATTGTGATTAGCGACGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCAGAAAGGGACTGGAGGaacc                                                                                             >  1:290543/1‑248 (MQ=255)                                                                                                                                                                     acagtgagttcacgctatacgCCGATTTCGATGTTTCGCGTCTTACCGAGCTCAACGCCCCAAATATCATATTAGTCACGCCCTTTATTGTGATTAGCGTCGGGGCTGCGCGTGGTAT‑TGGCGTGGTGATAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAACAGACGCTGAAGCTg                                                                               >  2:95029/21‑251 (MQ=255)                                                                                                                                                                               aaCGTTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGccc                                                                     >  1:431291/1‑251 (MQ=255)                                                                                                                                                                                       aaGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTCT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGCACGGCGCTGACCCTTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGAGTGCATATACCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTAAAGCTTAATATGGCCCTGAAtttt                                                             >  2:557990/1‑250 (MQ=255)                                                                                                                                                                                        aGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGCACGGCGCTGACCATTGAACCCGGCACCGGGGTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAGGGGGACTGGCGGAAGCAGAACCTGGAGccggataatgccctgtatata                                                             >  2:110076/1‑230 (MQ=255)                                                                                                                                                                                                 cGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTAcc                                                    <  2:345456/250‑1 (MQ=255)                                                                                                                                                                                                  gATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCACGCCCTTTATTGTGATTAGCGTCAGGGATGCGCGTGGCAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCGGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCCGGGCGTGCATATCCCTTTACCCGAAGGGGCCTGGCGGAAGCAGAAACTGAAGCTGGATATGGCCCTGAATTTAAGCGGTGCCgg                                                  >  1:234160/1‑251 (MQ=255)                                                                                                                                                                                                            ggctgaGCGAGCTCAACGCGGCAAATATCATCTTAGGCAAGCCCTTTAGTGTGATTAGCGTCGGGGATGCGCGTGTTATGTGG‑GTGGTGAAAGCGCCAGAAGTTAACGGAACGGCGCTGACCATTGAGCCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGGATATCCCGTGATCCGAAGGGGTCTGTCGGAGGCAGAGCCTGAAGCTCAATATGGGCCTGAATGTGAGCGGTCCGGGCGCTCTTTCt                                        <  2:538537/246‑1 (MQ=255)                                                                                                                                                                                                              cTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCCTTGAACCCGGCACCGGGGTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGGAGCAGAACCTGAAACTGAATAATGCCCTGAATTTAAAGGGTACACGCGATCCTTCtgg                                      >  2:354458/1‑250 (MQ=255)                                                                                                                                                                                                                            cgcgCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCCTATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGGATATGGCCCTGAATTTTAGCGGTACCGGCGATCTTTCTTTGGGGCCTgccgggg                        >  2:155259/1‑250 (MQ=255)                                                                                                                                                                                                                                    atatCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGc                >  1:236351/1‑251 (MQ=255)                                                                                                                                                                                                                                      atcatcTTAGGCAAGCCCTTTATTGTGAGTAGCGTCGGAGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGTTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGa              <  2:170958/251‑1 (MQ=255)                                                                                                                                                                                                                                           cTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGCCGGGCGTAATAGCGAAATg          >  1:160885/1‑250 (MQ=255)                                                                                                                                                                                                                                                  aaGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTCT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGCCCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATTTCCCTTTACCCGAAGGGGCCTGGCGGCAGCCGAAACTGACGCTGAATGTGGCCCTGAATTTAAGCGGTACCGGGGGTCTTTCTTTGGTGCCTGCCGGGCGGAATAGCGAAATGACCTTAAc  >  2:611817/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       CATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGTTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAT‑TGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAAC  >  NC_000913/4638354‑4638845

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