Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,935,489	G→A	A419T (GCA→ACA)	fucP →	L‑fucose:H(+) symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,935,489	0	G	A	92.7%	111.1 / NA	41	A419T (GCA→ACA)	fucP	L‑fucose:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (12/26);  minor base C (3/0);  total (15/26)
Fisher's exact test for biased strand distribution p-value = 4.27e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTTCGCTGGCGGTCATGTGGGCTTAATAGCCCTGACTTTATGCAGCGCCTTTATGTCGATTCAGTACCCAACAATCTTCTCGCTGGGCATTAAGAATCTCGGCCAGGACACCAAATATGGTTCGTCCTTCATCGTTATGACCATTATTGGCGGCGGTATTGTCACTCCGGTCATGGGTTTTGTCAGTGACGCGGCGGGCAACATCCCCACTGCTGAACTGATCCCCGCACTCTGCTTCGCGGTCATCTTTATCTTTGCCCGTTTCCGTTC‑TCAAACGGCAACTAACTGAACATATTTTCCGAATAAAGTGAGGAATCTGTAATGAAAAAAATCAGCTTACCGAAAATTGGTATCCGCCCGGTTATTGACGGTCGTCGCATGGGTGTTCGTGAGTCGCTTGAAGAACAAACAATGAATATGGCGAAAGCTACGGCCGCACTGCTGACCGAGAA  >  NC_000913/2935262‑2935714                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                   ccTTCGCTGGCGGTCATGTGGGCTTAATAGCCCTGACTGTAGGCAGCGCCTTTATGTCGATTCAGTACCCAACAATCTTCTCGCTGGGCATTAAGAATCTCGGCCAGGACACCAAATATGGTTCGTCCTTCAGCGTTAGGACCAGTATTGGCGGCGGTATTGTCACTCCGGTCATGGGTTTTTTCAGTGACGCGGCTGGCAACAGCCCCACTGCTGAACTGATCCCCACACTCTGCTTGGCGCTCatcttt                                                                                                                                                                                                             <  2:186459/251‑1 (MQ=255)           cGGTCATGTGGGCTTAATAGCCCTGACTTTATGCAGCGCCTTTATGTCGATTCAGTACCCAACAATCTTCTCGCTGGGCATTAAGAATCTCGGCCAGGACACCAAATATGGTTCGTCCTTCATCGTTATGACCATTATTGGCGGCGGTATTGTCACTCCGGTCATGGGTTTTGTCAGTGACGCGGCGGGCAACATCCCCACTGCTGAACTGATCCCCACACTCTGCTTCGCGGTCATCTTTATCTTTGccc                                                                                                                                                                                                   <  1:417723/251‑1 (MQ=255)                aTGTGGGCTTAATAGCCCTGACTTTATGCAGCGCCTTTATGTCGATTCAGTACCCAACAATCTTCTCGCTGGGCATTAAGAATCTCGGCCAGGACACCAAATATGGTTCGTCCTTCATCGTTATGACCATTATTGGCGGCGGTATTGTCACTCCGGTCATGGGTTTTGTCAGTGACGCGGCGGGCAACATCCCCACTGCTGAACTGCTCCCCACAATCTGCTTCGCGGTCATCTTTATCTTTCCCCGTTTc                                                                                                                                                                                              >  1:50252/1‑251 (MQ=255)                  gtgGGCTTAATAGCCCTGACTTTATGCAGCGCCTTTATGTCGATTCAGTACCCAACAATCTTCTCGCTGGGCATTAAGAATCTCGGCCAGGACACCAAATATGGTTCGTCCTTCATCGTTATGACCATTATTGGCGGCGGTATTGTCACTCCGGTCATGGGTTTTGTCAGTGACGCGGCGGGCAACATCCCCACTGCTGAACTGATCCCCACACTCTGCTTCGCGGTCATCTTTATCTTTGCCCGTTTcc                                                                                                                                                                                             <  1:4333/250‑1 (MQ=255)                      gCTTAATAGGCCTGACTTTATGCAGCGCCTTTATGTCGATTCAGTACCCAACAATCTTCTCGCTGGGCATTAAGAATCTCGGCCAGGATACCAAATATGGTGCGTCCTTCATCGTTATGACCATTATTGGCGGCGGTATTGTCACTCCGGTCATGGGTTTTGTCAGTGACGCGGCGGGCAACATCCCCTCTGCTGAGCTGATCCCCACACTCTGCTTCGCGGTCATCGTTATCTTTGCCCGTTTCCGTtc‑t                                                                                                                                                                                       <  2:187047/251‑1 (MQ=255)                         tAATAGCCCTGACTTTATGCAGCGCCTTTATGTCGATTCAGTACCCAACAATCTTCTCGCTGGGCATTAAGAATCTCGGCCAGGACACCAAATATGGTTCGTCCTTCATCGTTATGACCATTATTGGCGGCGGTATTGTCACTCCGGTCATGGGTTTTGTCAGTGACGCGGCGGGCAACATCCCCACTGCTGAACTGATCCCCACACTCTGCTTCGCGGTCATCTTTATCTTTGCCCGTTTCCGTTC‑TCa                                                                                                                                                                                     <  1:230714/250‑1 (MQ=255)                          aaTAGCCCTGACTTTATGCAGCGCCTTTATGTCGATTCAGTACCCAACAATCTTCTCGCTGGGCATTAAGAATCTCGGCCAGGACACCAAATATGGTTCGTCCTTCATCGTTATGACCATTATTGGCGGCGGTATTGTCACTCCGGTCATGGGTTTTGTCAGTGACGCGGCGGGCAACATCCCCACTGCTGAACTGATCCCCACACTCTGCTTCGCGGTCATCTTTATCTTTGCCCGTTTCCGTTC‑TCaaa                                                                                                                                                                                   <  1:312438/251‑1 (MQ=255)                                 cTGACTTTATGCAGCGCCTTTATGTCGATTCAGTACCCAACAATCTTCTCGCTGGGCATTAAGAATCTCGGCCAGGACACCAAATATGGTTCGTCCTTCATCGTTATGACCATTATTGGCGGCGGTATTGTCACTCCGGTCTTGGGTTTTGTAAGTGACGCGGCGGGCAACATCCCCACTGCTGAACTGGTCCCCACAGTCTGATTCGCGGACCTCTTTATCTTTCCCCGTTTCCGGTC‑TCAAACGGCaac                                                                                                                                                                            >  2:317272/1‑251 (MQ=255)                                        tATGCAGCGCCTTTATGTCGATTCAGTACCCAACAATCTTCTCGCTGGGCATTAAGAATCTCGGCCAGGACACCAAATATGGTTCGTCCTTCATCGTTATGACCATTATTGGCGGCGGTATTGTCACTCCGGTCATGGGTTTTGTCAGTGACGCGGCGGGCAACATCCCCACTGCTGAACTGATCCCCACACTCTGCTTCGCGGTCATCTTTATCTTTGCCCGTTTCCGTTC‑TCAAACGGCAACTAACTg                                                                                                                                                                      <  1:302743/250‑1 (MQ=255)                                         atgtcgcgcCTGTATGTCGATGCATTACCCAACAATCGTCTTGCTGGGTATTAAGAAACTCGGCCAGGACACCAAATATGGTTAGTCCTTCATCGTTATGGTCATTATTGGCGGGCGTATTGTCACTCCGGTCATGGGTTTTGTCAGTGACGCGGCGTGCAACATCCCCACTGCTGAACTGATCCCCACTCGCTGCTTCGCGGTAATTTTTATCTTTGCCCGTTTCCGTTC‑TCAAACGGCAACTAACTGa                                                                                                                                                                     <  2:109223/245‑1 (MQ=255)                                                  cTTTATGTCGATTCAGTACCCAACAATCTTCTCGCTGGGCATTAAGAATCTCGGCCAGGACACCAAATATGGTTCGTCCTTCATCGTTATGACCATTATTGGCGGCGGTATTGTCACTCCGGTCATGGGTTTTGTCAGTGACGCGGCGGGCAACATCCCCACTGCTGAACTGATCCCCACACTCTGCTTCGCGGTCATCTTTATCTTTGCCCGTTTCCGTTC‑TCAAACGGCCACTAAATGAACATAttttt                                                                                                                                                           >  2:203326/1‑250 (MQ=255)                                                    ttATGTCGATTCAGTACCCAACAATCTTCTCGCTGGGCATTAAGAATCTCGGCCAGGACACCAAATATGGTTCGTCCTTCATCGTTATGACCATTATTGGCGGCGGTATTGTCACTCCGGTCATGGGTTTTGTCAGTGACGCGGCGGGCAACATCCCCACTGCTGAACTGATCCCCACACTCTGCTTCGCGGTCATCTTTATCTTTGCCCGTTTCCGTTC‑TCAAACGGCAACTAACTGAACATATTTTc                                                                                                                                                           <  1:198183/249‑1 (MQ=255)                                                           gattcaatacacacaaaTCTTCTCGCTGGGCATTCAGCATCTCGGCCAGGACACCAAATACGGTCAGTCCTTAATCGATATGACCCTTATTGGCGGCGGTATTGTCACTCCGGTCATGGGTTTTGTTAGTGACGCGGCGGGCAACATCCCCACTGCTGAACTGTTCCCCCCACTCTGCTTCGCGGTCAACTTTATCTTTGCCCGTTTCCGTTC‑TAAAACGCCAAATAACTGAACATATTTTCCGAATaa                                                                                                                                                    >  2:399521/16‑249 (MQ=255)                                                              tCAGTACCCAACAATCTTCTCGCTGGGCATTAAGAATCTCGGCCAGGACACCAAATATGGTTCGTCCTTCATCGTTATGACCATTATTGGCGGCGGTAGTGTCACTCCGGTCATGGGTTTTGTCAGTGACGCGGCGGGCAGCATCCCCACTGCTTAACTGATCCCCACGCTCTGCTGCGCTGTCATCTTTATCTTTGCCCGTTTCCGTTC‑TCAAACGGCAACTAACTGAACATATTTTCCGAATAAAGTg                                                                                                                                                <  1:209890/250‑1 (MQ=255)                                                               cAGTACCCAACAATCTTCTCGCTGGGCATTAAGAATCTCGGCCAGGACACCAAATATGGTTCGTCCTTCATCGTTATGACCATTATTGGCGGCGGTATTGTCACTCCGGTCATGGGTTTTGTCAGTGACGCGGCGGGCAACATCCCCACTGCTGAACTGATCCCCCCACTCTGCTTCGCGGGCCTCTTTATCTTTGCCCGTTTCCGTTC‑TCAAACGGCCACTAACTGAACCTATTTTCCGAATAAtgtg                                                                                                                                                >  2:8595/1‑249 (MQ=255)                                                                aGTACCCAACAATCTTCTCGCTGGGCATTAAGAATCTCAGCCAGGACACCAAATATGGTTCGTCCTTCCTCGTTATGACCATTATTGGCGGCGGTAATGTCACTCCGGTCATGGGTTTTGTCAGTGAAGCAGCGGGCAAACTCCCCCCTGCTGGACTGATCCCCAAACTCTGCTTCGCGGTCCTCTTTTTCTTTGCCCGTTTCCCTTCCTCAA‑CGGCAACTTACCGACCATTTTTTCCCAATAAAGTGAgg                                                                                                                                             >  2:358498/1‑251 (MQ=255)                                                                       aacaaTCTTCTCGCTGGGCATTAAGAATCTCGGCCAGGACACCAAATATGGTTCGTCCTTCATCGTTATGACCATTATTGGCGGCGGTATTGTCACTCCGGTCATGGGTTTTGTCAGTGACGCGGCGGGCAACATCCCCACTGCTGAACTGATCCCCACACTCTGCTTCGCGGTCATCTTTATCTTTGCCCGTTTCCGTTC‑TCAAACGGCAACTAACTGAACATATTTTCCGAATAAAGTGAGGAATCTGt                                                                                                                                      >  2:78675/1‑251 (MQ=255)                                                                        acaaTCTTCTCGCTGGGCATTAAGAATCTCGGCCAGGACACCAAATATGGTTCGTCCTTCATCGTTATGACCATTATTGGCGGCGGTATTGTCACTCCGGTCATGGGTTTTGTCAGTGACGCGGCGGGCAACATCCCCACTGCTGAACTGATCCCCACACTCTGCTTCGCGGTCATCTTTATCTTTGCCCGTTTCCGTTC‑TCAAACGGCAACTAACTGAACATATTTTCCGAATAAAGTGAGGAATCTGTa                                                                                                                                     <  1:16420/251‑1 (MQ=255)                                                                               tctcGCTGGGCATTAAGAATCTCGGCCAGGACACCAAATATGGTTCGTCCTTCATCGTTATGACCATTATTGGCGGCGGTATTGTCACTCCGGTCATGGGTTTTGTCAGTGACGCGGCGGGCAACATCCCCACTGCTTAACTGATCCCCACACTCTGCGTCGCGGTCATCTTTATCTTTGCCCGTTTCCGTTC‑TCAAACGTCAACTAACTGAACATATTTTCCGAATAAAGTGAGGAATCTGTAATGaaa                                                                                                                               <  2:201999/250‑1 (MQ=255)                                                                                          aTTAATAATCTCTGCCTGGACACCAAATATGGTTCGTCCTTCATCGTTATGACCATTATTGGCGGCGTTATTGGCACTCCGGTCAGGGGTTGTGTCAGGGACGCGGCGGGCAACAGCCCCACTTCTGAACTGATCCCCACACTCGGCTACGCGATCATCTTGATCTTTGCCCGGTTCCGTGC‑TCAAACGGCAACTAACTGAACATAGTTGCCGGATAGAGTGAGGAAGCGGGAATTAAAAAAAttggctta                                                                                                                   <  2:368654/251‑8 (MQ=255)                                                                                               gAATCTCGGCCAGGACACCAAATATGGTTCGTCCTTCATCGTTATGACCATTATTGGCGGCGGTATTGTCACTCCGGTCATGGGTTTTGTCAGTGACGCGGCGGGCAACATCCCCACTGCTGAACTGATCCCCACACTCTGCTTCGCGGTCATCTTTATCTTTGCCCGTTTCCGTTC‑TCAAACGGCAACTAACTGAACATATTTTCCGAATAAAGGGAGGAATCTGTAATGAAAAAAAGCAGCTTACCGaa                                                                                                              <  1:105595/251‑1 (MQ=255)                                                                                                     cGGCCATGACGCCAAATACATTTCGTCCTTCATCGTTATGACCATTATTGGCGGCGGTATTGGCCCTCCGGTCATGGGTTTTGTCAGTGACGCGGCGGGCAACACCCCCACTGCTGAACTGATCCCCACACTCTGCTTCGCGGTCATCTTTATCTTTGCCCGTTTCCGTTC‑TCAAACGGCAACTAACTGAACAGATTTTCCTAATAAAGTGAGGAAACTGTAATGAAAAAAAACACCTTACCGAAAAtt                                                                                                          <  2:360153/249‑1 (MQ=255)                                                                                                         cAGGACACCAGATATGGTTCGTCCTTCATCGTTAGGATCATTATTGGCGGCGGTATTGTCACTCCGGTCATGGTGTTTGTCAGTGACGCGTCGGGCAACATCCCCACGGCTTAACTGATCCGTACTCTCTTCTTCGCGGTCATCTTTATCTTTGCCCTTTTCCGTTC‑TCAAACGGCaac                                                                                                                                                                            <  2:412860/179‑1 (MQ=255)                                                                                                         cAGGACACCAAATATGGTTCGTCCTTCATCGTTATGACCATTATTGGCGGCGGTATTGTCACTCCGGTCATGGGTTTTGTCAGTGGCGCGGCGGGGAACATCCCCACTGTTGCACTAATTTCCACACTATGCTTCGCGGTCATCTTTATCTTTGCCCGTTTCCGTTC‑TCAAACGGCaac                                                                                                                                                                            >  1:412860/1‑179 (MQ=255)                                                                                                               acCAAATATGGGTCGTCCTTCATCGTTATGACCATTATTGGCGGCGGTATTGTCACTCCGGTCCTGGGTTTTGTCAGTGACGCGGCGGGCAACATCCCCACTGCTGAACTGATCCCCACACTCTGCTTCGCGGTCATCTTTATCTTTGCCCGTTTCCGTTC‑TCAAACGGCAACTAACTGAACATATTTTCCGAATAAAGTGAGGAATCTGTAATGAAAAAAATCAGCTTACCGAAAATTGGTATCCGcccc                                                                                              >  2:158012/1‑250 (MQ=255)                                                                                                                      atGGTTCGTCCTTCATCGTTATGACCATTATTGGCGGCGGTATTGTCACTCCGGTCATGGGTTTTGTCAGTGACGCGGCGGGCAACATCCCCACTGCTGAACTGATCCCCACACTCTGCTTCGCGGTCATCTTTATCTTTGCCCGTTTCCGTTC‑TCAAACGGCAACTAACTGAACATATTTTCCGAATAAAGTGAGGAATCTGTAATGAAAAAAATCAGCTTACCGAAAATTGGTATCCGCCCGGTTATTg                                                                                       >  2:16683/1‑251 (MQ=255)                                                                                                                            cGTCCTTCATCGTTATGACCATTATTGGCGGCGGTATTGTCACTCCGGTCATGGGTTTTGTCAGTGACGCGGCGGGCAACATCCCCACTGCTGAACTGATCCCCACACTCTGCTTCGCGGTCATCTTTATCTTTGCCCGTTTCCGTTC‑TCAAACGGCAACTAACTGAACATATTTTCCGAATAAAGTGAGGAATCTGTAATGAAAAAAATCAGCTTACCGAAAATTGGTATCCGCCCGGTTATTGACGgtc                                                                                 <  1:396631/251‑1 (MQ=255)                                                                                                                                           ggACCATTATTGGCGGCGGTATTGTCACTCCGGTCATGGGTTTTGTCAGTGACGCGGCGGGCAACATCCCCACTGCTGAACTGATCCCCACACTCTGCTTCGCGGTCATCTTTATCTTTGCCCGTTTCCGTTC‑TCAAACGGCAACTAACTGAACATATTTTCCGAATAAAGTGAGGAATCTGTAATGAAAAAAATCAGCTTACCGAAAATTGGTATCCGCCCGGTTATTGACGGTCGTCGCATGGGTGTTc                                                                  <  1:203326/250‑1 (MQ=255)                                                                                                                                             aCCATTATGGGCGGCGGTATTGTCACTCCGGTCATGGGTTTTATCAGTGACGCGGCGGGCAACATCCCCACTGCTGAACTGATCCCCACACTCTGCTTCGCGGTCATCTTTATCTTTGCCCGTTTCCGTTC‑TCAAACGGCAACTAACTGAACATATTTTCCGAATAAAGTGAGGAATCTGTAATGAAAAAAATCAGCTTACCGAAAATTGGTATCCGCCCGGTTATTGACGGTCGTCGCATGGGTGTTCGt                                                                <  1:78675/251‑1 (MQ=255)                                                                                                                                                  tattGGCGGCGGTATTGTCACTCCGGTCATGGGGTTTGTCAGTGACGCGGCGGGCAAAATCCCCACTGCTGAACTGATCCCCACACTCTGCTTCGCGGTCATCTTTATCTTTGCCCGTTTCCGTTC‑TCAAACGGCAAATAACTGAACATATTTTCCGAATAAAGTGAGGAATCTGTAATGAAAAAAATCAGCTTACCGAAAATTGGTATCCGCCCGGATATTGACGGGCGTAGCATGGGTATTCGTGAggc                                                           >  1:252358/1‑249 (MQ=255)                                                                                                                                                   attGTCGGCGGTATTGTCACGCCGGTCATGGGTTTTGTCAGTGACGCGGCGGGCAACATCCCCACTGCTGAACTGATCCCCACACTCTGCTTCGCGGTCATCTTTATCTTTGCCCGTTTCCGTTC‑TCAAACGGCAACTAACTGAACATATTTTCCGAATAAAGTGAGGAATCTGTAATGAAAAAAATCAGCTTACCGAAAATTGGTATCCGCCCGGTTATTGACGGTCGTCGCATGGGTGTTCGTGAGTc                                                           <  2:216332/250‑1 (MQ=255)                                                                                                                                                   attGGCGGCGGTATTGTCCCTCCGGTCATGGGTTTTGTCAGTGACGCGGCGGGCAACATCCCCACTGCTGAACTGATCCCCACACTCTGCTTCGCGGTCATCTTTATCTTTGCCCGTTTCCGTTC‑TCAAACGGCAACTAACTGAACATATTTTCCGAATAAAGTGAGGAATCTGTAATGAAAAAAATCAGCTTACCGAAAATTGGTATCCGCCCGGTTATTGACGGTCGTCGCATGGGTGTTCGTGAGTc                                                           <  2:113148/250‑1 (MQ=255)                                                                                                                                                        cggcggTATTGTTACTACGGTCATGGGGTTTGTCAGTGCCGCGGCGGGCCACATCCCCCCTGCTGAACAGATCCCCCCACTCTGCTTCGCCGTCATCTTTATCTTTGCCCGTTTCCGTTC‑TCAAACGGCAACTAACTGAACATATTTTCCGAATAAAGTGAGGAATCTGTAATGAAAAAAATCAGCTTTCCGAAAATTGGTATCCGCCCGGTAATTGACGGGCGTCGCATGGGGGTTCGTGAGTCGCTTg                                                      >  2:80179/1‑250 (MQ=255)                                                                                                                                                          gcggTATTGTCACTCCGGTCATGGGTTTTGTCAGTGACGCGGCGGGCAACATCCCCACTGCTGAACTGATCCCCACACTCTGCTTCGCGGTCATCTTTATCTTTGCCCGTTTCCGTTC‑TCAAACGGCAACTAACTGAACATATTTTCCGAATAAAGTGAGGAATCTGTAATGAAAAAAATCAGCTTACCGAAAATTGGTATCCGCCCGGTTATTGACGGTCGTCGCATGGGTGTTCGTGAGTCGCTTgaag                                                   <  2:159040/251‑1 (MQ=255)                                                                                                                                                            ggTATTGGCACTCCGGTCATGGGTTTTGTCAGTGACGCGGCGGGCAACATCCCCACTGCTGAACTGATCCCCACACTCTGCTTCGCGGTCATCTTTATCTTTGCCCGTTTCCGGTC‑TCAAACGGCAACTAACTGAACATATTTTCCGAATAAAGTGAGGAATCTGTAATGAAAAAAATCAGCTTACCGAAAATTGGTATCCGCCCGGTTATTGACGGTCGTCGCATGGGTGTTCGTGAGTCGCTTgaagaa                                                 <  2:104407/251‑1 (MQ=255)                                                                                                                                                                     aCTCCGGTCATGGGTTTTGTCAGTGACGCGGCGGGCAACATCCCCACTGCTGAACTGATCCCCACGCTCTGCTTCGCGGTCATCGTTATCTTTGCCCGTTTCCGTTC‑TCAAACGGCAACTAACTGAACATATTTTCCGAATAAAGTGAGGAATCTGTAATGAAAAAAATCAGCTTACCGAAAAGTGGTATCCGCCCGGTGATTGACGGTCGTCGCATGGGTGTTCGTGAGTCGCTTGAAGAACAAACAATg                                        <  2:46287/251‑1 (MQ=255)                                                                                                                                                                        ccGGTCATGGGTTTTGTCAGTGACGCGGCGGGCAACATCCCCACTGCTGAACTGATCCCCACACTCTGCTTCGCGGTCATCTTTATCTTTGCCCGTTTCCGTTC‑TCAAACGGCAACTAACTGAACATATTTTCCGAATAAAGTGAGGAATCTGTAATGAAAAAAATCAGCTTACCGAAAATTGGTATCCGCCCGGTGATTGACGGTCGTCGCATGGGTGTTCGTGAGGCGCTTGAAGAGCAAACAATGa                                       <  2:397726/249‑1 (MQ=255)                                                                                                                                                                                        tCAGTGACGCGGCGGGCAACATCCCCACTGCTGAACTGATCCCCACACTCTCCTTCGCGGTCATCTTTCTCTTTGCCCGTTTCCGTTC‑TCAAACGGCAACTAACTGAACATATTTTCCGAATAAAGTGAGGAATCTGTAATGAACAAAATCAGCTTACCGAAAATTGGTATCCGCCCGGTTATTGACGGTCGTCGCATGGGGGTTCGTGAGTCGCTTGAAGAACAAACAATGAATATGGCGAAAGCTACgg                     >  1:24199/1‑251 (MQ=255)                                                                                                                                                                                          aGTGACGCGGCGGGCAACATCCCCACTGCTGAACTGATCCCCACACTCTGCTTCGCGGTCATCTTTATCTTTGCCCGTTTCCGTTC‑TCAAACGGCAACTAACTGAACATATTTTCCGAATAAAGTGAGGAATCTGTAATGAAAAAAATCAGCTTACCGAAAATTGGTATCCGCCCGGTTATTGACGGTCGTCGCATGGGTGTTCGTGAGTCGCTTGAAGAACAAACAATGCATATGGCGAAAGCTACGGcc                   >  2:287252/1‑251 (MQ=255)                                                                                                                                                                                                      ggCAACATCCCCACTGCTGTACTGATCCCCACACTCTGCGTCTCGGTCATCTTTATCTGTGCCCGTTTCCGTTC‑TCAAACGGCAACTAACTGAACATATGTTCCGAATAAAGTGAGGAATCTGTAATGAAAAAAATCAGCTTACCGAAAATTGTTATCCGCCCGGTTATTGACTGTCGTCGCATGGGTGTTCGTGAGTCGCTTTAAGAACAAACAATGAATATGGCGAAAGCTACGGCCGCACTGCTGAcc       <  1:416095/251‑1 (MQ=255)                                                                                                                                                                                                            aTCCCCACTGCTGAACTGATCCCCACACTCTGCTTCGCGGTCATCTTTATCTTTGCCCGTTTCCGTTC‑TCAAACGGCAACTAACTGAACATATTTTCCGAATAAAGTGAGGAATCTGTAATGAAAAAAATCAGCTTACCGACAATTGGTATCCGCCCGGTTATTGACGGTCGTCGCATGGGTGTTCGTGAGTCGCTTGAAGAACAAACACTGAATATGGCGAAAGCTACGGCCGAACTGCGGACCGAGaa  >  2:320636/1‑250 (MQ=255)                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                   CCTTCGCTGGCGGTCATGTGGGCTTAATAGCCCTGACTTTATGCAGCGCCTTTATGTCGATTCAGTACCCAACAATCTTCTCGCTGGGCATTAAGAATCTCGGCCAGGACACCAAATATGGTTCGTCCTTCATCGTTATGACCATTATTGGCGGCGGTATTGTCACTCCGGTCATGGGTTTTGTCAGTGACGCGGCGGGCAACATCCCCACTGCTGAACTGATCCCCGCACTCTGCTTCGCGGTCATCTTTATCTTTGCCCGTTTCCGTTC‑TCAAACGGCAACTAACTGAACATATTTTCCGAATAAAGTGAGGAATCTGTAATGAAAAAAATCAGCTTACCGAAAATTGGTATCCGCCCGGTTATTGACGGTCGTCGCATGGGTGTTCGTGAGTCGCTTGAAGAACAAACAATGAATATGGCGAAAGCTACGGCCGCACTGCTGACCGAGAA  >  NC_000913/2935262‑2935714

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