Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	292,019	C→T	L174L (CTA→TTA)	yagJ →	protein YagJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	292,019	0	C	T	98.0%	156.4 / ‑3.1	51	coding (520/732 nt)	yagJ	protein YagJ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (26/24);  minor base . (0/1);  total (26/25)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01

TTGTTGTGTGACGAAAGCGGGAACTATCCGGCTCTGTTGAACCATTGGGAAACCAAGGTTTTTGAGATCGAAACCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCCTAC‑ATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGCGAATCGAATCTGTCGCCGAAGTAAAGGGTAAATTACGAGTGTTAGATTTGAAACGGCAGGATGTGCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTT  >  NC_000913/291776‑292264                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       ttgttgTGTGACGAAAGCGGGAACTATCCGGCTCTGTTGAACCATTGGGAAACCAAGGTTTTTGAGATCGAAACCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTAC‑ATc                                                                                                                                                                                                                                                 <  2:222181/250‑1 (MQ=255)          gACGAAAGCGGGAACTATCCGGCTCTGTTTAACCATTGGGAAACCAAGGTTTTTGAGATCGAAACCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCTCTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTAC‑ATCTGGCTGATg                                                                                                                                                                                                                                        <  1:612126/250‑1 (MQ=255)                          aTCCGGCTCTGTTGAACCATTGGGAAACCAAGGTTTTTGAGATCGAAACCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTAC‑ATCTGGCTGATGCTTTGCCCAAACTgg                                                                                                                                                                                                                         <  1:388063/249‑1 (MQ=255)                                  ctGTTGAACCATTGGGAAACCAAGGTTTTTGAGATCGAAACCAAACGGGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTAC‑ATCTGGCTGATGCTTTGCCCAAACTGGAAGGACtt                                                                                                                                                                                                                 <  1:486853/249‑1 (MQ=255)                                    gTTGAACCATTGGGAAACCAAGGTTTTTGAGATCGAAACCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTAC‑ATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTgcgc                                                                                                                                                                                                             >  1:63487/1‑251 (MQ=255)                                    gTTGAACCATTGGGAAACCAAGGTTTTTGAGATCGAAACCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTAC‑ATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTgcg                                                                                                                                                                                                              <  1:93196/250‑1 (MQ=255)                                    gTTGAACCATTGGGAAACCAAGGTTTTTGAGATCGAAACCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACGGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTTGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGTCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTAC‑ATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTgc                                                                                                                                                                                                               <  1:262353/249‑1 (MQ=255)                                         aCCATGGGGAAACCAAGGTTTTTGAGATCGAAACCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGCACACAGGGCGGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGAGTGATCGTCGCGGTTGCCTGGTCTGTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGAGCCACATAGC‑TAGGATCTGGCTGATGCTTGGTCCAAACTGGAAGGACGGGCGCTatat                                                                                                                                                                                                        <  2:283385/251‑1 (MQ=255)                                                ggAAACCAAGGTTTTTGAGATCGAAACCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGACCAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTAC‑ATCTGGATGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAAc                                                                                                                                                                                                 >  2:581220/1‑251 (MQ=255)                                                  aaaCCAAGGTTTTTTAGATCGAAACCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTAC‑ATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAAcat                                                                                                                                                                                               <  1:481803/251‑1 (MQ=255)                                                              ttGAGATCGAAACCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTAC‑ATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCt                                                                                                                                                                                   >  1:229837/1‑251 (MQ=255)                                                                gagaTCGAAACCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGTTAGATCCACATAGCTTAC‑ATCTGGCTGATGCTTTGCCCAATCTGTAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTa                                                                                                                                                                                 <  2:459984/251‑1 (MQ=255)                                                                  gagcGAAAACAAACGTGAAGGATTTGCTTTCTTGTATCGTAATCCACAGTACACAGGGCAGTCGCCACTGGGAATCGCTTTTGTTTAAGCTGAGCAGTACAAGATTCTCCGTCTCTATTTTCTGTTCTTTGCCGAACAGTATGGCTAAATTCTTGTGGACCTGGTAGATCCACATAGCTTAC‑ATCTGGCTGATTCTTTGCCCAATCTTGAAGGTCTTGCGCCATATGCCTAACATCTTTCTGAGTCTTACa                                                                                                                                                                               <  2:422442/248‑1 (MQ=255)                                                                      gAAACCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTAC‑ATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATCCTGATGCTTACAGGcg                                                                                                                                                                           <  1:138382/251‑1 (MQ=255)                                                                               cGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTAC‑ATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGCgaatcgaatc                                                                                                                                                                  >  2:311151/1‑251 (MQ=255)                                                                                 tGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTCCACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTAC‑ATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGCGAATCGAATCTg                                                                                                                                                                >  1:182400/1‑251 (MQ=255)                                                                                 tGAAGGATTTGCTTTCTGGAATCGTAATCAACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTAC‑ATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGCGAATCGAATCTg                                                                                                                                                                >  1:489317/1‑251 (MQ=255)                                                                                  gAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTAC‑ATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGCGAATCGAATCTGt                                                                                                                                                               >  1:496803/1‑251 (MQ=255)                                                                                          tGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGTTTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGACCCACATAGCTTAC‑ATCTGGCTGATTCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTAAAGGCGAATCGAATCTGTCGCCGaa                                                                                                                                                        <  1:607322/250‑1 (MQ=255)                                                                                                  ggTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTAC‑ATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCGGATGCTTACAGGCGAATCGAATCTGTCGCCGAAGTAAAgg                                                                                                                                                 <  2:43341/249‑1 (MQ=255)                                                                                                   gTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTAC‑ATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCACTCTGATGCTTACAGGCGAATCGAATCTGTCGCCGAAGTAAAGGGt                                                                                                                                               >  1:8999/1‑250 (MQ=255)                                                                                                       cGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTAC‑ATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGCGAATCGAATCTGTCGCCGAAGTAAAGGGTAAAtt                                                                                                                                          <  1:515285/251‑1 (MQ=255)                                                                                                               acaGTACCCAGGGCAGTCGACACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCACAATGGTTGTGGACCTGGTAGATCCACATAGCTTAC‑ATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGCGAATCGAATCTGTCGCCGAAGTAAAGGGAAAATTACGAgtgt                                                                                                                                  >  2:329105/1‑251 (MQ=255)                                                                                                                       caGGGCAGTCGTCACTGGGACTCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTAC‑ATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGCGAATCGAATCTGTCGCCGAAGTAAAGGGTAAATTACGAGTGTTAGATTTg                                                                                                                          >  1:480391/1‑251 (MQ=255)                                                                                                                                 gtcACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTAC‑ATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGCGAATCGAATCTGTCGCCGAAGTAAAGGGTAAATTACGAGTGTTAGATTTGAAACGGCAgg                                                                                                                >  1:458205/1‑251 (MQ=255)                                                                                                                                   cACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATGTCCTGTTCTGTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTAC‑ATCTGGCTGATGCTTTGCCCAAACGGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCGTACAGGCGAATCGAATCTGTCGCCGAAGGAAAGTGTAAATGACGAGTGTTAGAGTCGAAACGGCAGGAt                                                                                                              <  2:623502/251‑1 (MQ=255)                                                                                                                                    aCTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTAC‑ATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGCGAATCGAATCTGTCGCCGAAGTAAAGGGTAAATTACGAGTGTTAGATTTGAAACGGCAGGAtg                                                                                                             <  2:469553/251‑1 (MQ=255)                                                                                                                                                    gTTGACGCTGAACACTACCCGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAATTGGTTGTGGCCCTGGTAGATCCACATAGCTTAC‑ATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATCTGCCGAACATCATTCTGATGCTTACAGGCGAATCGACTCTGTCGCCGACGTAAAGGGTAAATTACGAGTGTTAGATTTGAAACGGCAGGATGTGCAGGAAGCTGTTGc                                                                                             >  2:266435/1‑251 (MQ=255)                                                                                                                                                      tGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGCGGACCTGGTAGATCCACATAGCTTAC‑ATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGCGAATCGAATCTGTCGCCGAAGTAAAGGGTAAATTACGAGTGTTAGATTTGAAACGGCAGGATGTGCAGGCTGCTGTTGCt                                                                                            >  1:624368/1‑250 (MQ=255)                                                                                                                                                                         aTTGTTCGTTCCGATTTCCTGGTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTAC‑ATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGCGAATCGAATCTGTCGCCGAAGTAAAGGGTAAATTACGAGTGTTAGATTTGAAACGGCAGGATGTGCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAAc                                                                        <  2:341077/251‑1 (MQ=255)                                                                                                                                                                               cGTCCCGATTTCCTGTTCTTTGCCGAACAGGTTGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTAC‑ATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACCGGCGAATCGAATCTGTCGCCGAAGTAAAGGGTAAATTACGAGTGTTAGATTTGAACCGGAAGGATGTGCAGGATGCTGTTGCTACCGCTGAAAATGACGAAAAGTTAtt                                                                  >  2:481825/1‑251 (MQ=255)                                                                                                                                                                                gtaccGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTAC‑ATCTGGCTGATGCTTTGCCCAAACTGGCAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGCGAATCGAATCTGTAGCCGAAGTAAAGGGTAAATTACGAGTGTTAGATTTGAAACGGCAGGATGTGCAGGAAGCTGTTGATACCGCTGAAAATGCAGAAAAGTTAttt                                                                 >  2:419049/4‑251 (MQ=255)                                                                                                                                                                                 tCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTAC‑ATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGCGAATCGAATCTGTCGCCGAAGTAAAGGGTAAATTACGAGTGTTAGATTTGAAACGGCAGGATGTGCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTAttt                                                                 >  1:475207/1‑250 (MQ=255)                                                                                                                                                                                    cGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTAC‑ATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGCGAATCGAATCTGTCGCCGAAGTAAAGGGTAAATTACGAGTGTTAGATTTGAAACGGCAGGATGTGCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACTTTATTtagta                                                             <  1:188743/251‑1 (MQ=255)                                                                                                                                                                                      aTTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTAC‑ATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGCGAATCGAATCTGTCGCCGAAGTAAAGGGTAAATTACGAGTGTTAGATTTGAAACGGCAGGATGTGCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGc                                                           >  1:461542/1‑251 (MQ=255)                                                                                                                                                                                      aTTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTAC‑ATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGCGAATCGAATCTGTCGCCGAAGTAAAGGGTAAATTACGAGTGTTAGATTTGAAACGGCAGGATGTGCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGc                                                           <  1:456553/251‑1 (MQ=255)                                                                                                                                                                                      aTTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTAC‑ATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGCGAATCGAATCTGTCGCCGAAGTAAAGGGTAAATTACGAGTGTTAGATTTGAAACGGCAGGATGTGCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGc                                                           <  1:581220/251‑1 (MQ=255)                                                                                                                                                                                          ccTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTAC‑ATCTGGCTGATGCGTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGCGAATCGAATCGGTCGCCGAAGTAAAGGGTAAATTACGAGTGTTAGATTTGAAACGGCAGGATGTGCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGAc                                                       <  2:172912/251‑1 (MQ=255)                                                                                                                                                                                                      cGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCTTAC‑ATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGCGAATCGAATCTGTCGCCGAAGTAAAGGGTAAATTACGAGTGTTAGATTTGAAACGGCAGGATGTGCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACt                                           >  1:73373/1‑251 (MQ=255)                                                                                                                                                                                                           aGGATGGCAACATGGTTGTGGACCTGGTAGATCCACATAGCTTAC‑ATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGCGAATCGAATCTGTCGCCGAAGTAAAGGGTAAATTACGAGTGTTAGATTTGAAACGGCAGGATGTGCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATc                                        >  2:493277/1‑249 (MQ=255)                                                                                                                                                                                                                           tgtgGACCTGGTAGATCCACATAGCTTAC‑ATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGCGAATCGAATCTGTCGCCGAAGTAAAGGGTAAATTACGAGTGTTAGATTTGAAACGGCAGGATGTGCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAAtt                      >  2:317306/1‑251 (MQ=255)                                                                                                                                                                                                                            gtgGACCTGGTAGATCCACATAGCTTAC‑ATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGCACATCATTCTGGTGCTTACAGGCGAATCGAATCTGTCGCCGAAGTAAAGGGTAAATTACGAGTGTTAGATTTGAAAAGACAGGAAGTGCAGGATGCTGTTGCTACCGATGAAAATGCAGAAACGTTATTTAGTAGCGGTCTTGCTGGTGACTATCAGTAATCTATAGAAAATTg                     >  1:629117/1‑251 (MQ=255)                                                                                                                                                                                                                            gtgGACCTGGTAGATCCACATAGCTTAC‑ATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTTATGCTTACAGGCGAATCGAATCTGTCGCCGAAGTAAAGGGTAAATTACGAGTGTTAGATTTGAAACGGCAGGATGTGCAGGATGCTGTTGCTACCGTTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAAttt                     <  2:469666/251‑2 (MQ=255)                                                                                                                                                                                                                              ggACCTGGTAGATCCACATAGCTTAC‑ATCTGGCTGTTGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGCGAATCGAATCTGTCGCCGAAGTAAAGGGTAAATTACGAGTGTTAGATTTGAAACGGCAGGATGTGCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGCCTTGCTGATGACTATCAGTACTCTATAGCAAATTGCg                   >  2:538653/1‑251 (MQ=255)                                                                                                                                                                                                                                    ggTAGATCCACATAGCTTAC‑ATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGCGAATCGAATCTGTCGCCGAAGTAAAGGGTAAATTACGAGTGTTAGATTTGAAACGGCAGGATGTGCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTAGTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACg              <  2:45776/250‑1 (MQ=255)                                                                                                                                                                                                                                     gTAGATGCACCTAGCTTAC‑ATCTGGCTGATGCTTCGCCCAAACTGGAAGGACTTCCGCTATATGCCGAAACTAATTCTGATGCTTAAAGGCGAATCGAATCTGTCGCCGAGGAAAAGGGTAAATTACGAGTGTTAGATTTGAAACGGCAGGATGTGCAGGATGCTGTTGCTACCGCTGCAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGCTGACTATCAGTAATCTCTAGAAAATTGCGTAACgc             >  2:599713/1‑249 (MQ=255)                                                                                                                                                                                                                                        gATCCACATAGCTTAC‑ATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGCGAATCGAATCTGTCGCCGAAGTAAAGGGTAAATTACGAGTGTTAGATTTGAAACGGCAGGATGTGCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACt         >  1:4098/1‑251 (MQ=255)                                                                                                                                                                                                                                           ccACATAGCTTAC‑CTCTGGCTGATGCTTTGCCCAACCTGGAAGGACTTGCGCTATATGCCGACCATCATTCTGATGCTTACAGGCGAATCGCATCTGTCGCCGCAGTCAAGGGTAAATTACGAGTGTTAGATTTGAACCGGCAGGATGTGCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGc       >  2:268619/1‑250 (MQ=255)                                                                                                                                                                                                                                           ccACATAGCTTAC‑ATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGCGAATCGAATCTGTCGCCGAAGTAAAGGGTAAATTACGAGTGTTAGATTTGAAACGGCAGGATGTGCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGc       >  1:32434/1‑250 (MQ=255)                                                                                                                                                                                                                                            cacaTAGCTTAC‑ATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGTCGAATCGAATCTGTCGCCGAAGTAAAGGGTAAATTACGAGTGTTAGATTTGAAACGGCAGGATGTGCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCAt     <  2:240842/251‑1 (MQ=255)                                                                                                                                                                                                                                               aTAGCTTAC‑ATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGCGAATCGAATCTGTCGCCGAAGTAAAGGGTAAATTACGAGTGTTAGATTTGAAACGGCAGGATGTGCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAAGTGCGGAACGAAACTGCAtttt  <  2:437081/251‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       TTGTTGTGTGACGAAAGCGGGAACTATCCGGCTCTGTTGAACCATTGGGAAACCAAGGTTTTTGAGATCGAAACCAAACGTGAAGGATTTGCTTTCTGGTATCGTAATCCACAGTACACAGGGCAGTCGTCACTGGGAATCGCTTATGTTGAAGCTGAACAGTACAAGATTGTTCGTCCCGATTTCCTGTTCTTTGCCGAACAGGATGGCAAAATGGTTGTGGACCTGGTAGATCCACATAGCCTAC‑ATCTGGCTGATGCTTTGCCCAAACTGGAAGGACTTGCGCTATATGCCGAACATCATTCTGATGCTTACAGGCGAATCGAATCTGTCGCCGAAGTAAAGGGTAAATTACGAGTGTTAGATTTGAAACGGCAGGATGTGCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTT  >  NC_000913/291776‑292264

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