Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,788,754	G→A	L218L (TTG→CTA)	yibQ →	divergent polysaccharide deacetylase domain‑containing protein YibQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,788,754	0	G	A	95.6%	135.3 / NA	45	L218L (TTG→CTA)	yibQ	divergent polysaccharide deacetylase domain‑containing protein YibQ
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (23/20);  minor base C (2/0);  total (25/20)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

AGGCGCTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGCACTGGCGTGAAGGTGATCAAACGGAAGGTATTCCTCGACGATTCGCAAAATGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGACCTGGCGCGTCGCAACGGTTCGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCA‑GGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCA  >  NC_000913/3788519‑3788996                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                    aGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTtat                                                                                                                                                                                                                                      >  1:151309/1‑251 (MQ=255)       tGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCCACGGTACGCGTGCTCCAACAGATGGTTTATAACCTg                                                                                                                                                                                                                                >  1:309026/1‑251 (MQ=255)             gTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCCTTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTgccgcc                                                                                                                                                                                                                           >  2:319765/1‑250 (MQ=255)                       ttACTTCCTCGACAGCGTCACCATCGGAAATACCCAGGCGATGCGCGCCGCGCACGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCACCGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGc                                                                                                                                                                                                                >  2:51842/1‑251 (MQ=255)                       ttACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTa                                                                                                                                                                                                                   <  1:82789/248‑1 (MQ=255)                          cTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGCTGCGCGCCGCGCAAGGTACTGGCGTGACGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGAGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTgg                                                                                                                                                                                                             >  2:116614/1‑251 (MQ=255)                                cGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGaa                                                                                                                                                                                                         >  2:57433/1‑249 (MQ=255)                                 gACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCACGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGCCATTACGCTGGTGaacgc                                                                                                                                                                                                      >  1:225001/1‑248 (MQ=255)                                  aCAGCGTAACCATCGGTAATACCCAAGCGCTGCGCGCCGCGCACGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGCCGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGGACGCGTGCTACAACAGATGGTTTCTAACCTGCCGCCAGAAATTACGCTGGTGAAAGcc                                                                                                                                                                                                     >  2:307353/1‑251 (MQ=255)                                             aTCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCt                                                                                                                                                                                          >  1:148962/1‑251 (MQ=255)                                                ggTAATACCCAGGCGATGAGCGCCGCGCAAGGCACTGGCGTGACGGTGATCAAACGGAAGGTGTTCCTCGACGATTTGCCGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCACCCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTCCAACAGATGGTTTATAACCTGCCGCCCGACATTCCGCTGGTGAAAGCCAGCAGCTTGCtga                                                                                                                                                                                        >  2:275752/1‑250 (MQ=255)                                                  tAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCtgaatg                                                                                                                                                                                     >  1:334818/1‑251 (MQ=255)                                                                          gcAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCA‑GGTTGATACTTCTacac                                                                                                                                                            >  1:230097/1‑251 (MQ=255)                                                                                    ggCGTGAAGGTGAACAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCACCCTCCCCCTTCAACGGTACGCGTGCTACAACAGAGGGTTTATAACCGGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCATGCTGAAGGAACCGCAAGGTT‑ATACTTCTACACCGCCGaaaaa                                                                                                                                                  <  2:285170/251‑1 (MQ=255)                                                                                     gCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCA‑GGTTGATACTTCTACACCGCCGAAAAAc                                                                                                                                                 <  1:216387/251‑1 (MQ=255)                                                                                        ggAAGGGGAACAAACGGAAGGTGTTCCCCGGCGAATCGCAGAACCAAAGGGACCTCCGTGTGGAAATTAAACGCGCAAATGAACTGGCGGGGCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCA‑GGTTGATACTTCTACACCGCCGAAAAACGCt                                                                                                                                              <  1:260241/250‑1 (MQ=255)                                                                                                aTCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCA‑GGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTg                                                                                                                                       >  1:65446/1‑250 (MQ=255)                                                                                                      cGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCA‑GGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCAcc                                                                                                                                <  1:57433/251‑1 (MQ=255)                                                                                                       ggAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCCGCAGCTTGCTGAATGAACCGCA‑GGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACcg                                                                                                                               >  1:163254/1‑251 (MQ=255)                                                                                                                 ccTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACACCAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCA‑GGTTGATACTTCTCCACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAAtccgt                                                                                                                     >  1:240472/1‑251 (MQ=255)                                                                                                                     gacgaTTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCG‑GGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAAtccgttccg                                                                                                                 >  1:12818/1‑251 (MQ=255)                                                                                                                            cgcagaacgAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCA‑GGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGt                                                                                                           <  1:198229/242‑1 (MQ=255)                                                                                                                              cagaacgAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCtg                                                                                                                                                                                         <  2:166462/165‑1 (MQ=255)                                                                                                                              cagaacgAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCtg                                                                                                                                                                                         >  1:166462/7‑171 (MQ=255)                                                                                                                                 aacgAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCA‑GGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCt                                                                                                     <  1:51842/248‑1 (MQ=255)                                                                                                                                       gCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCAGTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCA‑GGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTGCCGTGGCGTGAAGCTGTGCaa                                                                                               <  2:230097/251‑1 (MQ=255)                                                                                                                                             aTCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAg                                                                                                                                                                                                                                               >  1:104043/1‑102 (MQ=255)                                                                                                                                             aTCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCA‑GGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCga                                                                                          >  1:30684/1‑250 (MQ=255)                                                                                                                                             aTCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTCCAACGAGTGTCATTGGGCATCCTCACCCTTCAACGGTACCGGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCTCA‑GGTTGATACTGCGACACCGCCGAATAACGCGGTGCCTTATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCgaa                                                                                         <  1:56351/251‑1 (MQ=255)                                                                                                                                                    tgCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCCGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCA‑GGTTGATCCTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAGGAAACCg                                                                                  >  1:430196/1‑251 (MQ=255)                                                                                                                                                            aaTCGCGCAATTGATCTGGCGCGACGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCA‑GGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAACCCGATAGAAc                                                                           <  2:242862/250‑1 (MQ=255)                                                                                                                                                                   cAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCA‑GGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTAt                                                                   >  1:157827/1‑251 (MQ=255)                                                                                                                                                                   aaattgaacgggcgctacGCAACGGTTCAACCATTGCCATTGGGCATCCGCACCCTTCAACGGTACGCGTGCTACAACAGAGGGTTTATAACCTGCCGCCAGACATTACTCTGTTGAAAGCCAGCAGCTTGCTGAATGAACCGCA‑GGTTGATACTTCTACTCCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCTTTCCGTGGCGTGAAGCTGTGCAATCCGAAGAAACCGATAGAACCTGTTTAt                                                                   <  2:397569/234‑1 (MQ=255)                                                                                                                                                                             ggCGCGTCGCAACGGTTCAACCATTGCCATTGGGCCTCCTCACCCTTCGACGGTACGCGTGCTACAAAAGATGGTTTATAACCTGCCGCCAGACATTACTCTGGTGAAAGCCAGCAGCTGGCTGAATGCACCGCA‑GGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTTCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGct                                                         <  2:101628/251‑1 (MQ=255)                                                                                                                                                                             ggCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCA‑GGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAAt                                                             <  2:411961/247‑1 (MQ=255)                                                                                                                                                                                  gTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCA‑GGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCttt                                                    >  1:123377/1‑251 (MQ=255)                                                                                                                                                                                       aaCGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCA‑GGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGt                                               >  1:93561/1‑251 (MQ=255)                                                                                                                                                                                                     tgccattgGGCATCCTCCCCCTTCAACGGTACGCGTGCTAAAACAGAGGGGTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCA‑GGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGGGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATGAAGCGAAAGc                                  <  2:197604/250‑1 (MQ=255)                                                                                                                                                                                                       ccattgGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCA‑GGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATc                               <  1:266222/251‑1 (MQ=255)                                                                                                                                                                                                             ggCATCCTCACCCTTCAACGGTCCGCGAGCTACAACAGCTGGTTTATAACCTGCCGCCCGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCA‑GGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAACCCGAAGAAACCGATAGACCCTGTTTATGCTACTCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGcc                           >  2:178189/1‑249 (MQ=255)                                                                                                                                                                                                                  ccTCACCCTTCAACGGTACGCGTGCGACAACAGATGGGTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCA‑GGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGGGCAAACCTAAGAAACCGATAGAACCAGTGTATGCTGATCGCCTCTTTGAAGTACTAAGCGAAAGGATCAGCCAGAGCa                     <  2:299311/250‑1 (MQ=255)                                                                                                                                                                                                                      aCCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCA‑GGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTg                >  2:134679/1‑251 (MQ=255)                                                                                                                                                                                                                        ccttcAACGGGACGCGTGCTACAACAGATGGTTTATAACCGGCCGCCAGACATGACGCTGGTGAAAGCCAGCAGATGGCGGAATGAACCGCA‑GGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCTTTCCGTGGCGTGAAGCTGTGTAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCGTTGAAGTATTAAGCGAAAGCATCAGCCAGATCACGCTGAt              <  2:364967/247‑1 (MQ=255)                                                                                                                                                                                                                                   aCGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTTCTGAATGAACCGCA‑GGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTc    <  1:347276/250‑1 (MQ=255)                                                                                                                                                                                                                                    cgcgTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCA‑GGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCa  <  1:134679/251‑1 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                    AGGCGCTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGCACTGGCGTGAAGGTGATCAAACGGAAGGTATTCCTCGACGATTCGCAAAATGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGACCTGGCGCGTCGCAACGGTTCGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCA‑GGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCGTGGCGTGAAGCTGTGCAAACCGAAGAAACCGATAGAACCTGTTTATGCTAATCGCTTCTTTGAAGTATTAAGCGAAAGCATCAGCCAGAGCACGCTGATCGTTTACTTCCA  >  NC_000913/3788519‑3788996

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