Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,610,440	C→T	E81K (GAA→AAA)	uxaB ←	tagaturonate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,610,440	0	C	T	97.1%	184.0 / NA	69	E81K (GAA→AAA)	uxaB	tagaturonate reductase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base T (30/36);  minor base A (1/1);  total (31/38)
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

TTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTCCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCTCATT‑CAGGAGATCG‑ATTTGCCAGTCAACAAAGGCGCGCAGGAAGTTACCTTCACCAAACTGAATGATACGTTCTGGATACTGTGCACCGGGAAAATCGCGACGATTTAGTGTTTTCACAA  >  NC_000913/1610193‑1610682                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     ttATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGCATCAGACGCGCGTCGCTGACCGCTTtccc                                                                                                                                                                                                                                                   >  1:795615/1‑251 (MQ=255)   aTCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACTGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCt                                                                                                                                                                                                                                                 >  1:130548/1‑251 (MQ=255)     ccAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTACCGGAAAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCTTCGCTGACCGCTTACCCCTtc                                                                                                                                                                                                                                               >  1:561406/1‑251 (MQ=255)     ccAGCGCACCGTTGAAAGGTCGAAAGCGTTCGAACGACAGACGGGTCATTTTTGCCGGATAGCTTACCGCTGGTGCGTCATCGAATTGATCGCCCGCATGGTAACTGTTACCCGCTTTGGTGGTGGTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGTGTCGCTGACCGCTTACCCCtt                                                                                                                                                                                                                                                <  2:716081/250‑1 (MQ=255)         cgcACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCCCGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTtctctt                                                                                                                                                                                                                                           >  2:36354/1‑249 (MQ=255)            aCCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTc                                                                                                                                                                                                                                         >  1:56441/1‑250 (MQ=255)            aCCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGATTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCCTTc                                                                                                                                                                                                                                         >  1:593342/1‑250 (MQ=255)                   aaaTGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACg                                                                                                                                                                                                                                 >  1:613332/1‑251 (MQ=255)                            aaGCGATCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGAGTCATCGAATTTATCGCCCGCTTGGTAGCTGACACCCGCTTCGGTGGAGTTAGAGAAAACAAAGCGCATTTCCGGGTGTTGCGCCAGTTTCAGGAATTAATCGATTTCACTGTAGACGCTGATTTCACGATTAACAGATCGAATCAAACGCGCGTCACTGCCCGCTTTCCCCTTCTCATTCAGGCCCCGGATAAtggt                                                                                                                                                                                                                        >  2:548306/1‑251 (MQ=255)                             aGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAAtggtg                                                                                                                                                                                                                       <  1:3448/251‑1 (MQ=255)                                               gggTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCAAGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTATCCCCTTCTCATTCAGGCAACGGATAATGGTGGTGTACAGAACATCCTGc                                                                                                                                                                                                     >  1:443453/1‑251 (MQ=255)                                                   cAGTTTTGCCTGATAGCTTACCGCGGGCGCGTCATCGAATTTATCGCCCGCATGTTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTg                                                                                                                                                                                                  <  2:172585/250‑1 (MQ=255)                                                    aGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCt                                                                                                                                                                                                <  1:386634/251‑1 (MQ=255)                                                                    ttACCGCTGGCGCGTCATCGAATTTATCGCCCGCCTGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGATGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGGGTACAGCACATCCTGCGTGCTCAGTGACGGCgggga                                                                                                                                                                                 >  2:1813/1‑248 (MQ=255)                                                                    ttACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCATTGACGGCGGGa                                                                                                                                                                                  >  2:350403/1‑249 (MQ=255)                                                                                    aTCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGAc                                                                                                                                                                <  1:534938/251‑1 (MQ=255)                                                                                        aaTTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGaa                                                                                                                                                             >  1:438111/1‑250 (MQ=255)                                                                                         aTTTATCGACCGCATGGTAGCTGCTACCCGCTTCGGTGGTGTTAGAGAAACCAAAGCGCATTTCCGGGTTGTGCTCCAGTCTCCGGAATTACTAGTATTCACTGTAGACGCTAATTTCACGATTAACCGAGCGAATCCGACGCGCGTCGCTGACCGATTTACCCTTCTCTTTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGAGTGCTCAGAGGCGGCGGGGAAGTAGTTTCactaggaggtaa                                                                                                                                                            >  2:545086/1‑239 (MQ=255)                                                                                             aTCGCCCGCATGGTAGCTGATCACCTCGTCGGTGGTGTTAGAGTAAACAATGCGCATTTCCGGGTTGTGTGCCAGGTTCAGGAATTCATCGTATTCCCTGTAGACTCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCTTTCAGGTCACTGATAATGGTTGTTTACAGACCATCCTGCGTGCACAGTGACGGCGGTAATGAAGTTTCAATCGGACGACCAACGa                                                                                                                                                       <  2:600527/251‑1 (MQ=255)                                                                                                    ccATGGTAGCTGATACCCTCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGTGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAACGACCACGcc                                                                                                                                                <  2:100559/250‑1 (MQ=255)                                                                                                        ggTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGGTTCAATCGGGCGAACACCGACCACGCCAGa                                                                                                                                             >  1:257345/1‑250 (MQ=255)                                                                                                         gTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAACGACCACGCCAGAAt                                                                                                                                           <  1:183588/251‑1 (MQ=255)                                                                                                             cggATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCCATCGGACGCACACCGGCCCCGCCAGacatca                                                                                                                                        >  2:739849/3‑245 (MQ=255)                                                                                                             cTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGCTAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAACGACCACGCCAGAATTc                                                                                                                                         >  1:697041/1‑249 (MQ=255)                                                                                                                aTACCCGCGTGGGTTGTGTTAGATAAAACAAAGCGCAGTGCCGGGTTGGGCGCCAGTTTCAGGAAGTCATCGTAGTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCTCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGAGAATGGTGGTGTACAGACCATCCTTCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAACGACCACGCCGGAATTCAGATc                                                                                                                                    <  2:279357/251‑1 (MQ=255)                                                                                                                 taccctcatCGGTGGTGTTAGAAAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGCCGTTGATTTCACGATTAACCGAGCGAATCCGATGCGCGTCGCTGTCTGCTTTCCCCTTCTCATTCCGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGTACTGAAGTTTCAATCGGACGAACAACGACCACGCCAGAAGTCAGATc                                                                                                                                    <  2:693191/242‑1 (MQ=255)                                                                                                                   cccGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGCTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTCCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAACGACCACGCCAGAATTCAGATCGgt                                                                                                                                 <  1:150245/251‑1 (MQ=255)                                                                                                                    ccGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAACGACCACGCCAGAATTCAGATCGgtg                                                                                                                                <  1:14212/251‑1 (MQ=255)                                                                                                                            ggtgtattAGAGAAAACAAAGCGCATTTCAGGGTTGTGAGCCGGTTTAAGGAATTCATCTTATTCCATGTAGACGCTGATTCCACGATTAACCGAGCGAATCAGACGTGCGTCGCTGACCGCTTTCCCCATCTCACTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCCCAGTGACGGCGGGAATGAAGTTTAAATCGGACGAACAACGTCCACGCCAGAATTCAGATCGTTGTGCTCAtg                                                                                                                        <  2:737834/245‑2 (MQ=255)                                                                                                                               gtgtTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCTCATT‑CAg                                                                                                                    <  1:739849/251‑1 (MQ=255)                                                                                                                                    agagAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCGTCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCTCATT‑CAGgatat                                                                                                               <  2:839110/251‑4 (MQ=255)                                                                                                                                       gaAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCTCGATTAACCGAGCGAATCAGACGCGCGTCGCTGTCCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCTCATT‑CAGGAGATCG‑a                                                                                                           <  1:554857/251‑1 (MQ=255)                                                                                                                                        aaaaCAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACACCGACCCCGCCAGAATTCAGATCGGTGTGCTCATT‑CAGGGGATCG‑At                                                                                                          >  2:407945/1‑251 (MQ=255)                                                                                                                                             aaaGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCTCATT‑CAGGAGATCG‑ATTTGcc                                                                                                     >  1:176921/1‑251 (MQ=255)                                                                                                                                                   cATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCTCATT‑CAGGAGATCG‑ATTTGCCAGTca                                                                                                >  1:292606/1‑250 (MQ=255)                                                                                                                                                      ttCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCTCATT‑CAGGAGATCG‑ATTTGCCAGTcaaca                                                                                             <  2:305365/250‑1 (MQ=255)                                                                                                                                                      ttCCGAGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGAACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAACGCCCACACCAAAATTAAGATAGGTGAGCTCATT‑AAgggaaacgatttggccgtccaca                                                                                              >  2:823334/1‑229 (MQ=255)                                                                                                                                                         cGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGGATGAAGTTTCAATCGGACGAAAAACGACCACGCCAGAATTCAGATCGGTGTGCTCATT‑CAGGAGATCG‑ATTTGCCCGTCAACAAAgg                                                                                         >  1:565721/1‑251 (MQ=255)                                                                                                                                                         cGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGATTCAATAGGGAGAAACACAACCACGCCAGAATTGCGATCGGCGGGCTCATT‑CAGGAGAACG‑GTTTGCCAGTCAACAAAgg                                                                                         >  2:822153/1‑251 (MQ=255)                                                                                                                                                          gggTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTGCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCTCATT‑CAGGAGATCG‑ATTTGCCAGTCAACAAAGgc                                                                                        <  2:292606/251‑1 (MQ=255)                                                                                                                                                                 gcgcCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCTCATT‑CAGGAGATCG‑ATTTGCCAGTCAACAAAGGCGCGCAg                                                                                  <  1:36354/250‑1 (MQ=255)                                                                                                                                                                  cgcCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAACGa                                                                                                                                                       <  1:694488/182‑1 (MQ=255)                                                                                                                                                                         ttCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCTCATT‑CAGGAGATCG‑ATTTGCCAGTCAACAAAGGCGCGCAGGAAGTTAcc                                                                         >  1:197714/1‑251 (MQ=255)                                                                                                                                                                               aaTTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCATACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGTCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAATGACGGCGGGAATGAAGTTTCAATCGGACGAACAACTACCACGCCCGAATTCAGATCGGTGTGCTTATTTCTGAAG‑TCG‑ATTTGCCAATCAACAAAGGCGCGCAGGAAATTACCTTCAcc                                                                   >  1:522441/1‑251 (MQ=255)                                                                                                                                                                                      cGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCTCATT‑CAGGAGATCG‑ATTTGCCAGTCAACAAAGGCGCGCAGGAAGTTACCTTCACCAAACTg                                                             <  1:407945/250‑1 (MQ=255)                                                                                                                                                                                       gTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGCCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTACGTGCTCAGTGACGGCGGGGATGAAGTTTCACTCGATCCAACACCGACCACACCCGAATTCCGATCGGCGTGGTCATT‑AATGCCATCG‑ATTTGCCAGGCACCAAAGCCGCGCAGGAAATTCGCTTCACCAAATTGa                                                            >  2:234771/1‑250 (MQ=255)                                                                                                                                                                                            cACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCTCATT‑CAGGAGATCG‑ATTTGCCAGTCAACAAAGGCGCGCAGGAAGTTACCTTCACCAAACTGAATGAt                                                       <  2:176921/250‑1 (MQ=255)                                                                                                                                                                                             aCTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCTCATT‑CAGGAGATCG‑ATTTGCCAGTCAAAAAAGGCGCGCAGGAAGTTACCTTCACCAAACTGAATGaaac                                                     >  1:136858/1‑248 (MQ=255)                                                                                                                                                                                             aCTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCTCATT‑CAGGAGATAG‑ATTTGCCAGTCAACAAAGGCGCGCAGGAAGTTACCTTCACCAAACTGAATGATAc                                                     <  1:350403/251‑1 (MQ=255)                                                                                                                                                                                                  aGACGCTAATTTCCCGATTAACCGAGCGAATCTGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCTCTGATAATGGTTTTTTACATACCATCCTGCGTGCTCAGTGACGGCGGGAATGAATTTTCAAGCGGACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCTCATT‑CAGTAGATCG‑ATTTGCCAGTCAACAAAGGCGCGCAGGAAGTTACCTTCACCAAACTGAATGATACGTTCt                                                <  1:521290/251‑1 (MQ=255)                                                                                                                                                                                                       cTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAACGACCCCGCCAGAATTCAGATCGGGGCGCTCATT‑CAGGAGATCG‑ATTTGCCCGTCAACAAAGGCGCGCAGGGAGTTACCTTCACCAAACTGAATGATACGTTCTGGAt                                            >  1:732690/1‑250 (MQ=255)                                                                                                                                                                                                         gATTTCACGATGAACCGAGCGAAGCATCCGCTCGCTGATGACCGCTGTCTCCTTCTCATTCAGGTCACGTATAATGTTGGTGAACAGATCATCCTCCTTGCTCAGTGACGGTAGGAATGAAGTTTCAATCGGACGAACAACGACCACGCCAGAATGCAGATCGGTGTGCACATT‑CAGGAGAACG‑ATTTGCCAGTCAACAAGGGCGCGCAGGAAGGTACCTTCATCTAACTGAATGATACGTTTTGGAt                                            <  2:778976/248‑1 (MQ=255)                                                                                                                                                                                                               aCGATTAAGCGAGCGAATCTGACGCGCGTCTCTGACCGATTTCACCTTCGCAGTCAGGTCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGAGAATGAAGTTTCAATCGGACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCTCAAT‑CAGGAGATCG‑ATTTGCCAGTCAACAAAGGCGCGCAGGAAGTTCCCTTCACCAAACTGAATGATACGTTCTGGATACTGGGCAc                                   <  2:178003/251‑1 (MQ=255)                                                                                                                                                                                                                 gATTAACCGAGCGAATCAGACGAGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCACTCGGACGAACAACGACCACCCCCGAAATCAAATCGGTGTGCTCATT‑CCGGAGATCG‑GTTTGCCCGTCAACAAAAGAGCGCAGGGAGTTACCTTCACCACACTGGATGATACATTCTGGGTACTTTGCCCCg                                 >  2:312804/1‑251 (MQ=255)                                                                                                                                                                                                                  aTTAACCGAGCGAACCAGACGCCCGGCGCGGACCTCTTTCACCTTCTCATTCAGTCCACTGGTAAGGGTGGATTCCAGACCATCCTGCGTGCTCAGTTACGGCGGGAATGAAGTGTCAATCGTACGAACTACGGCCGCGCCAGAATTCAGATCGGTGTGCTCATT‑CAGGAGATCG‑AGTTGTGAGTCAACTAAGTCGCACAGGGAGTGATCTACTTCAAACTTAATGATTCCTTGTGGATATTGTGCATCgg                                <  2:832897/251‑1 (MQ=255)                                                                                                                                                                                                                     acccGAGCGAACCCGACGAGCGGTGCTGACCGCTTTCCCCTTCTCAGTCAGGGCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGATGGCGGGAATGAAGTTTCAATCGGACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCTCAGT‑CAGGAGATCG‑ATTTGCCAGTCAACAAAGGCGCGCAGGAAGTTACCTGCACCAAACTGAATGAGACGTTCTGGTTACTGTGCACCGGGaa                             <  2:778721/249‑1 (MQ=255)                                                                                                                                                                                                                     aaCCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAACGACCACGCCAGAATTCAGATCGGGGGGCTCATT‑CAGGAGATCG‑ATTTGCCAGTCAACAAAGGCGCGCCGGCAGTTACCTTCACCAAACTGAATGATACGGTCTGGGTACTGTGCACCGGGa                              >  1:68246/1‑250 (MQ=255)                                                                                                                                                                                                                     aaCCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGAGAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCTCATT‑CAGGAGATCG‑ATTTGCCAGTCAACAAAGGCGCGCAGTAAGTTACCTTCACCAAACTGAATGATACGTTCTGGATACTGTGCACCGGGaa                             <  2:68246/251‑1 (MQ=255)                                                                                                                                                                                                                      aCCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAACGACCACGc                                                                                                                                                 <  2:156274/136‑1 (MQ=255)                                                                                                                                                                                                                      aCCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAACGACCACGc                                                                                                                                                 >  1:156274/1‑136 (MQ=255)                                                                                                                                                                                                                       ccGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGCATGAAGGTTCCATCGGACGGACAAAGACCACCCCAGAATTCTGATCGGGGGTCGCAAT‑CAGGAGATCG‑GTTTGCCAGTCAACAAAAGCGCGCAGGAAGTTACATTCACCAAACTGAAAGATCCGTTTTGGAAACTGAGTACCGGGaaa                            >  1:535480/1‑250 (MQ=255)                                                                                                                                                                                                                         gAGCGAATCAGACGCGCGTCGCGGACCGCTTTCCCCTTCGCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCTCATT‑CAGGAGATCG‑ATTTGCCAGTCAACAAAGGCGCGCAGGAAGTTACCTTCACCAAACTGAATGATACGTTCTGGATACTGTGCACCGGGAAAATc                         <  1:673659/251‑1 (MQ=255)                                                                                                                                                                                                                              cagcAGACGCTCGTCGTCGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAACGACCACGCCAGAATTCCGATCGGTGTTCTCATT‑CAGGAGATCG‑ATTTGCCAGTCAACAAAGGCGCGCAGGAAGTTACCTTCACCAACCTGAATGATACGTTATGGATACTGTGCACCGGGAAAATCGcgac                    <  2:156688/248‑1 (MQ=255)                                                                                                                                                                                                                                tCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAAAGACCACGCCAGAATTCAGCTTGGTGTGCTCATT‑CAAGAGCACGCATTT‑CCAGTCAACAAAAGAGCGGAGGAAGTTACCTTCACCAAAATGAATGATAAGTTTTGGATACTGTGCACCGGGAAAATCGcgaagc                  >  2:395325/1‑248 (MQ=255)                                                                                                                                                                                                                                        ggggctgtGACCGCTTTCCCCTTCCCATTCAGGCCCCGGATAATGGTGGTGTACAGACCAACCTGAGGGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAACGACCACGCCAGAATTCAGAGCGGTGTGCTCATT‑CAGGAGATCG‑ATTTGCCAGTCAACAAAGGCGCGCAGGAAGTTACCTTCACCAAACTGAATGATACGTTCTGGATACTGTGCACCGGGAAAATCGCGACGATTTAgtgt          <  2:697041/244‑1 (MQ=255)                                                                                                                                                                                                                                          gTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATTAAGTTTCAATCGGACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCTCATT‑CAGGAGATCG‑ATTTGCCAGTCAACAAAGGCGCGCAGGAAGTTGCCTTCACCAAACTGAATGATACGTTCTGGATACTGTGCACCGGGAAAATCGCTACGATTTAGTGttt        <  2:613332/251‑1 (MQ=255)                                                                                                                                                                                                                                            cGCTGACCGCTTTCCCCTTCTCAGCCAGGCCACGGATAATGGTAGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAACGTCCACGCCAGAATTCAGATCGGTGTGCTCATT‑CAGGAGATCG‑ATTTGCCAGTCAACAAAGGCGCTCAGGAAGTTACCTTCACCAAACTGAATGATACGTTCTGGAGACTGTGCACCGGGAAAATCGCGACGATT‑AGTGTTTTc      <  1:408493/250‑1 (MQ=255)                                                                                                                                                                                                                                             gCTGACCGCTTTCCCCTCCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCGCATT‑CAGGAGATCG‑ATTTGCCAGTCAACAAAGGCGCGCAGGAAGTTACCTTCACCAAACTGAATGATACGTTCTGGATACTGTGCACCGGGAAAATCGCGACTATTTAGTGTTTTca     <  2:565721/251‑1 (MQ=255)                                                                                                                                                                                                                                                gACCTCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCTCATT‑CAGGAGATCG‑ATTTGCCAGTCAACAAAGGCGCGCAGGAAGTTACCTTCACCAAACTGAATGAGACGTTCTGGATACTGTGCACCGGGAAAATCGCGACGATTTAGTGTTTTcac    <  1:402695/249‑1 (MQ=255)                                                                                                                                                                                                                                                gACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCTCATT‑CAGGAGATCG‑ATTTGCCAGTCAACAAAGGCGCGCAGGAAGTTACCTTCACCAAACTGAATGATACGTTCTGGATACTGTGCACCGGGAAAATCGCGACGATTTAGTGTTTTCACaa  <  1:143383/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     TTATCCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTCCCCCTTCTCATTCAGGCCACGGATAATGGTGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCTCATT‑CAGGAGATCG‑ATTTGCCAGTCAACAAAGGCGCGCAGGAAGTTACCTTCACCAAACTGAATGATACGTTCTGGATACTGTGCACCGGGAAAATCGCGACGATTTAGTGTTTTCACAA  >  NC_000913/1610193‑1610682

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