BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,315,800	A→G	6.8%	intergenic (‑30/+56)	yciF ← / ← yciG	putative rubrerythrin/ferritin‑like metal‑binding protein/KGG family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,315,800	0	A	G	6.8%	571.4 / 20.0	177	intergenic (‑30/+56)	yciF/yciG	putative rubrerythrin/ferritin‑like metal‑binding protein/KGG family protein
Reads supporting (aligned to +/- strand): ref base A (76/88); new base G (3/9); total (79/98)
Fisher's exact test for biased strand distribution p-value = 2.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.95e-01

Only 100 of 179 total aligned reads displayed.

TCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCGTTGCG  >  NC_000913/1315535‑1316053

tcaaTTGCCACACATTTCATGCGCTTAATTTTCAGGTTAGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTTCCTGAGTTTCATCGAGATGCGCATGAAAAGCCTGACTTAATTGTTCATTTGATGTTGCTCTGCCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATTTTCAATGGTCTTCATTTTCATATTTTTCTCAATTGAAATCCATTCGATCCACACGGcc                                                                                                                                                                                                                                                        >  4:333249/4‑273 (MQ=255)

         aaacctttCTTGCGCTTAATTTTCCGGTTCGATTCGGATGCCTCAACTTGATCAATACCTTCGTTATGTCCATGAGTTTCCGCGAGTTGCGGATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTTTTGCGAGTTGTGCCAGTGCCCGGGTTAATTGTTGTTCTGCGCTGTAGGTATCTGAAAGCATGTGAATAAATACATCATCGATGGTCTTCATATTCATATTTTTCTCCAGGGAAATCACTGCGAGCCAAACGGCCCGCAGATGCag                                                                                                                                                                                                                                             <  4:281176/270‑2 (MQ=255)

               ttCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCTCGCGGTGCGCATGACAAGCCTGACTTAATCTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCCGTGCCCGGGTTAATTGTTTTTCTGCGCCGTAGGTATCTGAAAGCAGGTGAATAAATACATTTTAATTGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGACATACGGCCCGCAGATGCACTTATca                                                                                                                                                                                                                                       >  2:160599/1‑274 (MQ=255)

               ttCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACTTCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCCCTGCGAGCCATACGGCCCGCAACTGAAATTATCt                                                                                                                                                                                                                                       >  4:275329/1‑275 (MQ=255)

                         tAAGTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCTCGAGGGGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTg                                                                                                                                                                                                                             <  2:85508/275‑1 (MQ=255)

                                         ttCGGACTCCACAAATGTATTAATACGTTCAAGCTGTCCATGAGTTTCCTCGAGTTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACg                                                                                                                                                                                                             <  4:173785/275‑1 (MQ=255)

                                                            tCAATACGTTCAATCTGTCCATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGa                                                                                                                                                                                                                                                  >  2:113928/1‑219 (MQ=255)

                                                                              ccATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAAt                                                                                                                                                                                                       <  2:143866/244‑1 (MQ=255)

                                                                                       tCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTTAATCAGGCTTTGCCTGATTTATTACCACCGTTTTGTTGAACGc                                                                                                                                                               >  4:104634/1‑275 (MQ=255)

                                                                                                                     aaTTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCa                                                                                                                                                                                                                                                               <  6:122322/149‑1 (MQ=255)

                                                                                                                       tttttCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATa                                                                                                                                                                                                                                                             >  5:848872/1‑149 (MQ=255)

                                                                                                                       tttttCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATa                                                                                                                                                                                                                                                             <  6:182058/149‑1 (MQ=255)

                                                                                                                         tttCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGTGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACg                                                                                                                                                                                                                                                           >  5:327271/1‑149 (MQ=255)

                                                                                                                              tttGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTATTCTCCAGTGAAATCACTGCGAGCCATACGGCCCg                                                                                                                                                                                                                                                      >  8:952366/1‑149 (MQ=255)

                                                                                                                              tttGATGTTGCTCTTGCGAGTTTTGCCAGTGCACGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGGCATCATATTCATATTTTTCTCCAGTGAAATCACAGCGAGCCATACGGCCCg                                                                                                                                                                                                                                                      >  8:912898/1‑149 (MQ=255)

                                                                                                                                 gATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGGGCCGTATGGCTCGCAg                                                                                                                                                                                                                                                   <  7:208561/149‑1 (MQ=255)

                                                                                                                                 gATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAg                                                                                                                                                                                                                                                   >  8:599762/1‑149 (MQ=255)

                                                                                                                                  aTGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGGGCCGTATGGCTCGCAgt                                                                                                                                                                                                                                                  <  7:1340543/149‑2 (MQ=255)

                                                                                                                                    gTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGGGCCGTATGGCTCGCAgtga                                                                                                                                                                                                                                                <  7:1269856/149‑4 (MQ=255)

                                                                                                                                    gTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATg                                                                                                                                                                                                                                                <  8:1319145/149‑1 (MQ=255)

                                                                                                                                     ttGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGGGCCGTATGGCTCGCAgtgat                                                                                                                                                                                                                                               <  6:967256/149‑5 (MQ=255)

                                                                                                                                     ttGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGGGCCGTATGGCTCGCAgtgat                                                                                                                                                                                                                                               <  6:1453502/149‑5 (MQ=255)

                                                                                                                                           ttGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCaa                                                                                                                                                                                                                                             >  6:506301/1‑145 (MQ=255)

                                                                                                                                           ttGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTAt                                                                                                                                                                                                                                         <  6:1322305/149‑1 (MQ=255)

                                                                                                                                            tGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCACTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATc                                                                                                                                                                                                                                        <  7:952366/149‑1 (MQ=255)

                                                                                                                                               gAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCttt                                                                                                                                                                                                                                     <  7:256711/149‑1 (MQ=255)

                                                                                                                                                 gTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTAt                                                                                                                                                                                                                                   <  5:723431/149‑1 (MQ=255)

                                                                                                                                                   tttGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCa                                                                                                                                                                                                                                 >  5:948388/1‑149 (MQ=255)

                                                                                                                                                   tttGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCa                                                                                                                                                                                                                                 <  7:773313/149‑1 (MQ=255)

                                                                                                                                                       ccAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTg                                                                                                                                                                                                                             >  7:721927/1‑149 (MQ=255)

                                                                                                                                                       ccAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTg                                                                                                                                                                                                                             >  5:1475338/1‑149 (MQ=255)

                                                                                                                                                          gTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGGGCCGTATGGCTCGCAgt                                                                                                                                                                                                                                                  <  6:893456/125‑2 (MQ=255)

                                                                                                                                                          gTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGGGCCGTATGGCTCGCAgt                                                                                                                                                                                                                                                  >  5:893456/1‑124 (MQ=255)

                                                                                                                                                          gTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGata                                                                                                                                                                                                                          >  5:356893/1‑149 (MQ=255)

                                                                                                                                                          gTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGata                                                                                                                                                                                                                          <  7:688107/149‑1 (MQ=255)

                                                                                                                                                               cGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATtaat                                                                                                                                                                                                                     <  6:446809/149‑1 (MQ=255)

                                                                                                                                                               cGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATtaat                                                                                                                                                                                                                     >  5:1195339/1‑149 (MQ=255)

                                                                                                                                                                gggTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTTATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATtaata                                                                                                                                                                                                                    >  7:926873/1‑149 (MQ=255)

                                                                                                                                                                gggTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATtaata                                                                                                                                                                                                                    >  8:304712/1‑149 (MQ=255)

                                                                                                                                                                 ggTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATtaataa                                                                                                                                                                                                                   >  8:382571/1‑149 (MQ=255)

                                                                                                                                                                   ttAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAAtt                                                                                                                                                                                                                 <  7:615306/149‑1 (MQ=255)

                                                                                                                                                                      aTTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAAc                                                                                                                                                                                                              <  5:221247/149‑1 (MQ=255)

                                                                                                                                                                        tGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACgg                                                                                                                                                                                                            <  7:1402773/149‑1 (MQ=255)

                                                                                                                                                                                tGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAg                                                                                                                                                                                                    <  5:510088/149‑1 (MQ=255)

                                                                                                                                                                                   gcTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGAt                                                                                                                                                                                                 <  7:382571/149‑1 (MQ=255)

                                                                                                                                                                                    cTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCGGATGCGCGTTTCGGATCATTTTTAAAATTACCGCTCTTATGTTGACCGCCTTTACGGCCTGCGTCGGATACATtct                                                                  >  4:121212/1‑275 (MQ=255)

                                                                                                                                                                                       tAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGc                                                                                                                                                                                             >  5:1112364/1‑149 (MQ=255)

                                                                                                                                                                                         ggTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCt                                                                                                                                                                                           >  7:1424536/1‑149 (MQ=255)

                                                                                                                                                                                           tATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGa                                                                                                                                                                                         >  5:968238/1‑149 (MQ=255)

                                                                                                                                                                                            aTCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGAt                                                                                                                                                                                        <  5:356697/149‑1 (MQ=255)

                                                                                                                                                                                             tCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGAtt                                                                                                                                                                                       <  7:1077421/149‑1 (MQ=255)

                                                                                                                                                                                              cTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGAttt                                                                                                                                                                                      >  7:1253888/1‑149 (MQ=255)

                                                                                                                                                                                                 aaaGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTa                                                                                                                                                                                                                                    >  1:31065/1‑100 (MQ=255)

                                                                                                                                                                                                  aaGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATttatta                                                                                                                                                                                  >  5:1236849/1‑149 (MQ=255)

                                                                                                                                                                                                  aaGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATttatta                                                                                                                                                                                  <  8:119271/149‑1 (MQ=255)

                                                                                                                                                                                                   aGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTac                                                                                                                                                                                 <  8:347523/149‑1 (MQ=255)

                                                                                                                                                                                                    gCAGGTGAATAAATACATCTTCAATGGCCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTacc                                                                                                                                                                                >  6:374085/1‑149 (MQ=255)

                                                                                                                                                                                                     cAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAg                                                                                                                                                                                                                                <  3:31065/100‑1 (MQ=255)

                                                                                                                                                                                                     cAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTacca                                                                                                                                                                               <  6:968238/149‑1 (MQ=255)

                                                                                                                                                                                                            ataaataCATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCttt                                                                                                                                                                        <  8:810902/149‑1 (MQ=255)

                                                                                                                                                                                                             taaataCATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTttg                                                                                                                                                                       >  5:301035/1‑149 (MQ=255)

                                                                                                                                                                                                              aaataCATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTCAATCAGGATTTGCCTGATTTATTACCACCGCTttgt                                                                                                                                                                      >  7:1292760/1‑149 (MQ=255)

                                                                                                                                                                                                              aaataCATCTTAAAGTGTCTTCATATTCATATTTTTATCCAGTAAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGaa                                        <  2:236103/275‑1 (MQ=255)

                                                                                                                                                                                                                ataCATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTttgttg                                                                                                                                                                    >  8:502927/1‑149 (MQ=255)

                                                                                                                                                                                                                 taCATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGa                                                                                                                                                                   <  7:947703/149‑1 (MQ=255)

                                                                                                                                                                                                                    aTCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCg                                                                                                                                                                >  7:471809/1‑149 (MQ=255)

                                                                                                                                                                                                                       ttCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATttatta                                                                                                                                                                                  <  6:208835/128‑1 (MQ=255)

                                                                                                                                                                                                                       ttCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCt                                                                                                                                                             >  5:1307902/1‑149 (MQ=255)

                                                                                                                                                                                                                       ttCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCt                                                                                                                                                             <  6:656982/149‑1 (MQ=255)

                                                                                                                                                                                                                       ttCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCAGATTTATTACCACCGCTTTGTTGACCGCCt                                                                                                                                                             >  5:959487/1‑149 (MQ=255)

                                                                                                                                                                                                                         cAATGGTCTTCATATTCATATTTTTCTCCATTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGc                                                                                                                                                                                             <  1:55334/115‑1 (MQ=255)

                                                                                                                                                                                                                          aaTGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCtttt                                                                                                                                                          <  5:742409/149‑1 (MQ=255)

                                                                                                                                                                                                                           aTGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCttttt                                                                                                                                                         <  8:471809/149‑1 (MQ=255)

                                                                                                                                                                                                                           aTGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCttttt                                                                                                                                                         <  8:1442270/149‑1 (MQ=255)

                                                                                                                                                                                                                              gTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGAt                                                                                                                                                                                        <  8:578276/115‑1 (MQ=255)

                                                                                                                                                                                                                              gTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGAt                                                                                                                                                                                        >  7:578276/1‑115 (MQ=255)

                                                                                                                                                                                                                                 ttcatattcataTTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTCCCGCCGCTATGCTGCCCGCCTTTACGGCCTGCGTCGGATTCCTTCTCCCGGTCTTCGGCGGAATTTCCGGAACCACCACGATGtt                           >  2:320004/1‑269 (MQ=255)

                                                                                                                                                                                                                                   catattcatattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCt                                                                                                                                                 <  8:171071/149‑1 (MQ=255)

                                                                                                                                                                                                                                    atattcatattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCtt                                                                                                                                                >  7:1442521/1‑149 (MQ=255)

                                                                                                                                                                                                                                      attcatattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCa                                                                                                                                              >  7:828283/1‑149 (MQ=255)

                                                                                                                                                                                                                                           tattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATgc                                                                                                                                         >  6:828315/1‑149 (MQ=255)

                                                                                                                                                                                                                                           tattTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATgc                                                                                                                                         >  5:684183/1‑149 (MQ=255)

                                                                                                                                                                                                                                             tttttCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGAtt                                                                                                                                                                                                >  2:374216/1‑92 (MQ=255)

                                                                                                                                                                                                                                             tttttCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATgcgc                                                                                                                                       <  6:1307902/149‑1 (MQ=255)

                                                                                                                                                                                                                                             tttttCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATgcgc                                                                                                                                       <  6:500487/149‑1 (MQ=255)

                                                                                                                                                                                                                                                 tctcCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGcc                                                                                                                                                                                            <  4:374215/92‑1 (MQ=255)

                                                                                                                                                                                                                                                  ctcCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGc                                                                                                                                  <  8:1057848/149‑1 (MQ=255)

                                                                                                                                                                                                                                                  ctcCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGcc                                                                                            >  2:289502/1‑187 (MQ=255)

                                                                                                                                                                                                                                                    ccatttaaaaCACAGCGACCCATACGGCCCGCAAATGAAATAAACTTTATCGTCGGTAATTAATAATTAACGTTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCcgttgcg  <  4:294529/269‑1 (MQ=255)

                                                                                                                                                                                                                                                       gTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGa                                                                                                                               <  7:11204/147‑1 (MQ=255)

                                                                                                                                                                                                                                                          aaaTCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCAtt                                                                                                                          <  6:301035/149‑1 (MQ=255)

                                                                                                                                                                                                                                                          aaaTCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCAtt                                                                                                                          <  8:16124/149‑1 (MQ=255)

                                                                                                                                                                                                                                                           aaTCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCAttt                                                                                                                         >  5:184044/1‑149 (MQ=255)

                                                                                                                                                                                                                                                             tCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCtt                                                                                                                                                <  6:432561/124‑1 (MQ=255)

                                                                                                                                                                                                                                                             tCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCtt                                                                                                                                                >  5:432561/1‑124 (MQ=255)

                                                                                                                                                                                                                                                               tcTGCGGGCCGTATGGCTCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTttgtt                                                                                                                                                                     <  8:1358860/100‑1 (MQ=255)

                                                                                                                                                                                                                                                               tcTGCGGGCCGTATGGCTCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTttgtt                                                                                                                                                                     >  7:1358860/2‑101 (MQ=255)

                                                                                                                                                                                                                                                                 tGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTaaaa                                                                                                                   >  8:376595/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                   cGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAAtt                                                                                                                 >  6:1162405/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                     aGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTAc                                                                                                               <  6:334953/149‑1 (MQ=255)

TCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCGTTGCG  >  NC_000913/1315535‑1316053

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