Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,939,297	G→T	G137V (GGG→GTG)	fucU →	L‑fucose mutarotase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,939,297	0	G	T	97.9%	160.2 / NA	47	G137V (GGG→GTG)	fucU	L‑fucose mutarotase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base T (21/25);  minor base A (1/0);  total (22/25)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

GATTATCCCGTTATTTGAACTGGACAGTTATGCACCGCCGCTGGTGATGATGGCGGCGGTAGAAGGTGACACTCTCGATCCTGAAGTAGAACGACGTTACCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCATTCGTCGCGATCTCAATGAATTACAGACGCAGGGTAAAATTCTGC  >  NC_000913/2939055‑2939515                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                           gATTATCCCGTTATTTGAACTGGACAGTTATGCACCGCCGCTGGTGATGATGGCGGCGGTAGAAGGTGACACTCTCGATCCTGAAGTAGAACGACGTTACCGTAATGCGCTTTCACTACAAGCACCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAccac                                                                                                                                                                                                                    >  1:158247/1‑251 (MQ=255)    tATCCCGTTATTTGAACTGGACAGTTATGCACCGCCGCTGGTGATGATGGCGGCGGTAGAAGGTGACACTCTCGATCCTGAAGTAGAACGACGTTACCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCg                                                                                                                                                                                                                  <  2:114449/250‑1 (MQ=255)      tCCCGTTATTTGAACTGGACAGTTATGCACCGCCGCTGGTGATGATGGCGGCGGTAGAAGGTGACACTCTCGATCCTGAAGTAGAACGACGTTACCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTa                                                                                                                                                                                                                >  2:287998/1‑250 (MQ=255)       cccGTTATTTGAACTGGACAGTTATGCACCGCCGCTGGTGATGATGGCGGCGGTAGAAGGGGACACTCTCGATCCTGAAGTAGAACGACGTTACCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGCGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAAt                                                                                                                                                                                                              >  2:243721/1‑251 (MQ=255)          gTTATTTGAACTGGACAGTTATGCCCCGCCGCTGGTGATGATGGCGGCGGTAGAAGGTGACACTCTCGATCCTGAAGTAGAACGACGTTACCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAAtctc                                                                                                                                                                                                           <  2:276491/251‑1 (MQ=255)               ttGAACTGGACAGTTATGCACCGCCGCTGGTGATGATGGCGGCGGTAGAAGGTGACACTCTCGATCCTGAAGTAGAACGACGTTACCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATAc                                                                                                                                                                                                       <  1:63973/250‑1 (MQ=255)                tGAACTGGACAGTTATGCACCGCCGCTGGTGATGATGGCGGCGGTAGAAGGTGACACTCTCGATCCTGAAGTAGAACGACGTTACCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATAcc                                                                                                                                                                                                      <  1:364171/250‑1 (MQ=255)                    cTGGACAGTTATGCACCGCCGCTGGTGATGATGGCGGCGGTAGAAGGTGACACTCTCGATCCTGAAGTAGAACGACGTTACCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTAc                                                                                                                                                                                                 >  1:59855/1‑251 (MQ=255)                    cTGGACAGTTATGCACCGCCGCTGGTGATGATGGCGGCGGTAGAAGGTGACACTCTCGATCCTGAAGTAGAACGACGTTACCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTAc                                                                                                                                                                                                 >  2:141795/1‑251 (MQ=255)                            ttATGCACCGCCGCTGGTGATGATGGCGGCGGTAGAAGGTGACACTCTCGATCCTGAAGTAGAACGACGTTACCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCAt                                                                                                                                                                                         >  1:208468/1‑251 (MQ=255)                             tATGCACCGCCGCTGGTGATGATGGCGGCGGTAGAAGGTGACACTCTCGATCCTGAAGTAGAACGACGTTACCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATg                                                                                                                                                                                        >  1:13426/1‑251 (MQ=255)                             tATGCACCGCCGCTGGTGATGATGGCGGCGGTAGAAGGTGACACTCTCGATCCTGAAGTAGAACGACGTTACCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATg                                                                                                                                                                                        >  1:291488/1‑251 (MQ=255)                               tGCACCGCCGCTGGTGATGATGGCGGCGGTAGAAGGTGACACTCTCGATCCTGAAGTAGAACGACGTTACCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGAc                                                                                                                                                                                      >  1:171054/1‑251 (MQ=255)                                gaaccactgcaggagatgatGGCGGCGGTAGGAGGTGCCACGCTCGATACCGAAGTAGCACGACGATACCGTAATGCGCTTTCAATACAAGCCACGAGGCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTATGCGATAGTTATCACAGGCAAACGAGCGAAGTACGGGAATATTCTTTTCACAAAAGTGATAACCCAGTTATCTCATAACGGTACGcgcgtataacc                                                                                                                                                                                     >  2:349091/15‑241 (MQ=255)                                gCACCGCCGCTGGTGATGATGGCGGCGGTAGAAGGGGCCACTCTCGATCCTGAAGTAGAACGACGTTACCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACTGGCGCAAAAAGCCTTTGCGCTCGTTATCACAGGCGAACGCGCGAAGTACGGGACTATTCTTTTAAAAAAAGTGGCAACACCGTAATCTCATACCGGTCCGCCCGCATGAcg                                                                                                                                                                                     >  2:349119/1‑251 (MQ=255)                                                   gggcgtcGGTAGAAGGTGACACTCTCGATCCTGAAGTAGAACGACGTTACCGTAATGCGCTTTCACTTCAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAatgat                                                                                                                                                                  <  1:341577/250‑1 (MQ=255)                                                    ggcggcGGTAGAAGGTGACACTCTCGATCCTGAAGTAGAACGACGTTACCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAatgatg                                                                                                                                                                 <  1:137305/251‑1 (MQ=255)                                                                     acacTCTCGATCCTGAAGTAGAACGACGTTACCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGaaa                                                                                                                                                <  2:13426/251‑1 (MQ=255)                                                                      cacTCTCGATCCTGAAGTAGAACGACGTTACCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGaaa                                                                                                                                                >  2:136949/1‑250 (MQ=255)                                                                                ccTGAAGTAGAACGACGTTACCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATAGGAAAGCGGCACGCCa                                                                                                                                     <  1:270798/251‑1 (MQ=255)                                                                                          aaCGACGTTACCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATa                                                                                                                           <  2:291488/251‑1 (MQ=255)                                                                                                   aCCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACctg                                                                                                                  <  1:15025/251‑1 (MQ=255)                                                                                                   aCCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTGTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCGTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACct                                                                                                                   <  1:243721/250‑1 (MQ=255)                                                                                                       tAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACatgctgc                                                                                                              >  2:208802/1‑251 (MQ=255)                                                                                                         aTGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACctgctgctg                                                                                                            >  2:46010/1‑251 (MQ=255)                                                                                                          ggcgcTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGGGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACTGTAATCTCATACCGGTACGCCCGTGTGTCTCGGGCGTTTAGCGAAGGATGGGGTGAAAAATATTAAAGCGGCTCGGTAGCAGGCTATTGTCGACCAGCTGCTGa                                                                                                           <  2:35618/250‑1 (MQ=255)                                                                                                                    cTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATAcc                                                                                                   >  1:90712/1‑249 (MQ=255)                                                                                                                                                    cAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTa                                                                  >  1:13135/1‑250 (MQ=255)                                                                                                                                                         gTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTc                                                            <  1:54293/251‑1 (MQ=255)                                                                                                                                                          ttttGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGc                                                                                                                                              >  2:14426/1‑168 (MQ=255)                                                                                                                                                          ttttGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGc                                                                                                                                              <  1:14426/168‑1 (MQ=255)                                                                                                                                                            ttGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGt                                                         >  2:187434/1‑251 (MQ=255)                                                                                                                                                                gTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTaa                                                       >  2:354912/1‑249 (MQ=255)                                                                                                                                                                     aTGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACCGCTAAAGGTCAGTAAAGAAAcc                                                >  2:362102/1‑251 (MQ=255)                                                                                                                                                                        aaCGGGCGCAAAAAGCCTTTGCGATCATTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATTTCAGACCGGTAAGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAACGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCATGACCACGGAAGCGCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCAtt                                             <  2:289760/251‑1 (MQ=255)                                                                                                                                                                            ggCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCATtcgtc                                         <  1:195979/251‑1 (MQ=255)                                                                                                                                                                            ggCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCATtcgtc                                         <  1:14320/251‑1 (MQ=255)                                                                                                                                                                             gcgcAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCATtcgtcg                                        <  1:205660/251‑1 (MQ=255)                                                                                                                                                                               gcAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCATTCGTcgcg                                      >  1:169389/1‑251 (MQ=255)                                                                                                                                                                                          ttGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCATTCGTCGCGATCTCAATg                             <  2:208468/249‑1 (MQ=255)                                                                                                                                                                                               aTCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCATTCGTCGCGATCTCAATGAATTACa                      <  1:129997/251‑1 (MQ=255)                                                                                                                                                                                                               aaCGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAATGTCAGTAAAGAAACCATTCGTCGCGATCTCAATGAATTACAGACGCAGGGTAAAAtt      <  2:347427/251‑1 (MQ=255)                                                                                                                                                                                                                 cGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCATTCGTCGCGATCTCAATGAATTACAGACGCAGGGTAAAATTCt    <  1:17425/251‑1 (MQ=255)                                                                                                                                                                                                                  gAGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCATTCGTCGCGATCTCAATGAATTACAGACGCAGGGTAAAATTCTg   <  2:299256/251‑1 (MQ=255)                                                                                                                                                                                                                   aGCGAAGTACGGGAATATTCTTTTAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCATTCGTCGCGATCTCAATGAATTACAGACGCAGGGTAAAATTCTgc  <  1:208802/251‑1 (MQ=255)                                                                                                                                                                                                                                          tAAAAAAAGTGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCATTCGTcgcg                                      <  1:147315/192‑1 (MQ=255)                                                                                                                                                                                                                                          tAAAAAAAGAGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCATTCGTcgcg                                      >  2:147315/1‑192 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                           GATTATCCCGTTATTTGAACTGGACAGTTATGCACCGCCGCTGGTGATGATGGCGGCGGTAGAAGGTGACACTCTCGATCCTGAAGTAGAACGACGTTACCGTAATGCGCTTTCACTACAAGCCCCGTGTCCTGACATCATCCGCATCAATCGTTTTGCGTTTTATGAACGGGCGCAAAAAGCCTTTGCGATCGTTATCACAGGCGAACGAGCGAAGTACGGGAATATTCTTTTAAAAAAAGGGGTAACACCGTAATCTCATACCGGTACGCCCGCATGACGCGGGCGGTTATCGAATGATGGGGTGAAAAATATGAAAGCGGCACGCCAGCAAGCGATAGTCGACCTGCTGCTGAACCATACCAGCCTGACCACGGAAGCTCTCTCTGAACAGCTAAAGGTCAGTAAAGAAACCATTCGTCGCGATCTCAATGAATTACAGACGCAGGGTAAAATTCTGC  >  NC_000913/2939055‑2939515

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