Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,070,014	A→G	L114P (CTG→CCG)	rutF ←	FMN reductase RutF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,070,014	0	A	G	97.5%	133.4 / NA	40	L114P (CTG→CCG)	rutF	FMN reductase RutF
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (23/16);  total (24/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

CGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCAG‑AGCCTCTTCCAGTTGCGGGCATCCGGTCACCCCCGTCTGCCAGCGGGCGGCGGCAAAGCGGTGTTCCATGGGCGTTTTGCCGCCAAAAAGGTTTGAAAGCGGCTCCTGCCCGGCGCTAAGTGTATTTACACACAGCGTTCGATTTTCATTGAATGCCGGCCAGACGGACGCCCCACGATTCAGGCACACCAGTAATGTGGGCGGCGTATCGGTCACACTGCAGACGGC  >  NC_000913/1069771‑1070243                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                       gggccgtaCCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCGG‑AGCctc                                                                                                                                                                                                                                <  2:294271/244‑1 (MQ=255) cGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCGG‑AGCctc                                                                                                                                                                                                                                <  2:59268/251‑1 (MQ=255)          cACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGGCAAACGAAACCGG‑AGCCTCTTCCAGTTg                                                                                                                                                                                                                       <  2:260104/251‑1 (MQ=255)             gccgcTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCGG‑AGCCTCTTCCAGTTGCgg                                                                                                                                                                                                                    >  2:112366/1‑251 (MQ=255)             gccgcTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCGG‑AGCCTCTTCCAGTTGCgg                                                                                                                                                                                                                    >  1:66353/1‑251 (MQ=255)                             gTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCGG‑AGCCTCTTCCAGTTGCGGGCATCCGGTCAccccc                                                                                                                                                                                                    >  2:210089/1‑251 (MQ=255)                                            cAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCGG‑AGCCTCTTCCAGTTGCGGGCATCCGGTCACCCCCGTCTGCCAgcgggcg                                                                                                                                                                                     >  1:179865/1‑251 (MQ=255)                                               tGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGAGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAACCAGAATGTCGTGGGTGCAGACGCTCACCACCTGGCTGATACGGCAGTCAAACGCAACCGG‑AGCCTCTTCAAGTTGCGGGCATCGGGTCACCCCCGTCTGCCAGCGggcggcg                                                                                                                                                                                  >  1:194680/1‑251 (MQ=255)                                                                  ccATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCGG‑AGCCTCTTCCAGTTGCGGGCATCCGGTCACCCCCGTCTGCCAGCGGGCGGCGGCAAAGCGGTGTTCCATgg                                                                                                                                                               <  2:210507/251‑1 (MQ=255)                                                                  ccATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGCGCGAATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCGG‑AGCCTCTTCCAGTTGCGGGCATCCGGTCACCCCCGTCTGCCAGCGGGCGGCGGCAAAGCGGCGTTCCATgg                                                                                                                                                               >  1:230978/1‑251 (MQ=255)                                                                      cAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCGG‑AGCCTCTTCCAGTTGCGGGCATCCGGTCACCCCCGTCTGCCAGCGGGCGGCGGCAAAGCGGTGTTCCATGGGCGt                                                                                                                                                           <  1:259406/251‑1 (MQ=255)                                                                                   aaGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCGG‑AGCCTCTTCCAGTTGCGGGCATCCGGTCACCCCCGTCTGCCAGCGGGCGGCGGCAAAGCGGTGTTCCATGGGCGTTTTGCCGCCaaa                                                                                                                                               >  1:363565/1‑250 (MQ=255)                                                                                                 gCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCGG‑AGCCTCTTCCAGTTGCGGGCATCCGGTCACCCCCGTCTGCCAGCGGGCGGCGGCAAAGCGGTGTTCCATGGGCGTTTTGCCGCCAAAAAGGTTTGAAAGCgg                                                                                                                                <  1:308765/251‑1 (MQ=255)                                                                                                    gggCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCGG‑AGCCTCTTCCAGTTGCGGGCATCCGGTCACCCCCGTCTGCCAGCGGGCGGCGGCAAAGCGGTGTTCCATGGGCGTTTTGCCGCCAAAAAGGTTTGAAAGCGGCTc                                                                                                                             <  1:132378/251‑1 (MQ=255)                                                                                                     ggCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCGG‑AGCCTCTTCCAGTTGCGGGCATCCGGTCACCCCCGTCTGCCAGCGGGCGGCGGCAAAGCGGTGTTCCATGGGCGTTTTGCCGCCAAAAAGGTTTGAAAGCGGCTcc                                                                                                                            >  1:65551/1‑251 (MQ=255)                                                                                                          aTCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCGG‑AGCCTCTTCCAGTTGCGGGCATCCGGTCACCCCCGTCTGCCAGCGGGCGGCGGCAAAGCGGTGTTCCATTGGCGTTTTGCCGCCAAAAAGGTTTGAAAGCGGCTCCTGccc                                                                                                                       >  1:282825/1‑251 (MQ=255)                                                                                                                            ggATCGAGCCACACCAGCCCGGAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGTTGCCGACGCTCACCCCCTGGCTGATACGGCAGTCAAACGAAACCGG‑AGCCTCTTCCAGTTGCGGGCATCCGGTCACCCCCGTCTGCCAGCGGGCGGCGGCAAAGCGGTGTTCCATGGGCGTTTTGCCGCAAAAAAGGTTTGAAAGCGGCTCCTGCCCGGCGCTAAGTGTATTTac                                                                                                     <  2:277038/250‑1 (MQ=255)                                                                                                                                    ccACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCGG‑AGCCTCTTCCAGTTGCGGGCATCCGGTCACCCCCGTCTGCCAGCGGGCGGCGGCAAAGCGGTGTTCCATGGGCGTTTTGCCGCCAAAAAGGTTTGAAAGCGGCTCCTGCCCGGCGCTAAGTGTATTTACACACAGCg                                                                                             >  1:288496/1‑251 (MQ=255)                                                                                                                                       cacCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCTACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCGG‑AGCCTCTTCCAGTTGCGGGCATCCGGTCACCCCCGTCTGCCAGCGGGCGGCGGCAAAGCGGTGTTCCATGGGCGTTTTGCCGCCAAAAAGGTTTGAAAGCGGCTCCTGCCCGGCGCTAAGTGTATTTACACACAGCGTTc                                                                                          <  1:180015/251‑1 (MQ=255)                                                                                                                                        acCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCGG‑AGCCTCTTCCAGTTGCGGGCATCCGGTCACCCCCGTCTGCCAGCGGGCGGCGGCAAAGCGGTGTTCCATGGGCGTTTTGCCGCCAAAAAGGTTTGAAAGCGGCTCCTGCCCGGCGCTAAGTGTATTTACACACAGCGTTCg                                                                                         >  1:125787/1‑251 (MQ=255)                                                                                                                                         ccAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCGG‑AGCCTCTTCCAGTTGCGGGCATCCGGTCACCCCCGTCTGCCAGCGGGCGGCGGCAAAGCGGTGTTCCATGGGCGGTTTGCCGCCAAAAAGGGTTGAAAGCGGCCCCTGCCCGGCGCTAAGCGTATTTACACAACGCGTTCg                                                                                         >  2:226914/1‑250 (MQ=255)                                                                                                                                                    ggtgtggtgtGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCGG‑AGCCTCTTCCAGTTGCGGGCATCCGGTCACCCCCGTCTGCCAGCGGGCGGCGGCAAAGCGGTGTTCCATGGGCGTTTTGCCGCCAAAAAGGTTTGAAAGCGGCTCCTGCCCGGCGCTCAGTGTATTTACACACAGCGTTCGATTTTCATTGa                                                                              >  1:85190/1‑250 (MQ=255)                                                                                                                                                      tgtggtgtgACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCGG‑AGCCTCTTCCAGTTGCGGGCATCCGGTCACCCCCGTCTGCCAGCGGGCGGCGGCAAAGCGGTGTTCCATGGGCGTTTTGCCGCCAAAAAGGTTTGAAAGCGGCTCCTGCCCGGCGCTAAGTGTATTTACACACAGCGTTCGATTTTCATTGAATg                                                                           >  1:286506/1‑251 (MQ=255)                                                                                                                                                      ggtggtgtgACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGGGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGGCAAACGAAACCGG‑AGCCTCTTCCAGTTGCGGGCATCCGGTCACCCCCGTCTGCCAGCGGGCGGCGGCAAAGCGGTGTGCCATGGGCGTTTTGCCGCCAAAAAGGTTTGAAAGCGGCTCCGGCCCGGCGCTAAGTGTATTTACACACAGCGTTCGATTTTCATTGAATg                                                                           <  2:194680/250‑1 (MQ=255)                                                                                                                                                                cGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCGG‑AGCCTCTTCCAGTTGCGGGCATCCGGTCACCCCCGACTGCCAGCGGGCGGCGGCAAAGCGGTGTTCCATGGGCGTTTTGCCGCCAAAAAGGTTTGAAAGCGGCTCCTGCCCGTCGCTAAGTGTATTTACACACATCGTTCGATTTTCATTGAATGCCGGCCAgac                                                                 <  1:330537/251‑1 (MQ=255)                                                                                                                                                                                        aaCAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCGG‑AGCCTCTTCCAGTTGCGGGCATCCGGTCACCCCCGTCTGCCAGCGGGCGGCGGCAAAGCGGTGTTCCATGGGCGTTTTGCCGCCAAAAAGGTTTGAAAGCGGCTCCTGCCCGGCGCTAAGTGTATTTACACACAGCGTTCGATTTTCATTGAATGCCGGCCAGACGGACGCCCCACGATTCAGGCacac                                         <  1:7703/251‑1 (MQ=255)                                                                                                                                                                                         aCAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCCAACGAAACCGG‑AGCCTCTTCCAGTTGCGGGCATCCGGTCACCCCCGTCTGCCAGCGGGCGGCGGCAAAGCGGTGTTCCATGGGCGTTTTGCCGCCAAAAAGGTTTGAAAGCGGCTCCTGCCCGGCGCTAAGTGTATTTACACACAGCGTTCGATTTTCATTGAATGCCGGCCAGACGGACGCCCCACGATTCAGGCACAcc                                        >  2:26229/1‑251 (MQ=255)                                                                                                                                                                                          cAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCGG‑AGCCTCTTCCAGTTGCGGGCATCCGGTCACCCCCGTCTGCCAGCGGGCGGCGGCAAAGCGGTGTTCCATGGGCGTTTTGCCGCCAAAAAGGTTTGAAAGCGGCTCCTGCCCGGCGCTAAGTGTATTTACACACAGCGTTCGATTTTCATTGAATGCCGGCCAGACGGACGCCCCACGATTCAGGCACACCa                                       <  1:386039/251‑1 (MQ=255)                                                                                                                                                                                             aaTGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCGG‑AGCCTCTTCCAGTTGCGGGCATCCGGTCACCCCCGTCTGCCAGCGGGCGGCGGCAAAGCGGTGTTCCATGGGCGTTTTGCCGCCAAAAAGGTTTGAAAGCGGCTCCTGCCCGGCGCTAAGTGTATTTACACACAGCGTTCGATTTTCATTGAATGCCGGCCAGACGGACGCCCCCCGATTCAGGCACACCAGt                                     >  2:151402/1‑250 (MQ=255)                                                                                                                                                                                               tGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCCAACGAAACCGG‑AGCCTCTTCCAGTTGCGGGCATCCGGTCACCCCCGTCTGCCAGCGGGCGGCGGCAAAGCGGTGTTCCATGGGCGTTTTGCCGCCAAAAAGGTTTGAAAGCGGCTCCTGCCCGGCGCTAAGTGTATTTACACACAGCGTTCGATTTTCATTGAATGCCGGCCAGACGGACGCCCCACGATTCAGGCACACCAGTAAt                                  >  2:267992/1‑251 (MQ=255)                                                                                                                                                                                                gTCGTGGGCGCCGACGCTCACAACCTGGCTGATACGGCAGTCAAACGCAAC‑AGCAGCCTCTTCCAGTTGCGGGCATCCGGCCACCCCCGTCTGCCAGCGGGCGGCGGCAAAGCGGTGCTCCATGGGCGTTTTGCCGCCAAAAAGGTTTGAAAGCGGCTCCTCCCCGGCGCTAAGTGTATTTACACACAGCGTTCGATTTTCATTGAATGCCGGCCAGACCGACGCCCCACGATTCCGGCACACCAGTAAtt                                 >  2:407836/1‑250 (MQ=255)                                                                                                                                                                                                   gTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCGG‑AGCCTCTTCCAGTTGCGGGCATCCGGTCACCCCCGTCTGCCAGCGGGCGGCGGCAAAGCGGTGTTCCATGGGCGTTTTGCCGCCAAAAAGGTTTGAAAGCGGCTCCTGCCCGGCGCTAAGTGTATTTACACACAGCGTTCGATTTTCATTGAATGCCGGCCAGACGGACGCCCCACGATTCAGGCACACCAGTAATGTgg                              >  2:80025/1‑251 (MQ=255)                                                                                                                                                                                                    tGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCGG‑AGCCTCTTCCAGTTGCGGGCATCCGGTCACCCCCGTCTGCCAGCGGGCGGCGGCAAAGCGGTGTTCCATGGGCGTTTTGCCGCCAAAAAGGTTTGAAAGCGGCTCCTGCCCGGCGCTAAGTGTATTTACACACAGCGTTCGATTTTCATTGAATGCCGGCCAGACGGACGCCCCACGATTCAGGCACACCAGTAATGTggg                             >  1:84686/1‑251 (MQ=255)                                                                                                                                                                                                    tGGGCGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCGG‑AGCCTCTTCCAGTTGCGGGCATCCGGTCACCCCCGTCTGCCAGCGGGCGGCGGCCAAGCGGTGTTCCATGGGCGTTTTGCCGCCAAAAAGGTTTGAAAGCGGCTCCTGCCCGGCGCTAAGTGTATTTACACACCGCGTTCGATTTTCATTGACTGCCGGCCCGCCGGACGCCCCACGATTCAGGCACACCAGTAATGTggg                             >  2:261761/1‑251 (MQ=255)                                                                                                                                                                                                          ccGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCGG‑AGCCTCTTCCAGTTGCGGGCATCCGGTCACCCCCGTCTGCCAGCGGGCGGCGGCAAAGCGGTGTTCCATGGGCGTTTTGCCGCCAAAAAGGTTTGAAAGCGGCTCCTGCCCGGCGCTAAGTGTATTTACACACAGCGTTCGATTTTCATTGAATGCCGGCCAGACGGACGCCCCACGATTCAGGCACACCAGTAATGTGGGCGGCGt                       >  1:205668/1‑251 (MQ=255)                                                                                                                                                                                                            gACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCGG‑AGCCTCTTCCAGTTGCGGGCATCCGGTCACCCCCGTCTGCCAGCGGGCGGCGGCAAAGCGGTGTTCCATGGGCGTTTTGCCGCCAAAAAGGTTTGAAAGCGGCTCCTGCCCGGCGCTAAGTGTATTTACACACAGCGTTAGATTTTCATTGAATGCCGGCCAGACGGACGCCCCACGATTCAGGCACACCAGTAATGTGGGCGGCGTAt                     >  2:15501/1‑251 (MQ=255)                                                                                                                                                                                                                          ggCTGATACGGCAGTCAAACGAAACCGG‑AGCCTCTTCCAGTTGCGGGCATCCGGTCACCCCCGTCTGCCAGCGGGCGGCGGCAAAGCGGTGTTCCATGGGCGTTTTGCCGCCAAAAAGGTTTGAAAGCGGCTCCTGCCCGGCGCTAAGTGTATTTACACACAGCGTTCGAGTTTCATTGAATGCCGGCCAGACGGACGCCCCACGATTCAGGCACACCAGTAATGTGGGGGGCGTATCGGTCACACTGCAg       <  1:353138/251‑1 (MQ=255)                                                                                                                                                                                                                          ggCTGATACGGCAGTCAAACGAAACCGG‑AGCCTCTTCCAGTTGCGGGCAGCCGGTCACCCCCGTCTGCCAGCGGGCGGCGGCAAAGCGGTGTTCCATGGGCGTTTTGCCGCCAAAAAGGTTTGAAAGCGGCTCCTGCCCGGCGCTAAGTGTATTTACACACAGCGTTCGATTTTCATTGAATGCCGGCCAGACGGACGCCCCACGATTCAGGCACACCAGTAATGTGGGCGGCGTATCGGTCACACTGc         <  2:65551/249‑1 (MQ=255)                                                                                                                                                                                                                               atatggCAGTCAAAAGAAACCGG‑AGCCTCTTCCAGTTGCGGGCATCCGGTCACCCCCGTCTGCCAGCGGGCGGCGGCAAAGCGGTGTTCCATGGGCGTTTTGCCGACAAAAAGGTTTGAAAGCGGCTCCTGCCCGGCGCTAAGTGTATTTACACACAGCGTTCGATTTTCATTGAATGCCGGCCAGACGGACGCCCCACGATTGAGGCACACCAGTAATGTGGGCGGCGTCTCGGTCACACTGCAGTCggg  >  2:196230/5‑250 (MQ=255)                                                                                                                                                                                                                               aTACGGCAGTCAACCGAAACCGG‑AGCCGCTTCCAGTTGCGGGCATCCGGGCACCCCCTTCTTCCAGCGGGCGGCGGCAAAGCGGTGTTCCATGGGCGTTGTGCCGCCAAAAAGGTTTGAAAGCGGCTCCTGCCCGGCGCTAAGTGTATTTACACACAGCGTTCGATTTTCATTGAATGCCGGCCAGACGGACGCCCCACGATTCAGGCACACCAGTAATGTGGGCGGCGTATCGGTCACACTGCAGACGgc  <  2:347276/251‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                       CGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCAG‑AGCCTCTTCCAGTTGCGGGCATCCGGTCACCCCCGTCTGCCAGCGGGCGGCGGCAAAGCGGTGTTCCATGGGCGTTTTGCCGCCAAAAAGGTTTGAAAGCGGCTCCTGCCCGGCGCTAAGTGTATTTACACACAGCGTTCGATTTTCATTGAATGCCGGCCAGACGGACGCCCCACGATTCAGGCACACCAGTAATGTGGGCGGCGTATCGGTCACACTGCAGACGGC  >  NC_000913/1069771‑1070243

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