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.5%	109.6 / NA	40	E81K (GAA→AAA)	uxaB	tagaturonate reductase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (16/23);  minor base G (0/1);  total (16/24)
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

CCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTCCCCCTTCTCATTCAGGCCACGGATAATGG‑TGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAAC‑GACCACGCCAGAATTCAGATCGGTGTGCTC‑ATTCAGGAGATCGATTTGCCAGTCAACAAAGGCGCGCAGGAAGTTACCTTC‑ACCAAACTGAATGATACGTTCTGGATACTGTGCACCGGGAAAATCGCGACGATTTAGTGTTTTC  >  NC_000913/1610197‑1610678                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   ccAGCGCACCGTTGCAATGGCTAAAGCGTTCGAACAGCAGACGGGTAAGGTTTGCCGGATAGCTTACCGCTGGCGCGACATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGGGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGCGCGAATAAGCCGCGCGTCGCTGACCGCTTTCCCCTtc                                                                                                                                                                                                                                             >  2:353389/1‑251 (MQ=255) ccAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCGTACCGCTGGCGCGTCATCGAATGTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATGTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCtt                                                                                                                                                                                                                                              <  2:87342/250‑1 (MQ=255)          cGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAgg                                                                                                                                                                                                                                    <  1:392406/251‑1 (MQ=255)           gTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGc                                                                                                                                                                                                                                   >  2:269856/1‑251 (MQ=255)                aaTGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACgg                                                                                                                                                                                                                              >  1:378816/1‑251 (MQ=255)                    gCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACACGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGACAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGAAGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAAGCCACGGATa                                                                                                                                                                                                                           >  1:449964/1‑250 (MQ=255)                        aagcggtcgagaggcagACGGGTCAGTTTTGCCGGATAGCTCACCTCTGGCGCGTCATAGAATTTATCGCCAGCCTGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACACAGCGCATTTCCGGGTTGTCCGCAAGTTTCAGGACTTCATCGTATTCGCTGAAGCAGCTACTTTCACGATTAACCGAGAGAAAACTACGCCAGTCGCTGCCCGCTCTCCCCTTCTCATACAGGCCACGGAGAAtgg‑t                                                                                                                                                                                                                     >  2:422790/14‑251 (MQ=255)                            gTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCACATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGg‑tggtg                                                                                                                                                                                                                 <  1:386334/251‑1 (MQ=255)                               cGAACAGCAGACGGGTCTGTTTTGCCGGATAGCGTACCGCTGTCGCGTCATCGAATTTATCGCCCGCATGGAAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTGCTCATTCAGGCCACGGATAATGG‑TGGTGTa                                                                                                                                                                                                               <  2:204145/250‑1 (MQ=255)                                       agACGGGGCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGG‑TGGTGTACAGACCATc                                                                                                                                                                                                      <  2:143699/251‑1 (MQ=255)                                                 gTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGG‑TGGTGTACAGACCATCCTGCGTGCTc                                                                                                                                                                                            <  2:378816/251‑1 (MQ=255)                                                     tGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTGCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGCTTTCACGATTAACCGAGCGAATCAGCCGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGCTAATGG‑TGGTGGACAGACCATCCTGCGTGCACAGTg                                                                                                                                                                                        >  2:169745/1‑251 (MQ=255)                                                           aTAGCTGACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGG‑TGGTGTACAGACCATCCTGCGTGCTCAGTGAcggcg                                                                                                                                                                                  <  1:110160/251‑1 (MQ=255)                                                             aGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCGCATTCAGGCCACGGATAATGG‑TGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCggg                                                                                                                                                                                <  2:378455/251‑1 (MQ=255)                                                             aGCTTACCGCTGGCGCGTCATCGAATTTAGCGCCCGCATGGTAGCTGATCCCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTACGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGG‑CGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCggg                                                                                                                                                                                <  1:67864/251‑1 (MQ=255)                                                               cTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGCTAATGG‑TGGTGTACAGACCATCCTGCGTGCTCAGTGCCGGCGGGa                                                                                                                                                                               >  2:236139/1‑250 (MQ=255)                                                                               cATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGCTAATGG‑TGGTGTACAGACCATCCTGCGTGCTCAGTGCCGGCGGGAATGAAGTTTCAATCggg                                                                                                                                                              >  1:411440/1‑250 (MQ=255)                                                                                        taacGCCCGCATGGTAGCTGATACCCGCTTCGGGGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGG‑TGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGaacaac                                                                                                                                                      <  2:377077/247‑1 (MQ=255)                                                                                              ccGCATGGTAGCTGAGACCCGCTGCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGTTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGGAGACGCTGATTTCACGATTAACCGATCGAATCAGACGCGCGTCGCTGACCGCGTTCCCCTTCTCATTCAGGCTACGGATAATGT‑TGGTGTACAGACCATCCTGCGTGCTCAGTGACCGCGGGAATGAAGTTGCAATCGGACGAACAAC‑GACCACg                                                                                                                                              <  2:68492/251‑1 (MQ=255)                                                                                                        gCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGG‑TGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAAC‑GACCACGCCAGAATTCa                                                                                                                                    <  1:241764/251‑1 (MQ=255)                                                                                                                            tgtTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGG‑TGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAAC‑GAACACGCCAACATTCAGATCGGTGTGGTC‑ATTCAgg                                                                                                               >  1:385766/1‑251 (MQ=255)                                                                                                                                 gagaAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGG‑TGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAAC‑GACCACGCCAGAATTCAGATCGGTGTGCTC‑ATTCAGGAGATc                                                                                                          <  1:169745/251‑1 (MQ=255)                                                                                                                                 gagaAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGG‑TGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAAC‑GACCACGCCAGAATTCAGATCGGTGTGCTC‑ATTCAGGAGATc                                                                                                          <  1:269856/251‑1 (MQ=255)                                                                                                                                                        gTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGCCCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGG‑TGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAAC‑GACCACGCCAGAATTCAGATCGGGGTGCTC‑ATTCAGGAGATCGATTTGCCAGTCAACAAAGgcgc                                                                                   >  1:240588/1‑251 (MQ=255)                                                                                                                                                          tgtgCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGG‑TGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAAC‑GACCACGCCAGAATTCAGATCGGTGTGCTC‑ATTCAGGAGATCGATTTGCCAGTCAACAAAGgcgcg                                                                                  <  1:89535/250‑1 (MQ=255)                                                                                                                                                              cgcCAGCTTCAGGAATTCATCGGCTTCACTGTAGGCGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCAGTCAGGCCACGGATAATGG‑TGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAAC‑GACCACGCCAGAATTCAGATCGGTGTGCTC‑ATTCAGGAGATCGATTTGCCAGTCAACAAAGGCGCGCAGGa                                                                             <  2:108620/251‑1 (MQ=255)                                                                                                                                                               gcCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGG‑TGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAAC‑GACCACGCCAGAATTCAAATCGGTGTGCTC‑ATTCCGGAGATCGATTTGCCAGTCAACAAAAGCGCGCAGGa                                                                             >  1:460996/1‑250 (MQ=255)                                                                                                                                                                     ttCAGGAATTCATCGTTGTCACTGTAGACGCTGATTTCACGATTAACCGAGAGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCATGCCACTGATAATGG‑TGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAAC‑GACCACGCCAGAATTCAGATCGGTGTGCTC‑ATTCAGGAGATCGATTTGCCAGTCAACAAAGTCGCGCAGGAAGTTAc                                                                       <  2:415444/250‑1 (MQ=255)                                                                                                                                                                       cAGGAATTCCTCGTATTCACTGTAGCCGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTGCCCCTTCTCATTCAGGCCACGGATAATGG‑TGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAAC‑GACCACGCCAGAATTCAGATCGGTGTGCGC‑ATTCAGGAGATCGATTTGCCAGTCAACAAAGGCGCGCAGGAAGTTACCtt                                                                    <  2:240588/251‑1 (MQ=255)                                                                                                                                                                                aTCGTATTCACTGTAGACGCTGATTTCACGATTAACCAAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGG‑TGGTGTACAGACCATCCTGCGTGCTCAGTGCCGGCGGGGATGCAGGTTTAATCGGCCGAACAAC‑CGCCCCCCCACAATTCAGATCGGTGTGc                                                                                                                         >  2:22818/1‑190 (MQ=255)                                                                                                                                                                                aTCGTATTCACTGTAGACGCTGATTTCACGATTAACCAAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCGTCTCATGCAGGCCACGGATAATGG‑TGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAAC‑GACCACGCCAGAATTCAGATCGGTGTGc                                                                                                                         <  1:22818/190‑1 (MQ=255)                                                                                                                                                                                         actgtatgagccggtttCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGG‑TGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAAC‑GACCACGCCAGAATTCAGATCGGTGTGCTC‑ATTCAGGAGATCGATTTGCCAGTCAACAAAGGCGCGCAGGAAGTTACCTTC‑ACCAAACTGAATGATa                                                  <  2:416599/236‑1 (MQ=255)                                                                                                                                                                                             gagcctcTGATTTCACGATTAACCGAGCTAATCCGACGTGCGGCGCTGACCGCCTTCCCCTCCTCATTCAGGCCACGGATAATGGGTGG‑GTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAAGCGGACTAACAAC‑GACCACGCCAGAATTCAGATCGGTGTGCAC‑AATCAGGAGATCGATTTGCCAGTCAACAAAGGCGCGCAGGAAGTTACCTGC‑CCCAATCTGAATGATACGTTc                                             <  2:449964/245‑1 (MQ=255)                                                                                                                                                                                                    tGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGG‑TGGTGTACAGACCATCCTGCGTGCTCAGTGCCGGCGGGACTGGAGGTTCAATCTGACGAACCCCCGAC‑ACGCCAGAATTCAGATCGGTGTGCTCCATT‑AGGAGATCGCTGTTCCCGTTCACCACGGCGCGCAGGAAGTTACCTTCCACCaacatattgactagttctgtgaac                                       >  2:382223/1‑229 (MQ=255)                                                                                                                                                                                                                 aaCCGAGCGAATCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGG‑TGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAAC‑GACCACGCCAGAATTCAGATCGGTGTGCTC‑ATTCAGGAGATCGATTTGCCAGTCAACAAAGGCGCGCAGGAAGTTACCTTC‑ACCAAACTGAATGATACGTTCTGGATACTGTGCACCGGGaa                         <  1:25647/251‑1 (MQ=255)                                                                                                                                                                                                                     cccgatctCAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGG‑TGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAAC‑GACCACGCCAGAATTCAGATCGGTGTGCTC‑ATTCAGGAGATCGATTTGCCAGTCAACAAAGGCGCGCAGGAAGTTACCTTC‑ACCAAACTGAATGATACGTTCTGGATACTGTGCACCGGGAAAATc                     <  1:103156/244‑1 (MQ=255)                                                                                                                                                                                                                             cAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGG‑TGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGGAGTTTCCATCGGACGAACAAC‑GACCACGCCAGAATTCAGATCGGTGTGCTC‑ATTCAGGCGATCGATTTGCCAGTCAACAAAGGCGCGCCGGACGTTACCTTC‑CCCAAAATGACTGAGTCGTTCTGGTTAATGGGCACCGGGaaaaac                     >  2:103156/1‑241 (MQ=255)                                                                                                                                                                                                                             cAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGG‑TGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGGTTCAATCGGACGAACAAC‑GACCACGCCAGAATTCAGATCGGTGGGCTC‑ATTCAGGAGATCGATTTGCCAGTCAACAAAGGCGCGCAGGAAGTTACCTTC‑ACCAAAATGAATGATACGTTCTGGATAATGTGCACCGGGTAAATCGCGACGAt             >  1:438513/1‑251 (MQ=255)                                                                                                                                                                                                                             cAGACGCGCGTCGCTGACCGCTTTCCCCTTCTCATTCAGGCCACGGATAATGG‑TGGGGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATAGGACGCACAAC‑GACCACGCCAGAATTCAGCTCGGTGGGCTC‑ATTCAGGAGATCGGTTTGCCAGTCAACAAAAGCGCGAAGGAAGTTAACCTC‑ACCAAACTGGATAATAAGTTCTGGGTACTGGGCACCGGGAAAAACGcgacgtt             >  2:263586/1‑249 (MQ=255)                                                                                                                                                                                                                                        cGCTGACCGCTTTCCCCTTCTCATTCAGGCCGCGGATAATGG‑TGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAAC‑GACCACGCCAGAATTCAGATCGGTGTGCTC‑ATTCAGGAGATCGATTTGCCAGTCAACAAAGGCGCGCAGGAAGTTACCTTC‑ACCAAACTGAATGATACGTTCTGGATACTGTGCACCGGGAAAATCGCGACGATTTAGTGTTTTc  <  1:353389/251‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   CCAGCGCACCGTTGAAATGGCTAAAGCGTTCGAACAGCAGACGGGTCAGTTTTGCCGGATAGCTTACCGCTGGCGCGTCATCGAATTTATCGCCCGCATGGTAGCTGATACCCGCTTCGGTGGTGTTAGAGAAAACAAAGCGCATTTCCGGGTTGTGCGCCAGTTTCAGGAATTCATCGTATTCACTGTAGACGCTGATTTCACGATTAACCGAGCGAATCAGACGCGCGTCGCTGACCGCTTCCCCCTTCTCATTCAGGCCACGGATAATGG‑TGGTGTACAGACCATCCTGCGTGCTCAGTGACGGCGGGAATGAAGTTTCAATCGGACGAACAAC‑GACCACGCCAGAATTCAGATCGGTGTGCTC‑ATTCAGGAGATCGATTTGCCAGTCAACAAAGGCGCGCAGGAAGTTACCTTC‑ACCAAACTGAATGATACGTTCTGGATACTGTGCACCGGGAAAATCGCGACGATTTAGTGTTTTC  >  NC_000913/1610197‑1610678

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