BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	607,379	G→T	7.9%	G95C (GGC→TGC)	ribA →	GTP cyclohydrolase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	607,379	0	G	T	7.9%	312.1 / 12.4	101	G95C (GGC→TGC)	ribA	GTP cyclohydrolase II
Reads supporting (aligned to +/- strand): ref base G (48/45); new base T (2/6); total (50/51)
Fisher's exact test for biased strand distribution p-value = 2.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.63e-01

Only 100 of 104 total aligned reads displayed.

GGTACTGGGGCGCCTGCATTCCGAGTGCCTGACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGCCATATGCCTGCCGATGCTCGAGCACCTGGGCGTCAAGT  >  NC_002947/607231‑607526

ggTACTGGGGCGCCTGCATTCCGAGTGCCTGACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTg                                                                                                                                                     >  4:266961/1‑149 (MQ=255)

ggTACTGGGGCGCCTGCATTCCGAGTGCCTGACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTg                                                                                                                                                     <  5:75007/149‑1 (MQ=255)

 gTACTGGGGCGCCTGCATTCCGAGTGCCTGACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTgg                                                                                                                                                    >  3:288403/1‑149 (MQ=255)

       gggCGCCTGTATTCCGAGTGCCTGACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCgg                                                                                                                                              >  4:209689/1‑149 (MQ=255)

           gcCTGCATTCCGAGTGCCTGACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGATGCCACGGCctt                                                                                                                                          >  1:377327/1‑148 (MQ=255)

            ccTGCATTCCGAGTGCCTGACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCctgc                                                                                                                                         >  3:166617/1‑149 (MQ=255)

            ccTGCATTCCGAGTGCCTGACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTCCGTCAG‑‑GAAGGCCGATGCCACGGCcttc                                                                                                                                         <  7:44893/149‑3 (MQ=255)

             cTGCATTCCGAGTTCCTGACCGGCGATGCAATGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGATGCCACGGCcttcc                                                                                                                                        <  5:75074/149‑4 (MQ=255)

             cTGCATTCCGAGTGCCTGACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTGGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGAGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGATGCCACGGCcttcc                                                                                                                                        <  4:57518/149‑4 (MQ=255)

             cTGCATTCCGAGTGCCTGACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGATGCCACGGCcttcc                                                                                                                                        <  2:377327/149‑4 (MQ=255)

             cTGCATTCCGAGTGCCTGACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGACTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGATGCCACGGCcttcc                                                                                                                                        <  7:201003/149‑4 (MQ=255)

              tGCATTCCGAGTGCCTGACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCctgct                                                                                                                                        >  6:254044/1‑148 (MQ=255)

              tGCATTCCGAGTGCCTGACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCACCGGCctgct                                                                                                                                        >  6:360799/1‑148 (MQ=255)

                cATTCCGAGTGCCTGACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGaa                                                                                                                                     <  2:362981/149‑1 (MQ=255)

                 aTTCCGAGTGCCTGACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAAc                                                                                                                                    <  8:274279/149‑1 (MQ=255)

                  ttCCGAGTGCCTGACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAAc                                                                                                                                    >  5:188227/1‑148 (MQ=255)

                  ttCCGAGTGCCTGACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACa                                                                                                                                   <  6:188227/149‑1 (MQ=255)

                      gAGTGCCTGACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAAc                                                                                                                                    >  3:365414/1‑144 (MQ=255)

                        gTGCCTGACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCACCGGCCTGCTGAACAAGATcc                                                                                                                             <  5:360799/149‑1 (MQ=255)

                            cTGACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCgcgc                                                                                                                         >  1:12455/1‑149 (MQ=255)

                              gACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGcc                                                                                                                        >  1:203975/1‑148 (MQ=255)

                                 cGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAATGCCGATGCCACGGCcttcct                                                                                                                                       >  6:5853/1‑126 (MQ=255)

                                 cGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTAc                                                                                                                     >  3:309769/1‑148 (MQ=255)

                                 cGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACg                                                                                                                    <  3:209689/149‑1 (MQ=255)

                                 cGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGATGCCACGGCcttcct                                                                                                                                       <  5:5853/130‑5 (MQ=255)

                                  ggCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTAc                                                                                                                     >  5:73842/1‑147 (MQ=255)

                                  ggCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGa                                                                                                                   <  3:218787/149‑1 (MQ=255)

                                   gCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAAATCCGCGCCTACGAg                                                                                                                  >  1:327702/1‑149 (MQ=255)

                                         cACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCa                                                                                                             <  8:41869/148‑1 (MQ=255)

                                         cACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAg                                                                                                            >  8:234668/1‑149 (MQ=255)

                                          aCTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAg                                                                                                            >  8:348069/1‑148 (MQ=255)

                                           cTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATcc                                                                                                                             <  3:134226/130‑1 (MQ=255)

                                           cTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATcc                                                                                                                             >  4:134226/1‑130 (MQ=255)

                                              ttCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATgg                                                                                                       <  2:327702/149‑1 (MQ=255)

                                                cagccagcGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATggc                                                                                                      <  8:190717/148‑1 (MQ=255)

                                                 agccagcGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGc                                                                                                              >  7:164127/1‑139 (MQ=255)

                                                 agccagcGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGc                                                                                                              <  8:164127/139‑1 (MQ=255)

                                                 agccagcGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATggcgg                                                                                                    >  3:232317/1‑149 (MQ=255)

                                                  gccagcGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATggcggc                                                                                                   >  3:194543/1‑149 (MQ=255)

                                                    cagcGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGgcgc                                                                                                 <  8:4377/149‑1 (MQ=255)

                                                     agcGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGc                                                                                                              >  2:272109/1‑135 (MQ=255)

                                                     agcGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGcc                                                                                                >  5:68779/1‑149 (MQ=255)

                                                     agcGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAACGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGc                                                                                                              <  1:272109/135‑1 (MQ=255)

                                                      gcgcTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGcc                                                                                                >  2:102688/1‑148 (MQ=255)

                                                              aCTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTc                                                                                       <  7:348069/149‑1 (MQ=255)

                                                              aCTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTc                                                                                       <  4:365414/149‑1 (MQ=255)

                                                               cTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTc                                                                                       >  1:347373/1‑148 (MQ=255)

                                                               cTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCg                                                                                      <  7:160965/149‑1 (MQ=255)

                                                                tGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGTCGAGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCg                                                                                      >  1:324399/1‑148 (MQ=255)

                                                                        gCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGc                                                                                  <  7:295936/144‑1 (MQ=255)

                                                                        gCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGc                                                                                  >  8:295936/1‑144 (MQ=255)

                                                                        gCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAAc                                                                              >  7:342063/1‑148 (MQ=255)

                                                                         cAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAAc                                                                              >  8:40006/1‑147 (MQ=255)

                                                                         cAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGa                                                                            <  7:272867/149‑1 (MQ=255)

                                                                            ttGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAgcg                                                                         >  3:50360/1‑149 (MQ=255)

                                                                              ggAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAGATGAAGGCCGTGCCATCGCCCTGAGGAACAAGAACCCCCCCTACGCGCTGAAGGATGGCGGCGCCCACACCGAccaacccaaccaccc                                                                         >  1:112166/1‑144 (MQ=255)

                                                                              ggAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGCAGCCAACGAGCGcc                                                                       >  6:124791/1‑149 (MQ=255)

                                                                                aaGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTg                                                                     <  4:50360/149‑1 (MQ=255)

                                                                                   ccGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCGGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGc                                                                  <  2:58574/149‑1 (MQ=255)

                                                                                    cgcgCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCt                                                                 >  7:284082/1‑149 (MQ=255)

                                                                                           cAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTg                                                           <  2:347373/148‑1 (MQ=255)

                                                                                           cAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGc                                                          >  1:108389/1‑149 (MQ=255)

                                                                                             ggCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCg                                                        <  2:203975/149‑1 (MQ=255)

                                                                                             ggCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCg                                                        <  6:73842/149‑1 (MQ=255)

                                                                                                cATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCg                                                              <  4:354857/140‑1 (MQ=255)

                                                                                                cATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCg                                                              >  3:354857/1‑140 (MQ=255)

                                                                                                  tCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAg                                                   >  7:293503/1‑149 (MQ=255)

                                                                                                    gCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCCTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGCGCTTCGCTGCCGACCAGCg                                                 >  4:252073/1‑149 (MQ=255)

                                                                                                      ccGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGt                                                >  1:247108/1‑148 (MQ=255)

                                                                                                       cgcgAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATggc                                                                                                      >  5:107954/1‑93 (MQ=255)

                                                                                                       cgcgAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATggc                                                                                                      <  6:107954/93‑1 (MQ=255)

                                                                                                       cgcgAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGa                                              <  1:7875/149‑1 (MQ=255)

                                                                                                                cGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTgg                                                                    >  8:189252/1‑118 (MQ=255)

                                                                                                                cGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTgg                                                                    <  7:189252/118‑1 (MQ=255)

                                                                                                                 gTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGCCata                                    <  4:194543/149‑1 (MQ=255)

                                                                                                                  tGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGCCatat                                   >  8:138020/1‑149 (MQ=255)

                                                                                                                    gggTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCt                                                            <  6:10704/122‑1 (MQ=255)

                                                                                                                      gtgtTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGCCATAtgcct                               <  3:266961/149‑1 (MQ=255)

                                                                                                                         ttgttgtACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGCCATAtgcctgc                             <  4:166617/148‑1 (MQ=255)

                                                                                                                            ttgtACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCTTCGAAGGCAACGAGCGCCTGGTCTTCGCTGCCGACCAGCGTGACTACGCCATATGCCTGCCGAt                         <  3:252073/149‑1 (MQ=255)

                                                                                                                            ttgtACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGCCATATGCCTGCCGa                          <  7:338696/148‑1 (MQ=255)

                                                                                                                            ttgtACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGCCATATGCCTGCCGAt                         >  7:61014/1‑149 (MQ=255)

                                                                                                                                 ccTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGCCATATGCCTGCCGATGCTc                     <  1:102688/148‑1 (MQ=255)

                                                                                                                                 ccTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGCCATATGCCTGCCGATGCTc                     <  5:254044/148‑1 (MQ=255)

                                                                                                                                 ccTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGCCATATGCCTGCCGATGCTCg                    <  2:324399/149‑1 (MQ=255)

                                                                                                                                  cTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCTGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGCCATATGCCTGCCGATGCTc                     >  2:41310/1‑147 (MQ=255)

                                                                                                                                    gCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGCCATATGCCTGCCGATGCTCGAg                  >  2:94884/1‑148 (MQ=255)

                                                                                                                                     cGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCg                                                              >  5:235829/1‑103 (MQ=255)

                                                                                                                                     cGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCg                                                              <  6:235829/103‑1 (MQ=255)

                                                                                                                                     cGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGCCATATGCCTGCCGATGCTCGAg                  <  8:293503/147‑1 (MQ=255)

                                                                                                                                      gTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGCCATATGCCCGCCGATGCTCGAGCAc               <  4:232317/149‑1 (MQ=255)

                                                                                                                                       tCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGCCATATGCCTGCCGATGCTCGAGCAcc              >  7:8783/1‑149 (MQ=255)

                                                                                                                                        cAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGCCATATGCCTGCCGATGCTCGAGCACCt             >  6:360651/1‑149 (MQ=255)

                                                                                                                                          g‑‑gAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGCCATATGCCTGCCGATGCTCGAGCACCTgg           >  8:349189/1‑149 (MQ=255)

                                                                                                                                              aaGGCGGTGGCATCGGCGTGCTGAACAAGATGCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGCCATATGCCTGCCGATGCTCGAGCACCTGGGc         <  2:124662/149‑1 (MQ=255)

                                                                                                                                              aaGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGCCATATGCCTGCCGATGCTCGAGCACCTGGGc         <  2:108389/149‑1 (MQ=255)

                                                                                                                                                 gCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGCCATATGCCTGCCGATGCTCGAGCACCTGGGCGTc      <  6:68779/149‑1 (MQ=255)

                                                                                                                                                 gCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGCCATATGCCTGCCGATGCTCGAGCACCTGGGCGTc      <  7:234668/149‑1 (MQ=255)

                                                                                                                                                   cGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGCCATATGCCTGCCGATGCTCGAGCACCTGGGCGTCaa    <  8:342063/149‑1 (MQ=255)

                                                                                                                                                     tGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGCCATATGCCTGCCGATGCTCGAGCACCTGGGCGTCAAGt  >  2:59639/1‑149 (MQ=255)

GGTACTGGGGCGCCTGCATTCCGAGTGCCTGACCGGCGATGCACTGTTCAGCCAGCGCTGCGACTGCGGCTCGCAGTTGGAAGCCGCGCTGCAGGCCATCGCCCGCGAAGGCCGTGGGGTGTTGTTGTACCTGCGTCAG‑‑GAAGGCCGTGGCATCGGCCTGCTGAACAAGATCCGCGCCTACGAGCTGCAGGATGGCGGCGCCGACACCGTCGAAGCCAACGAGCGCCTGGGCTTCGCTGCCGACCAGCGTGACTACGCCATATGCCTGCCGATGCTCGAGCACCTGGGCGTCAAGT  >  NC_002947/607231‑607526

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