Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	46,655	A→G	L283L (CTA→CTG)	yaaU →	putative transporter YaaU

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	46,655	0	A	G	98.2%	177.5 / NA	57	L283L (CTA→CTG)	yaaU	putative transporter YaaU
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (29/27);  minor base T (0/1);  total (29/28)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

AAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATCAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGGGTTGGCAAAAACGCGGCACTAGGGAATGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTAATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTG‑TATG‑CCTTTTTCTCTGGCGGGCCGGGTAATTTGCAGTGGCTCTATCCTAATGA  >  NC_000913/46413‑46888                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                            aaggggCGAGTAAAAGAGGGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGaacgt                                                                                                                                                                                                                                     <  2:484239/248‑4 (MQ=255)      gctaGTAAAAGAGTGCGAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTCCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGa                                                                                                                                                                                                                                <  1:224068/248‑1 (MQ=255)                     cGAAGAGATGATGATAAAACTGTTTGGCGAACCGTTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCACTGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCACAAACGCGGCGCTGGGGAACGTGATGATTAGCCTGTTCttt                                                                                                                                                                                                                  >  1:196033/1‑250 (MQ=255)                     cGAAAAGATGATGAAAAAACCGTTTTGCGAACCGGTGGCTTTCGATGAAAAGCGGCAGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCCTTTTCCTTTTGTTCTGTTTGTTGACGCCATCTGGACCTGCCATGTGATCCCAATGTTCTCCATTTACACCTTTGGCCCGCAAATCGTTAGTTTGTTGGGATTGAGTGTTGGCAAAAACGCGGCGCTGGGGTACGTGGTGAGTAGCCTGTTGttt                                                                                                                                                                                                                  >  2:632653/1‑250 (MQ=255)                      gAAGAGATGATGATAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAATCGCGGCGCTGGGGCACGTGGTGATTAGCCTGTTCTTTAt                                                                                                                                                                                                                >  2:341949/1‑251 (MQ=255)                            atgatgatAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCAATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGACTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCgg                                                                                                                                                                                                          <  1:4070/251‑1 (MQ=255)                               atgatAAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTg                                                                                                                                                                                                       >  1:90872/1‑251 (MQ=255)                                tgataaaccTGTTTGGCGAACCGGTGGCTATCGATCAAGAGCAGCCTCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTGTCCTTTTGATCCGCTTGTTGCCGACATCTGGACCTGCCAGGCGATCCTACTGCTCGCCATTTACACCCTTGGACCGAAAATCGTTGGTTTGTTGGGATTGGGCATTAGCAGAAAAGCGGAGTTGGGGTACGTGGGGTTTAGCCTGTTATTTATGCTCGGCTGt                                                                                                                                                                                                      >  2:515448/9‑251 (MQ=255)                                     aaaCTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTcc                                                                                                                                                                                                 >  2:118849/1‑251 (MQ=255)                                      aaCTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTccg                                                                                                                                                                                                <  1:329949/251‑1 (MQ=255)                                          gTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTccgccg                                                                                                                                                                                             <  1:218096/250‑1 (MQ=255)                                                         ggtcttctatgaatagacgcCGCAGCAAACCCTTTTGCGCGATCTCTTTATTCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCAGTTACACCTTTGGCCCGCAAATCGTTGGGTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTTGGGACCGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAaca                                                                                                                                                                             <  2:622451/233‑1 (MQ=255)                                                                      gagCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCt                                                                                                                                                                                                                         >  1:3976/1‑194 (MQ=255)                                                                          agcCGCCGCAAACCCGTTTTCGCGATCTGTTTAATAGCCGCCCTTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCCTTTCCACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGACCGTGGTGATTAg                                                                                                                                                                                                                            >  2:402800/1‑187 (MQ=255)                                                                          agcCGCAGCAAACCCGTTTTTGCGATCTGTTTAATCGCCGCCAGTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCCGCCAGGTGATCCCAATGTTCGCCATTTACCCCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAg                                                                                                                                                                                                                            <  1:402800/187‑1 (MQ=255)                                                                              gcagcaAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCAttgt                                                                                                                                                        >  1:237571/1‑251 (MQ=255)                                                                                    aaCCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTg                                                                                                                                                  >  1:243067/1‑251 (MQ=255)                                                                                    aaCCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGCGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAATAACTGCCGGACGGCGTCCATTGTTGATTg                                                                                                                                                  >  1:635782/1‑251 (MQ=255)                                                                                           tttCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTTTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCAGTGTTGATTGGCAGCt                                                                                                                                            <  1:580026/250‑1 (MQ=255)                                                                                                     tGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGACATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGGCCATTGTTGATTGGCAGCTTTGCCatgatg                                                                                                                                 >  2:584196/1‑251 (MQ=255)                                                                                                      gTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCat                                                                                                                                     <  1:178797/246‑1 (MQ=255)                                                                                                          aaTCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGCCCTGCCAGGTGAGCCCAATGTTCGCCATTTACACCTTGGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGAcgct                                                                                                                            <  1:616268/251‑1 (MQ=255)                                                                                                             cgccgcCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAacac                                                                                                                                                                            >  2:155807/1‑200 (MQ=255)                                                                                                             cgccgcCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAacac                                                                                                                                                                            <  1:155807/200‑1 (MQ=255)                                                                                                             cgccgcCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGAcgctggc                                                                                                                         >  1:377468/1‑251 (MQ=255)                                                                                                                 gccATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAACTCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGCCGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGAcgctggctctg                                                                                                                     >  2:565165/1‑251 (MQ=255)                                                                                                                      tttCCTTTTGTTCTGTTTGTTGCCGCCATCTGGGCCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAAGCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCgg                                                                                                                 <  1:142829/250‑1 (MQ=255)                                                                                                                       ttCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCgggt                                                                                                               >  1:460681/1‑249 (MQ=255)                                                                                                                          cTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTTCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGTCGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTg                                                                                                            >  1:182377/1‑251 (MQ=255)                                                                                                                          cTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTGTTg                                                                                                            >  1:244881/1‑251 (MQ=255)                                                                                                                               gTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTTGGGCt                                                                                                       >  1:209545/1‑251 (MQ=255)                                                                                                                                    gtttgttGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATcc                                                                                                  >  1:436444/1‑251 (MQ=255)                                                                                                                                    gtttgtgGCCCCCATCTGGACCTGCCAGTGGATCCCAATGTTCGCCATTTACTCCTGTGGCCCGCAAATCGTTGTGTTGTTGGGATTGGTTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTGCTTTAGGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGTCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATcc                                                                                                  <  2:237571/251‑1 (MQ=255)                                                                                                                                              gccATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGGGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGTTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCCAGATGACGCTGGCGATGACGGTTTTtggggtggttcacgtgatggg                                                                                        >  2:367092/1‑234 (MQ=255)                                                                                                                                               ccATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGGTCCCGGGTATgggg                                                                                       >  1:404272/1‑251 (MQ=255)                                                                                                                                                   cTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCACATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTACACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTgggggttgtgccgggtttggggatct                                                                                   >  2:555232/1‑229 (MQ=255)                                                                                                                                                    tggcccTGCCAGGTGATCCCAATGTTTGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTGGGTATTGGTTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGAt                                                                                                                                                                                           <  2:22434/142‑1 (MQ=255)                                                                                                                                                    tGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGAt                                                                                                                                                                                           >  1:22434/1‑146 (MQ=255)                                                                                                                                                    gggACCTGCCAGGGGATCCCAAGGTTCGCCATTTACGCCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATCCCGGATATGGGGATctg                                                                                  <  2:90185/250‑1 (MQ=255)                                                                                                                                                     ggACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATCCCGGATATGGGGATct                                                                                   <  1:118849/249‑1 (MQ=255)                                                                                                                                                      gACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGGATTCCGCCGATGCTGTGGTTCAACACTGCCGGAAGGCGTCCATTGTTGATTGGCAGCTTTGCCATGGTGAATCCGGAGCTGGCGGTTTTGGGGCtgggccacgacattgggggttggt                                                                                >  2:22015/1‑228 (MQ=255)                                                                                                                                                        ccTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATCCCGGATATGGGGATctggctg                                                                              <  1:493978/251‑1 (MQ=255)                                                                                                                                                               ggTGATCCGAATGTGCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATTCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGAGTGGCAGCTTTGCCATGATTACGCTGGCGCTGGCGGTTTTGGGGCTGATCCCGGATATGGGGATCTGGCTGGTAGTGa                                                                       <  1:646487/251‑1 (MQ=255)                                                                                                                                                                      ccAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCtt                                                                <  1:112850/251‑1 (MQ=255)                                                                                                                                                                          tGTTCGCCATTTACCCCTTTGGCCCGCAAATCGTTGGTTTGGTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGAcgctggc                                                                                                                         <  2:557159/190‑1 (MQ=255)                                                                                                                                                                          tGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGAcgctggc                                                                                                                         >  1:557159/1‑190 (MQ=255)                                                                                                                                                                           gTTCTCCAGTTACCCCTTTGGCCCGCAAATCGGTGGTTTTGGGGGATTGGGGGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCCCGGCTGTATTCCGCCGATGCTGTGGTTAAACCCTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATTACGCGGGCGCTGGCTGTTTTGGGGCTGATCCCGGATATGGGGATCTGGCTGGTAGTGATGTGCTTTGct                                                            <  2:218386/250‑2 (MQ=255)                                                                                                                                                                                 cATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTTCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTG‑TATg                                                    <  1:640553/251‑1 (MQ=255)                                                                                                                                                                                  aTTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTG‑TATG‑c                                                  <  1:555232/251‑1 (MQ=255)                                                                                                                                                                                      acacCTTTGGCCCGCAAATCGTTGGGTTGTTTGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTG‑TATG‑CCtt                                               <  2:171072/250‑1 (MQ=255)                                                                                                                                                                                                cccGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGTCGCTGGGGAACGTGGTGATGAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGTTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTG‑TATG‑CCTTTTTCTCTGGCg                                    <  1:410519/251‑1 (MQ=255)                                                                                                                                                                                                    aaaaTCGTGGGGTTGATGGGATGGGGTGTTGGCAAAAACGCGGCGCTGGGAAACGTGGTGATTAGCCTGTTCTTTATGCGCGGCTGTATGCCGCCGATGCTGTGGTTAAACACTGCCTGACGTCGTCCATTGTTGATTGGCAGCTGTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCGTATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTG‑TATG‑CCTTTGTCTCTGGCGGGcc                                <  2:644146/250‑1 (MQ=255)                                                                                                                                                                                                                 tgttgGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTG‑TATG‑CCTTTTTCTCTGGCGGGCCGGGTAATTTGCa                    <  1:565165/250‑1 (MQ=255)                                                                                                                                                                                                                  gttgGGATTGGGAGTTGGCAAAAACGCGGCGCTGGGGAACGCGGTGATTAGCCTGTTCTTTATGCTCGGATGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCTATGATGACGCTGGCGCTGGCGGCTTTGGGGCTGATACCGGGTATGGGGATCTCGCTGGGTGGGGTGGCCTTTGCGGAG‑TATGGCCTTTTT‑TCTCGCGGGACGGGGAATTTGCAgg                  >  1:43415/1‑250 (MQ=255)                                                                                                                                                                                                                        aTTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGGATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGGTGATCCCGGGTTTGGGGGACTGGCTGGGAGGGATGGCGCTTGCGGGGGTTTG‑C‑TTTTTTTCTGGCGGGCCgggggattttgcatgggctc            >  2:19161/1‑234 (MQ=255)                                                                                                                                                                                                                           gggtgTTGGCAAAAACGCGGCGCTGGGGATCGGGTTGATTAGCTTGTTCTTTATGCTCGGTTGTATTCCTCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACTCTGGCGCTGGCGGTTTTGGGGCTGATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTG‑TATG‑CCTTTTTCTCTGGCGTGCCGGGTAATTTGCAGTGGCTCTATc         <  2:560773/247‑1 (MQ=255)                                                                                                                                                                                                                                   gCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGGTGATCCCGGATATGGGGATCTGGCTGGAAGTGATGGCCTTTGCGGGG‑TATG‑CCTTTTTCTCTGGCGGGCCGGGGGATTTGCAGTGGCTCTATCCTAATgg  >  1:465329/1‑249 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                            AAAGGGCGAGTAAAAGAGTGCGAAGAGATGATGATCAAACTGTTTGGCGAACCGGTGGCTTTCGATGAAGAGCAGCCGCAGCAAACCCGTTTTCGCGATCTGTTTAATCGCCGCCATTTTCCTTTTGTTCTGTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGGGTTGGCAAAAACGCGGCACTAGGGAATGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTAATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTG‑TATG‑CCTTTTTCTCTGGCGGGCCGGGTAATTTGCAGTGGCTCTATCCTAATGA  >  NC_000913/46413‑46888

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