BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,577,657	A→G	9.6%	intergenic (‑65/+74)	gntK ← / ← gntR	gluconate kinase 2/d‑gluconate inducible gluconate regulon transcriptional repressor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,577,657	0	A	G	9.6%	543.4 / 34.5	174	intergenic (‑65/+74)	gntK/gntR	gluconate kinase 2/d‑gluconate inducible gluconate regulon transcriptional repressor
Reads supporting (aligned to +/- strand): ref base A (78/79); new base G (8/9); total (86/88)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

Only 100 of 175 total aligned reads displayed.

CCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGATTCGCCACGAATACGCGCCAGCAGGCGTTCAGCGCCAATACTGCCCATCCGCTCACGCGGCGTCAGCACGCTCGCAAGTCGTGGCTCCATCACCTGACCAATGTCATGACCGTGGAAACCGGCAATCGCCATATCGTCAGG  >  NC_000913/3577391‑3577931

   gCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCAACAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCAAAATACAAGAATTAAAAATGTTACGGGTAACTGTTACCGGTAAAATTGTCCAACCGGAAAata                                                                                                                                                                                                                                                                         >  4:204369/1‑275 (MQ=255)

             tttGCGATCGTCGTCATTCAGTGGTTCGCCAGACGTCATTTTTTCGCTATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAg                                                                                                                                                                                                                                                                  <  2:418272/272‑1 (MQ=255)

                 cGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATatc                                                                                                                                                                                                                                                           <  2:313769/275‑1 (MQ=255)

                     cgtcgtCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATCTCAGAAGCAATatccatc                                                                                                                                                                                                                                                       >  2:416270/1‑275 (MQ=255)

                            aTTCAGGGGGTCGCCAGACGCCATTTCTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATGAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGcc                                                                                                                                                                                                                                                 <  2:115871/274‑1 (MQ=255)

                             ttCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCtt                                                                                                                                                                                                                                               <  2:169908/275‑1 (MQ=255)

                                                              ttGCGCCGTGGATGGGGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTtgtg                                                                                                                                                                                                              <  4:272698/275‑1 (MQ=255)

                                                                                    cGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCg                                                                                                                                                                                          >  1:127292/1‑273 (MQ=255)

                                                                                         tCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCa                                                                                                                                                                                                                                                                      >  3:80191/1‑192 (MQ=255)

                                                                                             gAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAg                                                                                                                                                                                                                                                                  <  1:80191/192‑1 (MQ=255)

                                                                                                               gatcagcCACGTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACTGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACAtt                                                                                                                                                             <  4:177072/270‑1 (MQ=255)

                                                                                                                               ggCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCGGCTGGACa                                                                                                                                                                                                                                                                            >  6:281957/1‑148 (MQ=255)

                                                                                                                                           tttGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCGGCTGGACaa                                                                                                                                                                                                                                                                           >  6:189421/1‑137 (MQ=255)

                                                                                                                                                           aCGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAat                                                                                                                                                                                                                                                                          <  5:530479/122‑1 (MQ=255)

                                                                                                                                                           aCGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAat                                                                                                                                                                                                                                                                          >  6:530479/1‑122 (MQ=255)

                                                                                                                                                           aCGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGATTCGCCACGAATACGCGCCAGCAGGCGTTc                                                                                                                 <  2:102234/275‑1 (MQ=255)

                                                                                                                                                            cGCCAAGCAAGACGTAAATGGGGTGATCATGTTTAGTAGTGCTCAAATTGGTGCCCCCAAAATACAAGAATTAAAANTGTTACGGGTAAATGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGTAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGATTCGCCACGAATACGCGCCAGCAGGCGTTc                                                                                                                 <  3:127292/274‑1 (MQ=255)

                                                                                                                                                            ataaaagataGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCACCACAAAACAAAAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGATTCGCCACGAATACGCGCCAGCAGGCGTTc                                                                                                                 <  1:155589/265‑1 (MQ=255)

                                                                                                                                                             gCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTa                                                                                                                                                                                                                                              <  7:355687/148‑1 (MQ=255)

                                                                                                                                                              cccATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTaa                                                                                                                                                                                                                                             >  8:706710/1‑148 (MQ=255)

                                                                                                                                                                 aTCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCGGCTGGACaa                                                                                                                                                                                                                                                                           <  6:887619/115‑1 (MQ=255)

                                                                                                                                                                 aTCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCGGCTGGACaa                                                                                                                                                                                                                                                                           >  7:647831/1‑115 (MQ=255)

                                                                                                                                                                     aGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTa                                                                                                                                                                                                                                              <  8:746980/140‑1 (MQ=255)

                                                                                                                                                                          tAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCGGCTGGACaa                                                                                                                                                                                                                                                                           >  6:992566/1‑106 (MQ=255)

                                                                                                                                                                          tAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCGGCTGGACaa                                                                                                                                                                                                                                                                           <  5:992566/106‑1 (MQ=255)

                                                                                                                                                                           aaaTGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATc                                                                                                                                                                                                >  3:12868/1‑180 (MQ=255)

                                                                                                                                                                              tgtgGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAg                                                                                                                                                                                                                                            <  5:896750/133‑1 (MQ=255)

                                                                                                                                                                               gtggtgATCATGGTTAGTCGTGCTCAAAGTGGGGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGcc                                                                                                                                                                                            <  1:12868/180‑1 (MQ=255)

                                                                                                                                                                                 ggtgATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAat                                                                                                                                                                                                                                                                          >  7:737241/1‑100 (MQ=255)

                                                                                                                                                                                 ggtgATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTAc                                                                                                                                                                                                                          >  7:600909/1‑148 (MQ=255)

                                                                                                                                                                                  gtgATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACa                                                                                                                                                                                                                                                                            >  8:693102/1‑97 (MQ=255)

                                                                                                                                                                                   tgATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGtgtg                                                                                                                                                                                                                                >  6:346639/1‑140 (MQ=255)

                                                                                                                                                                                     aTCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCa                                                                                                                                                                                                                      >  6:1026111/1‑148 (MQ=255)

                                                                                                                                                                                       cATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAata                                                                                                                                                                                                                                                                         <  6:623705/95‑1 (MQ=255)

                                                                                                                                                                                           gTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTtg                                                                                                                                                                                                                <  5:1026111/148‑1 (MQ=255)

                                                                                                                                                                                            ttAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGtata                                                                                                                                                                                                                                      >  8:297191/1‑125 (MQ=255)

                                                                                                                                                                                               gTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGATTCGCCACGAATACGCGCCAGCAGGCTTTCAGCGCCAATACTGCCCATCCGCCCACACGGCGTCa                                                                              >  4:484515/1‑274 (MQ=255)

                                                                                                                                                                                                tCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGtgtg                                                                                                                                                                                                                                <  8:844458/127‑1 (MQ=255)

                                                                                                                                                                                                 tctGCTCAAAGTGTGTCCCCCAAAAGACCAGAATTACAAATTTTACGGGGAACTGTTATCGTGAAGATTAGCCAACCGGACAAGATCAGAAGAAAATTCAATGCATGCCTTAAGTGTAGAAGTGTGAGAGACTTCAAATGTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGAGTCACAGATTCGCCACGAATACGCGCCAGCAGGCGGTCAGCGCCAATACTGCCCATCCGCTCACGCGGCGTCAGc                                                                            <  4:390497/272‑1 (MQ=255)

                                                                                                                                                                                                    gCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTaa                                                                                                                                                                                                       <  7:700747/148‑1 (MQ=255)

                                                                                                                                                                                                    gCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGATTCGCCACGAATACGCGCCAGCAGGCGTTCAGCGCCAATACTGCCCATCCGCTCACGCGGCCTCAGCACg                                                                         >  2:478569/1‑274 (MQ=255)

                                                                                                                                                                                                         aaGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGa                                                                                                                                                                                                  >  7:974468/1‑148 (MQ=255)

                                                                                                                                                                                                          aGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAgtg                                                                                                                                                                                                                                  <  7:1018731/115‑1 (MQ=255)

                                                                                                                                                                                                          aGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAgtg                                                                                                                                                                                                                                  >  8:1018731/1‑115 (MQ=255)

                                                                                                                                                                                                                gCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTCCAAAtt                                                                                                                                                                                                                  <  7:127085/125‑1 (MQ=255)

                                                                                                                                                                                                                gCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTCCAAAtt                                                                                                                                                                                                                  >  8:127085/1‑125 (MQ=255)

                                                                                                                                                                                                                 cccccACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTa                                                                                                                                                                                                                                              <  8:576916/96‑1 (MQ=255)

                                                                                                                                                                                                                     cacaATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATg                                                                                                                                                                                                                                                   <  5:354034/87‑1 (MQ=255)

                                                                                                                                                                                                                     cacaATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATg                                                                                                                                                                                                                                                   >  6:354034/1‑87 (MQ=255)

                                                                                                                                                                                                                      acaATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGATTCGCCACGAATACGCCCCAGCAGGCTTTCAGCGCCAATCCTGCCCATCCGCTCACGCGGCCTCAGCCCGCTCCCAAGTCGTTGcccc                                                       >  4:471969/1‑274 (MQ=255)

                                                                                                                                                                                                                         aTACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTtg                                                                                                                                                                                                                <  7:30054/118‑1 (MQ=255)

                                                                                                                                                                                                                              aGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCTTCCATGCTTTAAGTGTATAAGtgtg                                                                                                                                                                                                                                >  8:612842/1‑97 (MQ=255)

                                                                                                                                                                                                                                  ttAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAgg                                                                                                                                                                                 <  5:815554/140‑1 (MQ=255)

                                                                                                                                                                                                                                  ttAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGATTCGCCACGAATACGCGCCAGCAGGCGTTCAGCGCCAATACTGCCCATCCGCCCACGCGGCTTCACCCCCCTCGCAAGTCGTGGCCCCCTCACCTGACCaa                                          >  4:422172/1‑275 (MQ=255)

                                                                                                                                                                                                                                  ttAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGATTCGCCACGAATACGCGCCAGCAGGAGTTCCGCGCCAATACTGCCCATCCGCTCACGCGGCGTCAGCCCGCTAGCACGCCGTGGCTCCATACCATGACCa                                           >  4:240860/1‑274 (MQ=255)

                                                                                                                                                                                                                                    aaaaaTGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGATTCGCCACGAATACGCGCCAGCAGGCGTTCAGCGCCAATACTGCCCATCCGCTCACGCGGCGTCAGCACGCCCGCAAGTCGTGGCTCCATCACCTGACCAATg                                        >  2:387251/1‑275 (MQ=255)

                                                                                                                                                                                                                                      aaaTGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATatccatc                                                                                                                                                                                                                                                       <  8:1096705/66‑1 (MQ=255)

                                                                                                                                                                                                                                      aaaTGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCaa                                                                                                                                                                                                                     <  5:73033/100‑1 (MQ=255)

                                                                                                                                                                                                                                      aaaTGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCaa                                                                                                                                                                                                                     >  6:73033/1‑100 (MQ=255)

                                                                                                                                                                                                                                       aaTGTTACGGGTAACTGTTATCGGTAACATTGTCCGGCTGGACAatg                                                                                                                                                                                                                                                                         <  8:918227/47‑2 (MQ=37)

                                                                                                                                                                                                                                        aTGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGcc                                                                                                                                                                                            <  6:685398/123‑1 (MQ=255)

                                                                                                                                                                                                                                        aTGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGcc                                                                                                                                                                                            >  5:685398/1‑123 (MQ=255)

                                                                                                                                                                                                                                           ttACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTaa                                                                                                                                                                                                       >  7:791045/1‑109 (MQ=255)

                                                                                                                                                                                                                                           ttACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTaa                                                                                                                                                                                                       <  8:791045/109‑1 (MQ=255)

                                                                                                                                                                                                                                                ggTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAg                                                                                                                                                                                                                                            >  8:155964/1‑67 (MQ=255)

                                                                                                                                                                                                                                                ggTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAg                                                                                                                                                                                                                                            <  7:155964/67‑1 (MQ=255)

                                                                                                                                                                                                                                                ggTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCgg                                                                                                                                                                                         <  7:32596/118‑1 (MQ=255)

                                                                                                                                                                                                                                                ggTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCgg                                                                                                                                                                                         >  8:32596/1‑118 (MQ=255)

                                                                                                                                                                                                                                                ggTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACa                                                                                                                                                               <  8:159825/144‑1 (MQ=255)

                                                                                                                                                                                                                                                ggTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACa                                                                                                                                                               >  7:159825/1‑144 (MQ=255)

                                                                                                                                                                                                                                                  tAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAAt                                                                                                                                                                                                     <  8:631597/104‑1 (MQ=255)

                                                                                                                                                                                                                                                  tAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAAt                                                                                                                                                                                                     >  7:631597/1‑104 (MQ=255)

                                                                                                                                                                                                                                                       gTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTtgtg                                                                                                                                                                                                              >  5:563731/1‑90 (MQ=255)

                                                                                                                                                                                                                                                         tATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATg                                                                                                                                                                                                                                                   >  4:558499/1‑51 (MQ=255)

                                                                                                                                                                                                                                                         tATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCaa                                                                                                                                                                                                                     >  6:57711/1‑81 (MQ=255)

                                                                                                                                                                                                                                                         tATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTcac                                                                                                                                                  >  6:110593/1‑148 (MQ=255)

                                                                                                                                                                                                                                                         tATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGATTCGCCACGAATAcgcg                                                                                                                              >  4:329167/1‑168 (MQ=255)

                                                                                                                                                                                                                                                          aTCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGATTCGCCACGAATACGCGCCAGCAGGCGTTCAGCGCCAATACTGCCCATCCGCTCACGCGGCGTCAGCACGCTCGCAAGTCGTGGCGCCATCACCTGACCAATGTCATGACCGTTGAAACCGGCaa                  >  2:266856/1‑275 (MQ=255)

                                                                                                                                                                                                                                                            cGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAgtgt                                                                                                                                                                                                                                 <  7:828107/66‑1 (MQ=255)

                                                                                                                                                                                                                                                             ggTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCtt                                                                                                                                                                                                                                               <  2:558500/51‑1 (MQ=255)

                                                                                                                                                                                                                                                              gTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTc                                                                                                                                                          >  5:284917/1‑135 (MQ=255)

                                                                                                                                                                                                                                                                aaCATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATcc                                                                                                                                                                                               >  5:187737/1‑96 (MQ=255)

                                                                                                                                                                                                                                                                aaCATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACa                                                                                                                                                                                    >  6:915678/1‑107 (MQ=255)

                                                                                                                                                                                                                                                                aaCATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACa                                                                                                                                                                                    <  5:13896/107‑1 (MQ=255)

                                                                                                                                                                                                                                                                  cATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCaaa                                                                                                                                                                                                                    >  3:26217/1‑73 (MQ=255)

                                                                                                                                                                                                                                                                   aTTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTaa                                                                                                                                                                                                       >  6:925345/1‑85 (MQ=255)

                                                                                                                                                                                                                                                                   aTTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTaa                                                                                                                                                                                                       <  5:925345/85‑1 (MQ=255)

                                                                                                                                                                                                                                                                    ttGTCCGGCTGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACCTTTTCGGTGTCACAGAtt                                                                                                                                             <  6:233684/142‑1 (MQ=255)

                                                                                                                                                                                                                                                                    ttGTCCGGCTGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGAtt                                                                                                                                             >  5:233684/1‑142 (MQ=255)

                                                                                                                                                                                                                                                                    ttGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGa                                                                                                                                                                                                  <  5:238439/89‑1 (MQ=255)

                                                                                                                                                                                                                                                                    ttGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTa                                                                                                                                                                  <  5:195619/121‑1 (MQ=255)

                                                                                                                                                                                                                                                                    ttGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTa                                                                                                                                                                  >  6:195619/1‑121 (MQ=255)

                                                                                                                                                                                                                                                                     tGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGATTCGCCACg                                                                                                                                      <  5:581416/148‑1 (MQ=255)

                                                                                                                                                                                                                                                                      gTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTtg                                                                                                                                                                                                                <  1:26217/73‑1 (MQ=255)

                                                                                                                                                                                                                                                                      gTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGa                                                                                                                                                                                                  >  7:889219/1‑87 (MQ=255)

                                                                                                                                                                                                                                                                       tCCGGCTGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGa                                                                                                                                                                                                  <  6:150612/86‑1 (MQ=255)

                                                                                                                                                                                                                                                                       tCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGATTCGCCACGAATACGCGCCAGCAGGCGTTCAGCGCCAATACTGCCCATCCGCTCACGCGGCGTCAGCACGCTCGCAAGTCGTTGCTCCATCACCTGACCAATGTCCTGCCCGTGGAAACCGGCAATCCCCATATCgg      >  4:4889/1‑273 (MQ=255)

                                                                                                                                                                                                                                                                        ccAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACa                                                                                                                                                                                    >  7:1139778/1‑99 (MQ=255)

                                                                                                                                                                                                                                                                         cAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGATTCGCCACGAATa                                                                                                                                  <  5:110593/148‑1 (MQ=255)

                                                                                                                                                                                                                                                                          aGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGATTCGCCACGAATACGCGCCAGCAGGCGTTCAGCGCCAATACTGCCCATCCGCTCACGCGGCGTCAGCACGCTCGCAAGTCGTGGCTCCATCACCTGACCAATGTCATGACCGTGGAAACCGGCAATCGCCATATCGTCAgg  <  2:205703/275‑1 (MQ=255)

CCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAGTTGATGCGCCACTTCACTGGCGACCGCAGATTTGCCGCTGCCCGATACGCCCATCAAGACGTAAATGTGGTGATCATGGTTAGTCGTGCTCAAAGTGGTGCCCCCACAATACAAGAATTAAAAATGTTACGGGTAACTGTTATCGGTAACATTGTCCAGCCGGACAATATCAGAAGCAATATCCATCCATGCCTTAAGTGTATAAGTGTGAGCTACTTCAAATTTGTGGGCTTAAATAGATCCGCCCGGTGACAAGGTGAAACCTAAATCTAACATTTTCGGTGTCACAGATTCGCCACGAATACGCGCCAGCAGGCGTTCAGCGCCAATACTGCCCATCCGCTCACGCGGCGTCAGCACGCTCGCAAGTCGTGGCTCCATCACCTGACCAATGTCATGACCGTGGAAACCGGCAATCGCCATATCGTCAGG  >  NC_000913/3577391‑3577931

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