Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,036,834	G→A	C29Y (TGT→TAT)	hyaF →	protein HyaF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,036,834	0	G	A	97.9%	149.9 / NA	48	C29Y (TGT→TAT)	hyaF	protein HyaF
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (28/19);  minor base C (1/0);  total (29/19)
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

AAGCCATCGGCGATCGCTTTGGCGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAACTATCGCGGCGTGCTGAATGGTATTCACCCGTGGGCGGAACTGATAAACCTGATGCGCGGGCTTGTCGAACCGCAGCAGGAGCGTGCCTCATGAGCGAAACTTTTTTCCATCTGCTGGGGCCAGGAACGCAACCG‑AACGATGACAGTTTCAGCATGAATCCACTGCCGATCACCTGTCAGGTGAATGATGAACC‑GAGTATGG‑CGGCCCTGG‑AGCAATGTGC‑TCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAACATCAACTAAGCGAACGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCTTCTCGGGGAAGGGGAAGTGTCAGTGCGCATTCAGCAGGCTGACGACAGTGAAAGTGAAATAC  >  NC_000913/1036595‑1037076                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                       aaGCCATCGGCGATCGCTTTGGCGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAACTATCGCGGCGTGCTGAATGGTATTCACCCGTGGGCGGAACTGATAAACCTGATGCGCGGGCTTGTCGAACCGCAGCAGGAGCGTGCCTCATGAGCGAAACTTTTTTCCATCTTCTGGGGCCAGGAACGCAACCG‑AACGATGACAGTTTCAGCATGAATCCACTGCCGATCACCTATCAGGTGAatg                                                                                                                                                                                                                                             >  2:109630/1‑251 (MQ=255)      aTCGGCGATCGCTTTGGCGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAACTATCGCGGCGTGCTGAATGGTATTCACCCGTGGGCGGAACTGATAAACCTGATGCGCGGGCTTGTCGAACCGCAGCAGGAGCGTGCCTCATGAGCGAAACTTTTTTCCATCTGCTGGGGCCAGGAACGCAACCG‑AACGATGACAGTTTCAGCATGAATCCACTGCCGATCACCTATCAGGTGAATGATGaa                                                                                                                                                                                                                                        <  1:428501/251‑1 (MQ=255)      aTCGGCGATCGCTTTGGCGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAACTATCGCGGCGTGCTGAATGGTATTCACCCGTGGGCGGAACTGATAAACCTAATGCGCGGGCTTGTCGAACCGCAGCAGGAGCGTGCCTCATGAGCGAAACTTTTTTCCATCTGCTGGGGCCAGGAACGCAACCG‑AAAGAAGACAGTTTCAGCATGAATCCACTGCCGATCACCTAACAGGAGAAtgatga                                                                                                                                                                                                                                         >  2:291662/1‑250 (MQ=255)         ggCGATCGCTTTGGCGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAACTATCGCGGCGTGCTGAATGGTATTCACCCGTGGGCGGAACTGATAAACCTGATGCGCGGGCTTGTCGAACCGCAGCAGGAGCGTGCCTCATGAGCGAAACTTTTTTCCATCTGCTGGGGCCAGGAACGCAACCG‑AACGATGACAGTTTAAGCATGAATCCACTGCCGATCACCTATCAGGGGCATGATGAACC‑g                                                                                                                                                                                                                                    >  2:326958/1‑251 (MQ=255)            gATCGCTTTGGCGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAACTATCGCGGCGTGCTGAATGGTATTCACCCGTGGGCGGAACTGATAAACCTGATGCGCGGGCTTGTCGAACCGCAGCAGGAGCGTGCCTCATGAGCGAAACTTTTTTCCATCTGCTGGGGCCAGGAACGCAACCG‑AACGATGACAGTTTCAGCATGAATCCACTGCCGATCACCTATCAGGTGAATGATGAACC‑GAGt                                                                                                                                                                                                                                 >  2:197886/1‑251 (MQ=255)               cGCTTTGGCTTTTTTGTCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAACTATCGCGGCGTGCTGAATGGTATTCACCCGTGGGCGGAACTGATAAACCTGATGCGCGGGCTTGTCGAACCGCAGCAGGAGCGTGCCTCATGAGCGAAACTTTTTTCCATCTGCTGGGGCCAGGAACGCAACCG‑AACGATGACAGTTTCAGCATGAATCCACTGCCGATCACCTATCAGGTGAATGATGAACC‑TAGTAt                                                                                                                                                                                                                               <  2:258517/250‑1 (MQ=255)                   tgggCGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAACTATCGCGGCGTGCTGAATGGTATTCACCCGTGGGCGGAACTGATAAACCTGATGCGCGGGCTTGTCGAACCGCAGCAGGAGCGTGCCTCATGAGCGAAACTTTTTTCCATCTGCTGGGGCCAGGAACGCAACCG‑AACGATGACAGTTTCAGCATGAATCCACTGCCGATCACCTATCAGGTGAATGATGAACC‑GAGTATgg‑cg                                                                                                                                                                                                                          <  2:347447/248‑1 (MQ=255)                          cgtgcgcgttCCTGCCAATTTAGTGTTTACCGGCGGAAACTATCGCGGCGTGCTGAATGGTATTCACCCGTGGGCGGAACTGATAAACCTGATGCGCGGGCTTGTCGAACCGCAGCAGGAGCGTGCCTCATGAGCGAAACTTTTTTCCATCTGCTGGGGCCAGGAACGCAACCG‑AACGATGACAGTTTCAGCATGAAGCCACTGCCGATCACCTATCAGGTGAATGATGAGCC‑GAGTATGG‑CGGCCCTgg                                                                                                                                                                                                                   <  2:369184/242‑1 (MQ=255)                                 tttCCTGCCACTTTGGTGTTTACCGGCGGAAACTATCGCGGCGTGCTGAATGTTTTTCACCCGTGGGCGGAACTGATAAACCTGATGCGCGTGCTTGTCGAACCGCCGCAGGAGCGTGCCTCATGAGCGAAACTTTTTCCCATCTGCTGGGGCCAGGAACGCAACCG‑AACGATGACAGTTTCAGCATGAATCTACTGCCGATCACCTATCAGGTGAATGATGAACC‑GAGTATGG‑CGGCCCTGG‑AGCAAtg                                                                                                                                                                                                           <  1:419740/250‑1 (MQ=255)                                              tAGTGTTTACCGGCGGAAACTATCGCGGCGTGCTGAATGGTATTCACCCGTGGGCGGAACTGATAAACCTGATGCGCGGGCTTGTCGAACCGCAGCAGGAGCGTGCCTCATGAGCGAAACTTTTTTCCATCTGCTGGGGCCAGGAACGCAACCG‑AACGATGACAGTTTCAGCATGAATCCACTGCCGATCACCTATCAGGTGAATGATGAACC‑GAGTATGG‑CGGCCCTGG‑AGCAATGTGC‑TCACAGCCCGc                                                                                                                                                                                            >  1:196600/1‑251 (MQ=255)                                                gtgtTTACCGGCGGAAACTATCGCGGCGTGCTGAATGGTATTCACCCGTGGGCTGAACTGATAAACCTGATGCGCGGTCTTGTCGAACCGCAGCAGGAGCGTGCCTCATGAGCGAAACTTTTTTCCATCTGCTGGGGCCAGGAACGCAACCG‑AACGATGACAGTTTCAGCATGAATCCACTGCCGATCACCTATCAGGTGAATGATGAACC‑GAGTATGG‑CGGCCCTGG‑AGCAATGTGC‑TCACAGCCCGCAg                                                                                                                                                                                          <  1:283124/251‑1 (MQ=255)                                                gtgnTTACCGGCGGAAACTATCGCGGCGTGCTGAATGGTATTCACCCGTGGGCGGAACTGATAAACCTGATGCGCGGGCTTGTCGAACCGCATCAGGAGCGTGCCTCATGAGCGAAACTTTTTTCCATCTGCTGGGGCCAGGAACGCAACCG‑AACGATGACAGTTGCAGCATGAATCCACTGCCGATCACCTATCAGGTGAATGATGAACC‑GAGTATGG‑CGGCCCTGG‑AGCAATGTGC‑TCACAGCCCGCAg                                                                                                                                                                                          <  2:419806/251‑1 (MQ=255)                                                 tgtTTACCGGCGGAAACTATCGCGGCGTGCTGAATGGTATTAACCCGTGGGCGGAACTGATAAACCTGATGCGCGGGCTTGTCGAACCGCAGCAGGAGCGTGCCTCATGAGCGAAACTTTTTTCCATCTGCTGGGGCCAGGAACGCAACCG‑AACGATGACAGTTTCAGCATGAATCCACTGCCGATCACCTCTCAGGTGAATGATGAACC‑GAGTATGG‑CGGCCCTGG‑AGAAAtgtggtccccacccgcagg                                                                                                                                                                                          >  2:67548/1‑236 (MQ=255)                                                               aaCTATCGCGGCGTGCTGAATGGTATTCACCCGTGGGCGGAACTGATAAACCTGATGCGCGGGCTTGTCGAACCGCAGCAGGAGCGTGCCTCATGAGCGAAACTTTTTTCCATCTGCTGGGGCCAGGAACGCAACCG‑AACGATGACAGTTTCAGCATGAATCCACTGCCGATCACCTATCAGGTGAATGATGAACC‑GAGTATGG‑CGGCCCTGG‑AGCAATGTGC‑TCACAGCCCGCAGGTGATTGCGCTGtt                                                                                                                                                                            >  2:209374/1‑250 (MQ=255)                                                               aaCTATCGCGGCGTGCTGAATGGTATTCACCCGTGGGCGGAACTGATAAACCTGATGCGCGGGCTTGTCGAACCGCAGCAGGAGCGTGCCTCATGAGCGAAACTTTTTTCCATCTGCTGGGGCCAGGAACGCAACCG‑AACGATGACAGTTTCAGCATGAATCCACTGCCGATCACCTATCAGGTGAATGATGAACC‑GAGTATGG‑CGGCCCTGG‑AGCAATGTGC‑TCACAGCCCGCAGGTGATTGCGCTGTTa                                                                                                                                                                           <  2:103754/251‑1 (MQ=255)                                                                aCTATCGCGGCGTGCTGAATGGTATTCACCCGTGGGCGGAACTGATAAACCTGATGCGCGGGCTTGTCGAACCGCAGCAGGAGCGTGCCTCATGAGCGAAACTTTTTTCCATCTGCTGGGGCCAGGAACGCAACCG‑AACGATGACAGTTTCAGCATGAATCCACTGCCGATCACCTATCAGGTGAATGATGAACC‑GAGTATGG‑CGGCCCTGG‑AGCAATGTGC‑TCACAGCCCGCAGGTGATTGCGCTGTTaa                                                                                                                                                                          <  1:219329/251‑1 (MQ=255)                                                                     cgcgGCGTGCTGAATGGTATTCACCCGTGGGCGGAACTGATAAACCTGATGCGCGGGCTTGTCGAACCGCAGCAGGAGCGTGCCTCATGAGCGAAACTTTTTTCCATCTGCTGGGGCCAGGGACGCAACCG‑AACGATGACAGTTTCAGCATGAATCCACTGCCGATCACCTATCAGGTGAATGCTGAACC‑GCGTATGG‑CGGCCCTGG‑AGCAATGTGC‑TCACAGCCCGCAGGTGATTGCGCTGTTAAACg                                                                                                                                                                       >  2:104514/1‑249 (MQ=255)                                                                      gcgncgTGCTGAATGGTATTCACCCGTGGGCGGAACTGATAAACCTGATGCGCGGGCTTGTCGAACCGCAGCAGGAGCGTGCCTCATGAGCGAAACTTTTTTCCATCTGCTGGGGCCAGGAACGCAACCG‑AACGATGACAGTTTCAGCATGAATCCACTGCCGATCACCTATCAGGTGAATGATGAACC‑GAGTATGG‑CGGCCCTGG‑AGCAATGTGC‑TCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGt                                                                                                                                                                    <  2:418907/251‑1 (MQ=255)                                                                       cggcgTGCTGAATGGTATTCTCCCGTGGGCGGAACTGATAAACCTGATGCGCGGGCTTGTCGAACCGCAGCAGGAGCGTGCCTCATGAGCGAAACTTTTTTCCATCTGCTGGGGCCAGGAACGCAACCG‑AACGATGACAGTTTCAGCATGAATCCACTGCCGATCTCCTATCAGGTTAATGATGAACC‑GAGTATGG‑CGGCCCTGG‑AGCAATGTGC‑TCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGtt                                                                                                                                                                   <  2:371338/251‑1 (MQ=255)                                                                       cggagagatggaTGGAATTCACCCGTGGGCGGAACTGATAAATGTCATGCGCGGGTTTGTCGAACCGCAGCAGGAGCGTGCCTAATGAGCGAAACTTTATTCCCTCTGCTGGGGCCAGGAACGCAACCG‑AACGATGACAGTATCAGCATGAATCCACTCCCGATCACCTATCAGGGGAATGATGAACC‑GAGTATGA‑CGGCCCTGG‑AGCCCTGTGA‑TACCAGCCCACAGGAAATAGCGCTGTTAAAGGAggt                                                                                                                                                                   >  2:294664/12‑249 (MQ=255)                                                                                gAATGGTATTCACCCGTGGGCGGAACTGATAAACCTGATGCGCGGGCTTGTCGAACCGCAGCAGGAGCGTGCCTCATGAGCGAAACTTTTTTCCATCTGCTGGGGCCAGGAACGCAACCG‑AACGATGACAGTTTCAGCATGAATCCACTGCCGATCACCTATCAGGTGAATGATGAACC‑GAGTATGG‑CGGCCCTGG‑AGCAATGTGC‑TCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAACATc                                                                                                                                                           >  1:236602/1‑250 (MQ=255)                                                                                gAATGGTATTCACCCGTGGGCGGAACTGATAAACCTGATGCGCGGGCTTGTCGAACCGCAGCAGGAGCGTGCCTCATGAGCGAAACTTTTTTCCATCTGCTGGGGACAGGAACGCAACCG‑AACGATGACAGTTTCAGCATGAATCCACTGCCGATCACCTATCAGGTGAATGATGAACC‑GAGTATGG‑CGGCCCTGG‑AGCAATGTGC‑TCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAACAt                                                                                                                                                            <  1:409914/249‑1 (MQ=255)                                                                                  aTGGTATTCACCCGTGGGCGGAACTGATAAACCTGATGCGCGGGCTTGTCGAACCGCAGCAGGAGCGTGCCTCATGAGCGAAACTTTTTTCCATCTGCTGGGGCCAGGAACGCAACCG‑AACGATGACAGTTTCAGCATGAATCCACTGCCGATCACCTATCAGGTGAATGATGAACC‑GAGTATGG‑CGGCCCTGG‑AGCAATGTGC‑TCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAACATCAAc                                                                                                                                                        <  2:246061/251‑1 (MQ=255)                                                                                    ggTATTCACCCGTGGGCGGAACTGATAAACCTGATGCGCGGGCTTGTCGAACCGCAGCAGGAGCGTGCCTCATGAGCGAAACTTTTTTCCATCTGCTGGGGCCAGGAACGCAACCG‑AACGATGACAGTTTCAGCATGAATCCACTGCCGATCACCTATCAGGTGAATGATGAACC‑GAGTATGG‑CGGCCCTGG‑AGCAATGTGC‑TCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAACATCAACt                                                                                                                                                       <  2:29085/250‑1 (MQ=255)                                                                                       aTTCACCCGTGGGCGGAACTGATAAACCTGATGCGCGGGCTTGTCGAACCGCAGCAGGAGCGTGCCTCATGAGCGAAACTTTTTTCCATCTGCTGGGGCCAGGAACGCCACCG‑AACGATGCCAGTTTCAGCATGAATCCACTGCCGATCACCTATCCGGTGAATGATGAACC‑GAGTATGG‑CGGCCCTGG‑CGCAATGTGC‑TCACATCCCGCAGGTGATTGGGCTGTTAAACGAGTTACAACATCAACTAAGc                                                                                                                                                   >  2:99142/1‑251 (MQ=255)                                                                                         tCACCCGTGGGCGGAACTGATAAACCTGATGCGCGGGCTTGTCGAACCGCAGCAGGAGCGTGCCTCATGAGCGAAACTTTTTTCCATCTGCTGGGGCCAGGAACGCAACCG‑AACGATGACAGTTTCAGCATGAATCCACTGCCGATCACCTATCAGGTGAATGATGAACC‑GAGTATGG‑CGGCCCTGG‑AGCAATGTGC‑TCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAACATCAACTAAGCg                                                                                                                                                  >  1:454489/1‑250 (MQ=255)                                                                                                  ggCGGAACTGATAAACCTGATGCGCGGGCTTGTCGAACCGCAGCAGGAGCGTGCCTCATGAGCGAAACTTTTTTCCATCTGCTGGGGCCAGGAACGCAACCG‑AACGATGACAGTTTCAGCATGAATCCACTGCCGATCACCTATCAGGTGAATGATGAACC‑GAGTATGG‑CGGCCCTGG‑AGCAATGTGC‑TCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAACATCAACTAAGCGAACGCCAAcc                                                                                                                                        <  1:411424/251‑1 (MQ=255)                                                                                                          tGATAAAACTGATGCGCGGGCTTGTCGAACCGCAGCAGGAGCGTGCATCATGAGCGAAACTTTTTTCCATCTGCTGGGGCCAGGAACGCAACCG‑AACGATGACAGTTTCAGCATGAATCAACTGCCGATCACCTATCAGGTGAATGATGAACC‑GAGTATGG‑CGGCCCTGG‑AGCAATGTGC‑TCACAGCCCGCAGGTGATTGCGCTGTTACACGAGTAACACCACCACCCAAGCGAACGCAAACAGCCGTTgg                                                                                                                                >  2:331664/1‑251 (MQ=255)                                                                                                                  cTGATGCGCGGGCTTGTCGAACCGCAGCATGAGCGTGCCTCATGAGCGAAACTTTTTTCCATCTGCTGGGGCCAGGAACGCAACCG‑AACGATGACAGTTTCAGCATGAATCCACTGCCGATCACCTATCAGGTGAATGATGAACC‑GAGTATGG‑CGGCCCTGG‑AGCAATGTGC‑TCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAACATCAACTAAGCGAACGCCAACCGCCGTTGGGCGAGGTg                                                                                                                        >  1:358862/1‑251 (MQ=255)                                                                                                                   tGATGCGCGGGCTTGTCGAACCGCAGCAGGAGCGTGCCTCATGAGAGAAACTTTTTTCCATCTGCTGGGGCCAGGAACGCAACCG‑AACGATGACAGGTTCAGCATGAATCCACTGCCGATCACCTATCAGGTGAATGATGAACC‑GAGTATGG‑CGGCCCTGG‑AGCAATGTGC‑TCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAACATCAACTAAGCGAACGCCAACCGCCGTTGGGCGAggggt                                                                                                                       >  2:427442/1‑248 (MQ=255)                                                                                                                        cgcgGGCTTGTCGAACCGCAGCAGGAGCGTGCCTCATGAGCGAAACTTTTTTCCATCTGCTGGGGCCAGGAACGCAACCG‑AACGATGACAGTTTCAGCATGAATCCACTGCCGATCACCTATCAGGTGAATGATGAACC‑GAGAATGG‑CGGCCCTGG‑AGCAATGTGC‑TCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAACATCAACTAAGCGAACGCCAACCGCCGTTGGGCGAGGTGCTGGc                                                                                                                   >  1:421263/1‑250 (MQ=255)                                                                                                                                                 aGCGTGCCTCATGAGCGAAACTTTTTTCCATCTGCTGGGGCCAGGAACGCAACCG‑AACGATGACAGTTTCAGCATGAATCCACTGCCGATCACCTATCAGGTGAATGATGAACC‑GAGTATGG‑CGGCCCTGG‑AGCAATGTGC‑TCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAACATCAACTAAGCGAACGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCcg                                                                                         >  1:436266/1‑251 (MQ=255)                                                                                                                                                  gcgtgccacgtgaggcaattttttttCCATCGGCGGGGGCCAGGACCCCACCAG‑AAGGATGACAGTTTCAGCATAAATCCCCGGCCAATCCCCTATCGGGTGAATGATGAACC‑GAGTGTGG‑CG‑CCCTGGGAGCAGTGTCC‑CCACAGCCCGCAGTTGATGGCGTTTTAAAACGAGTTACAACATCAACAAAGCGAACGCCACCCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCcga                                                                                        <  2:101592/231‑1 (MQ=255)                                                                                                                                                          aaTGAGCAAAATTTTTTCCCATGTGGGGGGGCCAGGAACG‑AACCGGAACGAGGACATTTTAAGAATGAACCCCCTGCCGATCACCTAACGGGTGAATGA‑GAACCCGAG‑ATGGGCGGCC‑TGG‑GGGAAGTTGCCTCCCACCCCGGGGGTGATGGCGTTGTAAAACGATTAAAAACATCAACAAGGCAAACGCAAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTc                                                                                <  1:336543/250‑1 (MQ=255)                                                                                                                                                                 gAAACTTTTTTCCATCTGCTGGGGCCAGGAACGCAACCG‑AACGATGACAGTTTCAGCATGAATCCACTGCCGATCACCTATCAGGTGAATGATGAACC‑GAGTATGG‑CGGCCCTGG‑AGCAATGTGC‑TCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAACATCAACTAAGCGAACGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAACTCTCAACGCCGACGATCGTCACTTTAt                                                                         >  2:40353/1‑251 (MQ=255)                                                                                                                                                                            ccATCTGCTGAGGCCAGGAACGCAACCG‑AACGATGACAGTTTCAGCATGAATCCACTGCCGATCACCTATCAGGTGAATGATGAACC‑GAGTATGG‑CGGCCCTGG‑AGCAATGTGC‑TCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAACATCAACTAAGCGAACGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCTtc                                                              >  2:425939/1‑251 (MQ=255)                                                                                                                                                                                    tgGGGCCAGGAACGCAACCG‑AACGATGACAGTTTCAGCATGTATCCACTGCCGATTACCTATCAGGTGAATGATGAACA‑GAGTATGG‑CGGCCCTGG‑AGCAATGTGC‑TCACAGCCCGCAGGGGATTGCGCTGTTAAACGAGTTACGACATCAACTAAGCGAACTCCAACCGCCGTTCGGTGCGTGGCCTTCATGCGTTCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATTCGCGTCTCggggg                                                       >  1:376248/1‑249 (MQ=255)                                                                                                                                                                                       ggCCAGGAACGCAACCG‑AACGATGACAGTTTCAGCATGAATCCACTGCCGATCACCTATCAGGTGAATGATGAACC‑GAGTATGG‑CGGCCCTGG‑AGCAATGTGC‑TCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAACATCAACTAAGCGAACGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCTTCTCggggaaggg                                                   >  2:224396/1‑251 (MQ=255)                                                                                                                                                                                        gcaagaaccGCAACCG‑AACGATGACAGTTTAAGCATGAACCCACTGCCGATCACCTATCAGGTGAATGATAAGCG‑GATTATGG‑CGGCCCGGG‑AGAATTGTGC‑TCACCGCCGGCAGGTGATGGCGTTGTTAAACGAGTAACAACATCAAAAAAGCGAACCCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCTTCTCggggaagggg                                                  <  1:170278/243‑1 (MQ=255)                                                                                                                                                                                               aCGCAACCG‑AACGATGACAGTTTCAGCATGAATCCACTGCCGATCACCTATCAGGTGAATGATGAACC‑GAGTATGG‑CGGCCCTGG‑AGCAATGTGC‑TCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAACATCAACTAAGCGAACGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCTTCTCGGGGAAGGGGAAGTGTc                                           >  2:349348/1‑251 (MQ=255)                                                                                                                                                                                                         aaCGATGACAGTTTCAGCATGAATCCACTGCCGATCACCTATCAGGTGAATGATGAACC‑GAGTATGG‑CGGCCCTGG‑AGCAATGTGC‑TCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAACATCAACTAAGCGAACGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCTTCTCGGGGAAGGGGAAGTGTCAGTGCGCAt                                  >  2:334934/1‑251 (MQ=255)                                                                                                                                                                                                              ggaatttttCAGCATGACCCCCCTGCCGTTCACCTACGAGGTGAAGGAGGACCC‑GAGTATGG‑CGGCCCGGG‑AGAAATGTCC‑CCACAGCCCGAGGGTGTTGGCGCTGTTAACCGAGTTACAACACCAACAAAGCGACCGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCTTCTCGGGGAAGGGGAAGTGTCAGTGCGCATTcagc                             <  1:158016/245‑1 (MQ=255)                                                                                                                                                                                                                          cATGAATCCACTGCCGATCACCTATCAGGTGAATGATGAACC‑GAGTATGG‑CGGCCCTGG‑CGCAATGTGC‑TCACAGCCCGCAGGTGCTTGCGCTGTTAAACGAGTTACAACATCAACTAAGCGAACGCCAACCGCAGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCCCAACGCCGACGATCGTCACTTTATCAATACGCTTCTCGGGGAAGGGGAAGTGTCAGGGCGCATTCAGCCGGCTgacgac                  >  1:293630/1‑250 (MQ=255)                                                                                                                                                                                                                           aTGAATCCACTGCCGATCACCTATCAGGTGAATGATGAACC‑GAGTATGG‑CGGCCCTGG‑AGCAATGTGC‑TCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAACATCAACTAAGCGAACGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCTTCTCGGGGAAGGGGAAGTGTCAGTGCGCATTCAGCAGGCTgacgac                  >  1:168499/1‑249 (MQ=255)                                                                                                                                                                                                                                    cTGCCGATCACCTATCAGGTGAATGATGAACC‑GAGTATGG‑CGGCCCTGG‑AGCAATGTGC‑TCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAACATCAACTAAGCGAACGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCTTCTCGGGGAAGGGGAAGTGTCAGTGCGCATTCAGCAGGCTGACGACagtgaaaggg        >  2:95782/1‑248 (MQ=255)                                                                                                                                                                                                                                     tGCCGATCACCTATCAGGTGAATGATAAACC‑GAGTATGG‑CGGCCCTGG‑AGCAATGTGC‑TCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAACATCAACTAAGCGAACGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCTTCTCGGGGCAGGGGAAGTGTCAGTGCGCATTCAGCAGGCTGACGCCagtgaaagtga       >  2:118105/1‑250 (MQ=255)                                                                                                                                                                                                                                       ccGATCACCTATCAGGTGAATGATGAACC‑GAGTATGG‑CGGCCCTGG‑AGCAATGGGC‑TCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAACATCAACTAAGCGAACGCCAACCGCCGTTGGGCGAGGGGCTGGCAGTCGATCTGTTAAATCTCCACGCCGACGATCGTCACTTTATCAATACGCTTCTCGGGGAAGGGGAAGTGTCAGTGCGCATTCAGCAGGCTGACGACAGTGAAAGTGAAAt    >  1:458404/1‑251 (MQ=255)                                                                                                                                                                                                                                           tCACCTATCAGGTGAATGATGAACC‑GAGTATGG‑CGGCCCTGG‑AGCAATGTGC‑TCACAGCCCGCAGGTGATTGCGCTGTTAAACGCGTTACAACATCAACTAAGCGAACGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCTTCTCGGGGAAGGGGAAGTGTCAGTGCGCATTCAGCAGGCTGACGACAGTGAAAGTGAAATAc  >  1:203372/1‑249 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                       AAGCCATCGGCGATCGCTTTGGCGTCTTTCGCTTTCCTGCCACTTTAGTGTTTACCGGCGGAAACTATCGCGGCGTGCTGAATGGTATTCACCCGTGGGCGGAACTGATAAACCTGATGCGCGGGCTTGTCGAACCGCAGCAGGAGCGTGCCTCATGAGCGAAACTTTTTTCCATCTGCTGGGGCCAGGAACGCAACCG‑AACGATGACAGTTTCAGCATGAATCCACTGCCGATCACCTGTCAGGTGAATGATGAACC‑GAGTATGG‑CGGCCCTGG‑AGCAATGTGC‑TCACAGCCCGCAGGTGATTGCGCTGTTAAACGAGTTACAACATCAACTAAGCGAACGCCAACCGCCGTTGGGCGAGGTGCTGGCAGTCGATCTGTTAAATCTCAACGCCGACGATCGTCACTTTATCAATACGCTTCTCGGGGAAGGGGAAGTGTCAGTGCGCATTCAGCAGGCTGACGACAGTGAAAGTGAAATAC  >  NC_000913/1036595‑1037076

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