BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	69,014	A→G	F345F (TTT→TTC)	araB ←	ribulokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	69,014	0	A	G	97.5%	131.2 / NA	40	F345F (TTT→TTC)	araB	ribulokinase
Reads supporting (aligned to +/- strand): ref base A (0/0); major base G (21/18); minor base T (0/1); total (21/19)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

GAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCAGGCGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCAGTTGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCCGAGTACGCGACCAAACCAGGCGTAGATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAACCTGACCGCAAATACCTTTAACTGCCCGCTCGCCAACGCTCTGTTTGTCGGCAATCAGAATGTCGCAGGTGGAAGTACCGATAACTTTTACCAGTGCGTTAGGCTGTGCGCCTGCGCCAACTGCGCCCATATGGCAGTCAAACGCGCCGCCGGAAATCACCACGCTTTCA  >  NC_000913/68768‑69260

gAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCAGTTGTTTCTGGCTGGCGTTGATTTGCGCTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCTGAGTACGCGACCGAAcc                                                                                                                                                                                                                                                    <  2:41746/251‑1 (MQ=255)

    agcagcGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCAGTTGTTTCTGGCTGGCGTTGATTTGCGCTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCTGAGTACGCGACCGAACCAGGc                                                                                                                                                                                                                                                <  2:234313/251‑1 (MQ=255)

         cGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCAGTTGTTTCTGGCTGGCGTTGATTTGCGCTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCTGAGTACGCGACCGAACCAGGCGTAGa                                                                                                                                                                                                                                           <  1:60662/251‑1 (MQ=255)

           gAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCAGTTGTTTCTGGCTGGCGTTGATTTGCGCTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCTGAGTACGCGACCGAACCAGGCGTAGat                                                                                                                                                                                                                                          >  1:197638/1‑250 (MQ=255)

                     aGCGAGGTTAAGATCGGTAATCACCCCTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCAGTTGTTTCTGGCTGGCGTTGATTTGCGCTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCTGAGTACGCGACCGAACCAGGCGTAGATATCACCAAACg                                                                                                                                                                                                                               <  1:134874/251‑1 (MQ=255)

                        gAGGTTAAGATCGGTAATCACCCCTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCAGTTGTTTCTGGCTGGCGTTGATTTGCGCTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCTGAGTACGCGACCGAACCAGGCGTAGATATCACCAAACGCCg                                                                                                                                                                                                                            >  1:65206/1‑251 (MQ=255)

                          ggTTAAGATCGGTAATCACCCCTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCAGTTGTTTCGGGCTGGCGTTGATTTGCGCTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCTGAGTACGCGACCGAACCAGGCGTAGATATCACCAAACGCCGAt                                                                                                                                                                                                                          <  1:203892/251‑1 (MQ=255)

                               aGATCGGTAATCACCCCTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCAGTTGTTTCTGGCTGGCGTTGATTTGCGCTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCTGAGTACGCGACCGAACCAGGCGTAGATATCACCAAACGCCGATTGGcc                                                                                                                                                                                                                     <  1:151060/251‑1 (MQ=255)

                                gATCGGTAATCACCCCTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCAGTTGTTTCTGGCTGGCGTTGATTTGCGCTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCTGAGTACGCGACCGAACCAGGCGTAGATATCACCAAACGCCGATTGGCCt                                                                                                                                                                                                                    <  2:197638/251‑1 (MQ=255)

                                       aaTCACCCCTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCAGTTGTTTCTGGCTGGCGTTGATTTGCGCTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCTGAGTACGCGACCGAACCAGGCGTAGATATCACCAAACGCCGATTGGCCTGCTTcc                                                                                                                                                                                                              <  1:264768/250‑1 (MQ=255)

                                        aTCACCCCTTTCAGGCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCAGTTGTTTCTGGCTGGCGTTGATTTGCGCTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCTGAGTACGCGACCGAACCAGGCGTAGATATCACCAAACGCCGATTGGCCTGCTTcc                                                                                                                                                                                                              >  1:9127/1‑249 (MQ=255)

                                                     ggCGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCAGTTGTTTCTGGCTGGCGTTGATTTGCGCTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCTGAGTACGCGACCGAACCAGGCGTAGATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATcc                                                                                                                                                                                                >  1:70768/1‑250 (MQ=255)

                                                       cGTTGGTTAGCGTTTGGCGAGCGACGACCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGGGACGGATTTTTGGCCCATGCCTCGGTCAGCGCCGGAAGCAGTTGTTTCTGGCTGGCGTTGATTTGCGCTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCTGAGTACGCGACCGAACCAGGCGTAGATCTCACAAAACGCCGATTGGCCTGCTTCCCGAACGATAAAGCCAgg                                                                                                                                                                                             >  2:202000/1‑251 (MQ=255)

                                                              tagcgtttggcgagcgacgaccGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCAGTTGTTTCTGGCTGGCGTTGATTTGCGCTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCTGAGTACGCGACCGAACCAGGCGTAGATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGcaccac                                                                                                                                                                                       <  1:188412/230‑1 (MQ=255)

                                                               agcgtttggcgagcgacgaccGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCAGTTGTTTCTGGCTGGCGTTGATTTGCGCTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCTGAGTACGCGACCGAACCAGGCGTAGATATAACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCAAAGcacccccc                                                                                                                                                                                     >  2:10732/20‑247 (MQ=255)

                                                                agcgttggcgagcgacgaccGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCAGTTGTTTCTGGCTGGCGTTGATTTGCGCTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCTGAGTACGCGACCGAACCAGGCGTAGATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGCt                                                                                                                                                                                    >  2:154511/19‑251 (MQ=255)

                                                                         gagcgacgaccGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCAGTTGTTTCTGGCTGGCGTTGATTTGCGCTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCTGAGTACGCGACCGAACCAGGCGTAGATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAAc                                                                                                                                                                           <  1:154511/242‑1 (MQ=255)

                                                                           gcgacgaccGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCAGTTGTTTCTGGCTGGCGTTGATTTGCGCTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCTGAGTACGCGACCGAACCAGGCGTAGATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAACCt                                                                                                                                                                         >  2:237226/8‑251 (MQ=255)

                                                                             gacgaccGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCAGTTGTTTCTGGCTGGCGTTGATTTGCGCTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCTGAGTACGCGACCGAACCAGGCGTAGATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAACCTGa                                                                                                                                                                       <  2:173460/246‑1 (MQ=255)

                                                                                        aaCCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCAGTTGTTTCTGGCTGGCGTTGATTTGCGCTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCTGAGTACGCGACCGAACCAGGCGTAGATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAACCTGACCGCAAATAcc                                                                                                                                                            <  1:32597/251‑1 (MQ=255)

                                                                                            aGTCGAGAACCACCGGCCGGGGAAACAGAGACGGATTTTTCGCCCATGCTTCGGTCAGCGCCGGAAGCAGTTGTTTCTGGCTGGCGTTGATTTGCGCTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCTGAGTAAGCGACCGAACCAGGCGGAGATATCACCAAACGCCGATTGGCCTGCTTCCCGCCCGATAAATCCAGGCACCACGCTGCCATCAACCTTACCGCAAATAACttt                                                                                                                                                         >  2:153864/1‑250 (MQ=255)

                                                                                                                       gagaCGGATTTTTGGCCCATGCTTCGGGCAGCGCCGGAAGCAGTTGTTTCTGGCTGGCGTTGATTTGCGCTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCTGAGTACGCGACCGAACCAGGCGTAGATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAACCTGACCGCAAATACCTTTAACTGCCCGCTCGCCAACGCTCTGTTTg                                                                                                                             <  1:167865/251‑1 (MQ=255)

                                                                                                                             gATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCAGTTGTTTCTGGCTGGCGTTGATTTGCGCTTTCAGTTCCGGCTGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCTGAGTACGCGACCGAACCAGGCGTAGATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAACCTGACCGCAAATACCTTTAACTGCCCGCTCGCCAACGCCCTGTTTGTCGGc                                                                                                                        >  2:137151/1‑250 (MQ=255)

                                                                                                                                    ggCCCATGCTTCGGTCAGCGCCGGAAGCAGTTGTTTCTGGCTGGCGTTGATTTGCGCTTTCAGTGCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCTGAGTACGCGACCGAACCAGGCGTAGATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAACCTGACCGCAAATACCTTTAACTGCCCGCTCGCCAACGCTCTGTTTGTCGGCAATCAGa                                                                                                                 <  2:41075/250‑1 (MQ=255)

                                                                                                                                       ccATGCTTCGGTCAGCGCCGGAAGCAGTTGTTTCTGGCTGGCGTTGATTTGCGCTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCTGAGTACGCGACCGAACCAGGCGTAGATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAACCTGACCGCAAATACCTTTAACTGCCCGCTCGCCAACGCTCTGTTTGTCGGCAATCAGAATGt                                                                                                             >  1:236589/1‑251 (MQ=255)

                                                                                                                                           gCTTCGGTCAGCGCCGGAAGCAGTTGTTTCTGGCTGGCGTTGATTTGCGCTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCTGAGTACGCGACCGAACCAGGCGTAGATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAACCTGACCGCAAATACCTTTAACTGCCCGCTCGCCAACGCTCTGTTTGTCGGCAATCAGAATGTCGc                                                                                                          >  1:95234/1‑250 (MQ=255)

                                                                                                                                              tCGGTCAGCGCCGGAAGCAGGTGTTTCTGGCTGGCGTTGATTTGCGCTTTCAGTTCCGGATGCTGGACGGCAAGCTGTTCCAGCGGCCAGCTGAGTACGCGACCGAACCAGGCGTAGATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAACCTGACCGCAAATACCGTTAACTGCCCGCTCGCCAACGCTCTGTTTGTCGGCAATCAGAATGTCGCAgg                                                                                                       <  2:71609/250‑1 (MQ=255)

                                                                                                                                                                 gttgttTCTGGCTGGCGTTGATTTGCGCTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCTGAGTACGCGACCGAACCAGGCGTAGATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAACCTGACCGCAAATACCTTTAACTGCCCGCTCGCCAACGCTCTGTTTGTCGGCAATCAGAATGTCGCAGGTGGAAGTACCGATAACtttt                                                                                   >  2:47677/1‑251 (MQ=255)

                                                                                                                                                                                    gATTTGCGCTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCTGAGTACGCGACCGAACCAGGCGTAGATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAACCTGACCGCAAATACCTTTAACTGCCCGCTCGCCAACGCTCTGTTTGTCGGCAATCAGAATGTCGCAGGTGGAAGTACCGATAACTTTTACCAGTGCGTTAGGCtgtg                                                                <  1:130607/251‑1 (MQ=255)

                                                                                                                                                                                        tGCGCTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCTGAGTACGCGACCGAACCAGGCGTAGATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAACCTGACCGCAAATACCTTTAACTGCCCGCTCGCCAACGCTCTGTTTGTCGGCAATCAGAATGTCGCAGGTGGAAGTACCGATAACTTTTACCAGTGCGTTAGGCTGtgcgcc                                                            >  2:103391/1‑251 (MQ=255)

                                                                                                                                                                                          cgctttCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCTGAGTACGCGACCGAACCAGGCGTAGATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAACCTGACCGCAAATACCTTTAACTGCCCGCTCGCCAACGCTCTGTTTGTCGGCAATCAGAATGTCGCAGGTGGAAGTACCGATAACTTTTACCAGTGCGTTAGGCTGtgcgcctg                                                          >  1:90341/4‑251 (MQ=255)

                                                                                                                                                                                           gctttCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCTGAGTACGCGACCGAACCAGGCGTAGATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAACCTGACCGCAAATACCTTTAACTGCCCGCTCGCCAACGCTCTGTTTGTCGGCAATCAGAATGTCGCAGGTGGAAGTACCGATAACTTTTACCAGTGCGTTAGGCTGtgcgcctgc                                                         >  1:18064/3‑251 (MQ=255)

                                                                                                                                                                                           gctttCAGTTCCGGATGCTGGGCGGCAAGATGTTCCAGCGGCCAGCTGAGTACGCGACCGAACCAGGCGTAGATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAACCTGACCGCAAATACCTTTAACTGCCCGCTCGCCAACGCTCTGTTTGTCGGCAATCAGAATGTCGCAGGTGGATGTACCGATAACTTTTACCAGTGCGTTAGGCTGtgcgcctgc                                                         >  1:74073/3‑251 (MQ=255)

                                                                                                                                                                                                cAGTTCCGGATGCTGGGCGGCAAGCTGTGCCAGCGGCCAGCTGAGTACGCGACCTAACCAGGCGTAGATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAGCCTGACCGCAAATACCTTTAACTGCCCGCTCGCCAACGCTCTGTTTGTCGGCAATCAGAATGTCGCAGGTGTAAGTACCGATAACTTTTACCAGTGCGTTAGGCTGTGCGCCTGCGCCa                                                     <  2:248795/250‑1 (MQ=255)

                                                                                                                                                                                                                  ggcAAGCTGTTCCAGCGGCCAGCTGAGTACGCGACCGAACCAGGCGTAGATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAACCTGACCGCAAATGCCTTTAACTGCCCGCTCGCCAACGCTCTGTTTGTCGGCAATCAGAATGTCGCAGGTGGAAGTACCGATAACTTTTACCAGTGCTTTAGGCTGTGCGCCTGCGCCAACTGCGCCCATATGGCAGt                                  <  2:186229/251‑1 (MQ=255)

                                                                                                                                                                                                                     aaGCTGTTCCAGCGGCCAGCTGAGTACGCGACCGAACCAGGCGTAGATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAACCTGACCGCAAATACCTTTAACTGCCCGCTCGCCAACGCTCTGTTTGTCGGCAATCAGAATGTCGCAGGTGGAAGTACCGATAACTTTTACCAGTGCGTTAGGCTGTGCGCCTGCGCCAACTGCGCCCATATGGCAGTCa                                >  2:39059/1‑250 (MQ=255)

                                                                                                                                                                                                                                   gCCAGCTGAGTACGCGACCGAACCAGGCGTAGATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAACCTGACCGCAAATACCTTTAACTGCCCGCTCGACAACGCTCTGTTTGTCGGCAATCAGAATGTCGCAGGTGGAAGTACCGATAACTTTTACCAGTGCGTTAGGCTGTGCGCCTGCGCCAACTGCGCCCATATGGCAGTCAAACGCGCCGCCGGaa                 <  2:62536/251‑1 (MQ=255)

                                                                                                                                                                                                                                      agctgAGTACGCGACCGAACAAGGCGTAGATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAACCTGACCGCAAATACCTTTAACTGCCCGCTCGCCAACGCTCTGTTTGGCGGCAATCAGAATGTCGCAGGCGGAAGTACCGATAACTTTTACCAGTGCGTTAGGCTGTGCGCCTGCGCCAACTGCGCCCATATGGCAGTCAAACGCGCCGCCGGAAATc              >  2:41145/5‑251 (MQ=255)

                                                                                                                                                                                                                                        ctgAGTACGCGACCGAACCAGGCGTAGATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAACCTGACCGCAAATACCTTTAACTGCCCGCTCGCCAACGCTCTGTTTGTCGGCAATCAGAATGTCGCAGGTGGAAGTACCGATAACTTTTACCAGTGCGTTAGGCTGTGCGCCTGCGCCAACTGCGCCCATATGGCAGTCAAACGCGCCGCCGGAAATcac            >  2:219271/3‑251 (MQ=255)

                                                                                                                                                                                                                                                  gACCGAACCAGGCGTAGATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAACCTGACCGCAAATACCTTTAACTGCCCGCTCGCCAACGCTCTGTTTGTCGGCAATCAGAATGTCGCAGGTGGAAGTACCGATAACTTTTACCAGTGCGTTAGGCTGTGCGCCTGCGCCAACTGCGCCCATATGGCAGTCAAACGCGCCGCCGGAAATCACCCCGCTTTcc  >  2:138188/1‑250 (MQ=255)

GAACAGCAGCGGAGCGTCGGTAGCGAGGTTAAGATCGGTAATCACCCCTTTCAGGCGTTGGTTAGCGTTCGGTGTGCGGCGGCCGTTAAACCAGTCGAGCACCACCGGCAGGTGATCCAGAGACGGATTTTTGGCCCATGCTTCGGTCAGCGCCGGAAGCAGTTGTTTCTGGCTGGCGTTGATTTGCGTTTTCAGTTCCGGATGCTGGGCGGCAAGCTGTTCCAGCGGCCAGCCGAGTACGCGACCAAACCAGGCGTAGATATCACCAAACGCCGATTGGCCTGCTTCCAGACCGATAAATCCAGGCACCACGCTGCCATCAACCTGACCGCAAATACCTTTAACTGCCCGCTCGCCAACGCTCTGTTTGTCGGCAATCAGAATGTCGCAGGTGGAAGTACCGATAACTTTTACCAGTGCGTTAGGCTGTGCGCCTGCGCCAACTGCGCCCATATGGCAGTCAAACGCGCCGCCGGAAATCACCACGCTTTCA  >  NC_000913/68768‑69260

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 13 ≤` `ATCG/ATCG` `< 24 ≤` `ATCG/ATCG` `< 35 ≤` `ATCG/ATCG` `< 39 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`