Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,317,129	G→A	Q32* (CAG→TAG)	trpA ←	tryptophan synthase subunit alpha

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,317,129	0	G	A	97.5%	120.4 / NA	40	Q32* (CAG→TAG)	trpA	tryptophan synthase subunit alpha
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (16/23);  minor base T (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

CATCAATGCCTTTGTTAAACACCAGATTGGCATACATCAACAGGCCAATGGGAATGGTCGGGTGTTTCTGGCGAATCAGTGCCAGCATTTCAAAA‑CATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCTTTATCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCC  >  NC_000913/1316886‑1317369                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 cctcaATGCCTTTGTTAAACACCAGATTGGCATACATCAACAGGCCAATGGGACTGGTCGAGCGTTTCTGGCGAATCAGTGCCAGCATTTCAAAA‑CATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACAGGCTTCAATTAGCGTATCGATCATTTTCAATGACTACTCAATg                                                                                                                                                                                                                                           >  1:62742/3‑251 (MQ=255)                    cacCAGATTGGCATACATCAACAGGCCAATGGGAATGGTCGGGGGTTTCTGGCGAATCAGTGCCAGCATTTCAAAA‑CATTGTGCCGGAGTAACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGGTACCTAACTCCAGCGCGACAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTACTCAATGCCCGGCTCACCGAGCGTg                                                                                                                                                                                                                         >  1:275492/1‑250 (MQ=255)                     acCAGATTGGCATACATCAACAGGCCAAGGGGAGTGGTCGGGTGTTTCGGGCGACTCAGTGCCAGCATTTCAAAA‑CATTGGGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACGCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGAc                                                                                                                                                                                                                       <  1:42823/251‑1 (MQ=255)                        aGATTGGCATACATCAACAGGCCAATGGGAATGGTCGGGTGTTTCTGGAGAATCAGTGCCAGCATTTCAAAA‑AATTGTGCCGGTGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGGGGGTCGGAGAAAGTGTTACATAACTCCAGCGCGACAGAAACAGCTTCAATTAGCGTATCGCTAATTTTCAATGACTACTCAAAGCCCGGATCAcggcgctggacga                                                                                                                                                                                                                     >  2:416662/1‑238 (MQ=255)                         gATTGGCATACATCAACAGGCCAATGGGAATGGTCGGGTGTTTCTGGCGAATCAGTGCCAGCATTTCAAAA‑CATGGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGaaa                                                                                                                                                                                                                   <  1:414006/251‑1 (MQ=255)                           ttGGCATACATCAACAGGCCAATGGGAATGGTCGGGTGTTTCTGGCGAATCAGTGCCAGCATTTCAAAA‑CATTGTGCCGGAGTCACCCCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTACTc                                                                                                                                                                                                                                               <  1:332552/221‑1 (MQ=255)                           ttGGCATACATCAACAGGCCAATGGGAATGGTCGGGTGTTTCTGGCGAATAAGTGCCAGCATTTCAAAA‑CATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTACTc                                                                                                                                                                                                                                               >  2:332552/1‑221 (MQ=255)                             ggCATACATCAACAGGCCAATGGGAATGGTCGGGTGTTTCTGGCGAATCAGTGCCAGCATTTCAAAA‑CATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATAGTCGGGCAATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTCCAATTCGCGTATCGATAATTTTCAATGACTACTCAATGCCCGGATAACCGAGCGttgagaaagga                                                                                                                                                                                                                >  2:115276/1‑240 (MQ=255)                                              cAATGGGAATGGTCGGGTGTTTCTGGCGAATCAGTGCCAGCATTTCAAAA‑CATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGTAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTg                                                                                                                                                                                              >  1:396091/1‑251 (MQ=255)                                                 tGGGAATGGTCGGGTGTTTCTGGCGAATCAGTGCCAGCATTTCAAAA‑CATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTgcgc                                                                                                                                                                                           >  1:184116/1‑251 (MQ=255)                                                 gggaattgggcGGGTGTTTCTGGCGAATCAGTGCCAGCATTCAAAAA‑AATTGTGCCGGAGTCCCCCCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCGTCAATTAGCGTAGCGATAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTgcgc                                                                                                                                                                                           <  2:299322/241‑1 (MQ=255)                                                   ggattggcggggTGTTTCTGGCGAATCAGTGCAAGCATTCCAAAC‑CATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGAGCCTTCTTTGCGCTc                                                                                                                                                                                         <  1:199020/242‑1 (MQ=255)                                                                             cAGTGCCAGC‑TTTCAAAAACATTGGCCCGGAGTCACCCCGGCCGCAAAGGCGCGCAGAGGGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCg                                                                                                                                                                <  2:164160/250‑1 (MQ=255)                                                                                          tCAAAA‑CATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCa                                                                                                                                                   <  1:121206/250‑1 (MQ=255)                                                                                          acAAAA‑CATTGTCCCGGAGTCCCCCCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCa                                                                                                                                                   <  2:184116/249‑1 (MQ=255)                                                                                                       gccgggcccaccccTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGTAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGt                                                                                                                                      <  2:396091/239‑1 (MQ=255)                                                                                                          ggAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGCTACCTACCTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCt                                                                                                                                   >  2:379762/1‑251 (MQ=255)                                                                                                               cacccccgCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCaa                                                                                                                              <  2:217679/244‑1 (MQ=255)                                                                                                                   ccTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAata                                                                                                                          <  1:42098/251‑1 (MQ=255)                                                                                                                    cTGCCGCAAAGGCGCGCAGAGAGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAatat                                                                                                                         >  2:289038/1‑251 (MQ=255)                                                                                                                      gCCGCAAATGCGCGCAGAGGGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTACTCAATGCCCGGATCACCGATCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATc                                                                                                                        <  1:93427/250‑1 (MQ=255)                                                                                                                          cAAAGGCGCGCAGAGTGGCGTTTTGAATCTTCGGGCCAACCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCATCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGaa                                                                                                                   <  1:46483/251‑1 (MQ=255)                                                                                                                            aacgcgcgcAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGCGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGGCTACTCAATGCCCGTATCACCGAGCGGGACGAAGGGAAGGAATGCGCCTTCTTTGCGCTCATTCAACTGGGAAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCTTGAACg                                                                                                                 >  2:296980/4‑251 (MQ=255)                                                                                                                              ggCGCGCAGAGTGGCGTGTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCATATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGt                                                                                                               <  1:287179/251‑1 (MQ=255)                                                                                                                                            ggttttagatcggctggCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCGCGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCTTTAt                                                                                                 <  2:332653/236‑1 (MQ=255)                                                                                                                                                 tGAATCGTCGGGCCATCCGCCAGTGGGTCGTAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCTTTATcgcc                                                                                             <  1:221254/250‑1 (MQ=255)                                                                                                                                                  gAATCGGCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTACTCAATGCCCGGAGCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCTTTATcgccg                                                                                            <  2:35250/250‑1 (MQ=255)                                                                                                                                                               aTCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCGAAATATCGTGAACGGTGAAGATGTCTTTATCGCCGCGACCGGAAAGGtt                                                                              <  1:388694/251‑1 (MQ=255)                                                                                                                                                                  cGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCTTTATCGCCGCGACCGGAAAGGTTAac                                                                           <  1:379762/251‑1 (MQ=255)                                                                                                                                                                       gTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCTTTATCGCCGCGACCGGAAAGGTTAAccacca                                                                      <  1:402777/251‑1 (MQ=255)                                                                                                                                                                                 gaAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCTTTATCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGctc                                                            >  1:143906/1‑251 (MQ=255)                                                                                                                                                                                   aGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATATTTTTCAATGACTTCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCTTTATCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCtt                                                          <  1:364477/251‑1 (MQ=255)                                                                                                                                                                                                   cAGCGCGTCAGCACCGGCTGCAATTAGCGTATCGATAATTTTCAATGACTACGCAATGCCCGGATCACCGAGCGTGACGAAAGGACCGAATGCGCCTTCTTTGCGCTCCTACAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTTAAGATGTCTTTAGCGCCGCGACCGGAAAGGTTAACCCCCAGTAGCTGCTCTTTATCCGGGTTTTcgcg                                          <  1:202776/251‑1 (MQ=255)                                                                                                                                                                                                            aGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCTTTATCGCCGCGACCGGAAAGGTTAACCCCCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCAtttt                                 >  2:277153/1‑251 (MQ=255)                                                                                                                                                                                                            aGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCTTTATCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCAtttt                                 <  1:103651/251‑1 (MQ=255)                                                                                                                                                                                                                      tCAATTAGCTGATCGATAATTTTCAATGACGACTCAATGCCCTGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTGCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCTTTAGCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCTGTAGCCGGGTTTTCGCGCATCATTTTCAACGCATg                        <  2:275492/250‑1 (MQ=255)                                                                                                                                                                                                                             gCGTATCGATAATTTTCAATGACTACTCAATGCCCGGATCACCGCGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCTTTATCGCCGCGACCGGAAAGGTTAACCACCAGTAGATGCTCTTTATCCGGGATTTCGCGCATCATTTTCAACGCATGGGCCCggg                >  1:377308/1‑251 (MQ=255)                                                                                                                                                                                                                                  tCGATAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCTTTATCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCTTGGGCCAGGGCGTgg           >  2:94362/1‑251 (MQ=255)                                                                                                                                                                                                                                     aTAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCTTTATCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTggg          >  2:375535/1‑249 (MQ=255)                                                                                                                                                                                                                                             cAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAAAGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCTTTATCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTcc  >  2:244747/1‑249 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 CATCAATGCCTTTGTTAAACACCAGATTGGCATACATCAACAGGCCAATGGGAATGGTCGGGTGTTTCTGGCGAATCAGTGCCAGCATTTCAAAA‑CATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCTTTATCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCAACGCATGGGCCAGGGCGTGGGAGGATTCC  >  NC_000913/1316886‑1317369

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