Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,591,074	G→A	L418L (CTG→TTG)	ugpB ←	sn‑glycerol 3‑phosphate ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,591,074	0	G	A	97.8%	134.3 / ‑4.4	45	L418L (CTG→TTG)	ugpB	sn‑glycerol 3‑phosphate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (21/23);  total (22/23)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.72e-01

CGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAGTGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATACGCTGCTTTGGTGATTGGCAGATAACCG  >  NC_000913/3590836‑3591294                                                                                                                                                                                                                                               |                                                                                                                                                                                                                             cGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTgtgg                                                                                                                                                                                                                  <  1:17415/251‑1 (MQ=255)  ggTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGTAACGCATCCGACATTAAACTGCTTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCCCAACGGCGGTATCCAATGCCTGCTggggt                                                                                                                                                                                                                 >  2:6655/1‑251 (MQ=255)    ggatgatgaGCTGAGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATGACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGgtgt                                                                                                                                                                                                               <  2:81077/250‑1 (MQ=255)       tgatgaGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGTTGCCAGGTCAGATAAGCGCAACGCATCCGACCTTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCACCGGCGGTATCCAATTCCTTCTgtggggtttt                                                                                                                                                                                                            >  2:42663/1‑245 (MQ=255)             gCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTAcc                                                                                                                                                                                                      <  2:136273/251‑1 (MQ=255)                   gcgcGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCca                                                                                                                                                                                                <  1:277486/251‑1 (MQ=255)                       gACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCacacc                                                                                                                                                                                            >  1:174748/1‑250 (MQ=255)                           agcagATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGctct                                                                                                                                                                                        <  1:278757/251‑1 (MQ=255)                                    ggCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCATCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCtt                                                                                                                                                                               <  1:226394/251‑1 (MQ=255)                                                  gCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCAcg                                                                                                                                                                 <  1:405755/251‑1 (MQ=255)                                                   gggAACACCGGACGGGATGATGACATTACGGTTAATTCCGGAACTGGGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCAGCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCAcgc                                                                                                                                                                <  2:248906/250‑1 (MQ=255)                                                           cGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGc                                                                                                                                                        <  1:480596/251‑1 (MQ=255)                                                                gggATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGCTTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGACCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCAt                                                                                                                                                   >  1:321957/1‑251 (MQ=255)                                                                gggATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCAt                                                                                                                                                   <  2:23600/251‑1 (MQ=255)                                                                  gatgatGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCGTCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATACCTGCTGTGGTTTCTTCTTACCGGTTCACACGCTCTCCAGCTCTTCCTCCACAATCCCGCGGATCTGCGGCATGt                                                                                                                                                 >  2:14189/1‑251 (MQ=255)                                                                        gaCATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGCAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGCTCTGCGGCATGTTGccc                                                                                                                                            >  1:51384/1‑250 (MQ=255)                                                                              aCGGTTAATTCCTGAACTGTGCCGGATGCGCTTAGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGc                                                                                                                                      >  1:377379/1‑250 (MQ=255)                                                                                    aaTTCCTGAACTGTGCAGGATGCGCTTCGCTTATCCGGCATACACGCGAAGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACAGGTCCACACGCTCTCCAACTCTTCATCCACAATCACGAGGATCTGCGGCATGTTGCCCTGACGCAACCCtt                                                                                                                               >  1:168065/1‑251 (MQ=255)                                                                                           gAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGaa                                                                                                                        <  2:218493/251‑1 (MQ=255)                                                                                           gAACTGTGCCGGATGCGCTTCGCTTATACGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTTTGGTGTCTTCTTACCGGTACACACGATCTCCAGCTCTTCAACCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCCGCCATTTggggac                                                                                                                        >  1:345817/1‑247 (MQ=255)                                                                                            accTGTGCCGGATTGGCTTCCCTTAGCTGTCCTACACGCGATGCCAGGGCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTGCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGaac                                                                                                                       <  2:321957/249‑1 (MQ=255)                                                                                                  gCCGGATGCGCTTCGCTTATCCGGCCTACCCGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTCaacggcaac                                                                                                                 <  2:102494/251‑1 (MQ=255)                                                                                                              tcgcttATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTAtt                                                                                                      <  1:273412/250‑1 (MQ=255)                                                                                                                 cttATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCttttttacc                                                                                                  >  2:74722/1‑244 (MQ=255)                                                                                                                             acacGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTGCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTc                                                                                      <  1:449123/251‑1 (MQ=255)                                                                                                                                                  aaGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGttt                                                                 <  1:471402/251‑1 (MQ=255)                                                                                                                                                  aaGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGttt                                                                 <  1:363617/251‑1 (MQ=255)                                                                                                                                                  aaGCGAAACGCACCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAAGTGATGCCCACGCCCACCGGCGGTAACCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCAGCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGATTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGttt                                                                 <  2:28655/251‑1 (MQ=255)                                                                                                                                                     cGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGGCGCCGCCCTTTGGTGAACGGCAACGGCGGCTTATTCCGCATCTGCCGCGCCGCGGGATCCGCCCCTGGGtttttt                                                              >  2:419390/1‑250 (MQ=255)                                                                                                                                                     cGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTtc                                                              <  1:74722/251‑1 (MQ=255)                                                                                                                                                                   cATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGCACGGCAACGGCGGCTTATTCAGCATCTGCCGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAACCCtc                                                >  1:457689/1‑250 (MQ=255)                                                                                                                                                                       acacTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCCGGGGTTTTCTCCTAAAAGCCCTGCTcc                                            >  1:289485/1‑250 (MQ=255)                                                                                                                                                                               tAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGt                                                                                 <  1:451692/206‑1 (MQ=255)                                                                                                                                                                                 aGACTTCGTCGATTTCGCAAAGCGTCGCAGCAACTTATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTTGTGTCTTCTTACCTGTCCACACGCTCTCCAGCTCTTCAGCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGTCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCTCTGGGTTTTTCTCATAAAAGCCCTGCTCACGggtcaggt                                  <  1:293369/251‑1 (MQ=255)                                                                                                                                                                                   aCTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGggtcaggtc                                 >  1:340490/1‑250 (MQ=255)                                                                                                                                                                                    cTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCCGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGCGAACGGCAACGGCGGCTTCTTCAGCATCTGACGCGTACCGGGATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTAAGGTCAt                               >  1:299305/1‑251 (MQ=255)                                                                                                                                                                                       cgtcgATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATATGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTCTTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATACg                            >  1:8321/1‑251 (MQ=255)                                                                                                                                                                                          cgATTTCTCAACGCGGCGCAGCAACTGATTTCCACGCTAAACCGCGGTATCCAATGCCCGCTGGGGTGTCTTCTTACCGGTCCACCCGCTCTCCAGCTCTTCATCCCCAATCTCGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGACCTGCAACGGCGGCTTTTTCAGCATCTGACGCGTCGCGACATCCGCCCCTGGGTATTTCTCCTAAAAACCCTGCGCCCGgggcaggagctaccctg                         >  2:226501/1‑241 (MQ=255)                                                                                                                                                                                              ttCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAGTGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGCCATACGCTGCttt                     >  1:181666/1‑251 (MQ=255)                                                                                                                                                                                                     aGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTCAGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCCGCATCTGACGCGACGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGCTCACGTCATACGCCGATTTGGTGAtt              >  1:63949/1‑251 (MQ=255)                                                                                                                                                                                                          cgcAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCCTCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAACCCTGCTCCCGGGTCAGGTCATACGCTGCTTTGGGGATTGGccg         >  2:374791/1‑249 (MQ=255)                                                                                                                                                                                                                 aCTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCCTAAAAGCCCTGCTCACAGGGCAGGTCATACGCTGCtttagttcggaaaggcggcggg  >  2:69583/1‑232 (MQ=255)                                                                                                                                                                                                                 aCTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATACGCTGCttt                     <  1:69583/232‑1 (MQ=255)                                                                                                                                                                                                                                         tATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTa                                                                                >  2:110635/1‑149 (MQ=255)                                                                                                                                                                                                                                         tATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTa                                                                                <  1:110635/149‑1 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                             CGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAGTGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCCTTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATACGCTGCTTTGGTGATTGGCAGATAACCG  >  NC_000913/3590836‑3591294

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