BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,321,518	T→G	5.3%	K407Q (AAG→CAG) ‡	PP_1154 ←	sensory box protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,321,518	0	T	G	5.3%	400.1 / 9.1	123	K407Q (AAG→CAG) ‡	PP_1154	sensory box protein
Reads supporting (aligned to +/- strand): ref base T (59/57); new base G (2/5); total (61/62)
Fisher's exact test for biased strand distribution p-value = 4.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 137 total aligned reads displayed.

GCGCAACAAGGCCTGCTCCAGGTGCCGGGCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCAGCCAAGGTGACAGCGGCGTGGTATCGCAG  >  NC_002947/1321369‑1321666

gcgcAACAAGGCCTGCTCCAGGTGCCGGGCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCtt                                                                                                                                                      <  1:5121253/150‑1 (MQ=255)

gcgcAACAAGGCCTGCTCCAGGTGCCGGGCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCtt                                                                                                                                                      <  1:3935141/150‑1 (MQ=255)

gcgcAACAAGGCCTGCTCCAGGTGCCGGGCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCtt                                                                                                                                                      <  2:127016/150‑1 (MQ=255)

gcgcAACAAGGCCTGCTCCAGGTGCCGGGCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCtt                                                                                                                                                      >  2:1656727/1‑150 (MQ=255)

gcgcAACAAGGCCTGCTCCAGGTGCCGGGCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCAGGAAACCGATGGCCtt                                                                                                                                                      >  1:2426597/1‑150 (MQ=255)

gcgcAACAAGGCCTGCTCCAGGTGCCGGGCGAGCTTGTTGCGGTGGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCtt                                                                                                                                                      <  2:4535383/150‑1 (MQ=255)

  gcAACAAGGCCTGCTCCAGGTGCCGGGCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGc                                                                                                                                                    <  2:3794332/150‑1 (MQ=255)

    aacaaGGCCTGCTCCAGGTGCCGGGCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGccg                                                                                                                                                  >  2:4816363/1‑149 (MQ=255)

    aacaaGGCCTGCTCCAGGTGCCGGGCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCa                                                                                                                                                  >  2:4816341/1‑150 (MQ=255)

    aacaaGGCCTGCCCCAGGTGCCGGGCGAGCTGGTTGCGGTTGGCCAGTCCGGTCAGCGGGTCGTCCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGGGATCTCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCa                                                                                                                                                  <  2:969549/150‑1 (MQ=255)

     acaaGGCCTGCTCCAGGTGCCGGGCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGTAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCAt                                                                                                                                                 >  1:809066/1‑150 (MQ=255)

       aaGGCCTGCTCCAGGTGCCGGGCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCAAGGCCATcg                                                                                                                                               <  2:2536610/150‑1 (MQ=255)

        aGGCCTGCTCCAGGTGCCGGGCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCAAGGCCATcgg                                                                                                                                              <  1:3986925/150‑2 (MQ=255)

           ccTGCTCCAGGTGCCGGGCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGcc                                                                                                                                           <  2:658545/150‑1 (MQ=255)

           ccTGCTCCAGGTGCCGGGCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGcc                                                                                                                                           >  2:1897673/1‑150 (MQ=255)

            cTGCTCCAGGTGCCGGGCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCa                                                                                                                                          <  2:1965858/150‑1 (MQ=255)

            cTGCTCCAGGTGCCGGGCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCAAGGCCATcggtttc                                                                                                                                          <  1:2752459/150‑6 (MQ=255)

            cTGCTCCAGGGGCCGGGCGAGCTTGTTGCGGTTGGCCAGTCCGGACAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCAAGGCCATcggtttc                                                                                                                                          <  1:194282/150‑6 (MQ=255)

             tGCTCCAGGTGCCGGGCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTACAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAg                                                                                                                                         >  2:4326850/1‑150 (MQ=255)

                         cGGGCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCgcgggcg                                                                                                                             <  2:2137488/150‑1 (MQ=255)

                         cGGGCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCgcgggcg                                                                                                                             >  2:5028886/1‑150 (MQ=255)

                          gggCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCgcgggcgg                                                                                                                            >  1:2313987/1‑150 (MQ=255)

                          gggCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCgcgggcgg                                                                                                                            >  2:3890699/1‑150 (MQ=255)

                           ggCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTTCCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGa                                                                                                                           <  1:4383042/150‑1 (MQ=255)

                            gCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAGCCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGat                                                                                                                          <  2:2011762/150‑1 (MQ=255)

                                 cTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCg                                                                                                                     >  1:1183054/1‑150 (MQ=255)

                                  ttgttgCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGTCTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCgg                                                                                                                    >  2:4254318/1‑150 (MQ=255)

                                  ttgttgCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCgg                                                                                                                    <  2:1702154/150‑1 (MQ=255)

                                  ttgttgCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCgg                                                                                                                    <  2:1702035/150‑1 (MQ=255)

                                  ttgttgCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCgg                                                                                                                    >  2:4252063/1‑150 (MQ=255)

                                    gttgCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCa                                                                                                                  <  1:3158544/150‑1 (MQ=255)

                                      tgCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGt                                                                                                                >  1:3850068/1‑150 (MQ=255)

                                      tgCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGt                                                                                                                >  1:624437/1‑150 (MQ=255)

                                       gCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTa                                                                                                               <  2:1537883/150‑1 (MQ=255)

                                       gCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTa                                                                                                               <  1:1663908/150‑1 (MQ=255)

                                        cGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAg                                                                                                              <  1:1115865/150‑1 (MQ=255)

                                        cGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAg                                                                                                              <  1:1115784/150‑1 (MQ=255)

                                                cAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCTCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCTGTACACCGGAAAACGATTGCCTTGCCATCGCGCCAGATAGCTCGGGCGGATATCCGGCAGTAGCGGTTCTg                                                                                                      <  1:135368/150‑1 (MQ=255)

                                                cAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTg                                                                                                      <  1:4966583/150‑1 (MQ=255)

                                                         cAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTg                                                                                             <  2:2637085/150‑1 (MQ=255)

                                                         cAGTGGGTCGTACAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTg                                                                                             <  1:4326850/150‑1 (MQ=255)

                                                           gTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGt                                                                                           <  2:2482591/150‑1 (MQ=255)

                                                              ggTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTcgc                                                                                        >  2:2757206/1‑150 (MQ=255)

                                                              ggTCGTGCAGCGAAAGGCGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTcgc                                                                                        >  2:3923707/1‑150 (MQ=255)

                                                              ggTCGTGCAGCGAAAGGCGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTcgc                                                                                        >  2:3924105/1‑150 (MQ=255)

                                                               gTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTcgcg                                                                                       >  2:2990818/1‑150 (MQ=255)

                                                                   tGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTaa                                                                                   <  1:4708935/150‑1 (MQ=255)

                                                                   tGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTaa                                                                                   <  2:5067930/150‑1 (MQ=255)

                                                                     cAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTGCTGGCCATGCAGGTCGCGGTAAGt                                                                                 >  2:3391910/1‑150 (MQ=255)

                                                                       gCGAAAGGTGGTGTATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCGATTCTGGTCGCGGTAAGTAc                                                                               >  2:4073541/1‑150 (MQ=255)

                                                                        cGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACa                                                                              >  2:5048447/1‑150 (MQ=255)

                                                                        cGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACa                                                                              >  2:5047285/1‑150 (MQ=255)

                                                                        cGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACa                                                                              <  2:3248375/150‑1 (MQ=255)

                                                                         gAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAg                                                                             >  1:165902/1‑150 (MQ=255)

                                                                                gCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGt                                                                      <  1:1767641/150‑1 (MQ=255)

                                                                                 cTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGtt                                                                     <  2:4590125/150‑1 (MQ=255)

                                                                                    aaTGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGc                                                                  <  1:3649680/150‑1 (MQ=255)

                                                                                    aaTGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGc                                                                  <  2:2679059/150‑1 (MQ=255)

                                                                                    aaTGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGc                                                                  <  1:1897673/150‑1 (MQ=255)

                                                                                           gcgTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATtgc                                                           >  1:2522689/1‑150 (MQ=255)

                                                                                                   gTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGcggctc                                                   <  1:788964/150‑1 (MQ=255)

                                                                                                     cGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGcggctcgg                                                 <  2:4337871/150‑1 (MQ=255)

                                                                                                     cGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGcggctcgg                                                 <  1:3421385/150‑1 (MQ=255)

                                                                                                       aCTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGcggctcggct                                               <  1:2702507/150‑1 (MQ=255)

                                                                                                       aCTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGcggctcggct                                               <  1:4398363/150‑1 (MQ=255)

                                                                                                        cTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTg                                              <  2:2872049/150‑1 (MQ=255)

                                                                                                          tcgtcgGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCa                                            <  2:233824/150‑1 (MQ=255)

                                                                                                            gtcgGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGt                                          <  1:9680/150‑1 (MQ=255)

                                                                                                             tcgGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTc                                         >  2:2827411/1‑150 (MQ=255)

                                                                                                                  gATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCCTTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGc                                    >  2:3131036/1‑150 (MQ=255)

                                                                                                                  gATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGc                                    >  2:4675163/1‑150 (MQ=255)

                                                                                                                  gATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGc                                    >  2:1574141/1‑150 (MQ=255)

                                                                                                                   atatCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCg                                   >  1:2065552/1‑150 (MQ=255)

                                                                                                                   atatCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCg                                   >  1:2065550/1‑150 (MQ=255)

                                                                                                                    tatCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGt                                  >  1:3182593/1‑150 (MQ=255)

                                                                                                                    tatCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGt                                  >  1:3182513/1‑150 (MQ=255)

                                                                                                                     atCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTc                                 >  2:5120758/1‑150 (MQ=255)

                                                                                                                       cGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCa                               <  2:652367/150‑1 (MQ=255)

                                                                                                                            gCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCAGCCaa                          <  1:4257018/150‑1 (MQ=255)

                                                                                                                             cGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCAGCCAAg                         <  1:2147752/150‑1 (MQ=255)

                                                                                                                              ggTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCAGCCAAgg                        >  2:1324060/1‑150 (MQ=255)

                                                                                                                              ggTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCAGCCAAgg                        >  1:2416585/1‑150 (MQ=255)

                                                                                                                                  ccccGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCAGCCAAGGTGAc                    >  1:3006470/1‑150 (MQ=255)

                                                                                                                                   cccGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCTGTCAGGGCGTCCAGCCAAGGTGACa                   >  1:2452760/1‑150 (MQ=255)

                                                                                                                                    ccGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTCAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCAGCCAAGGTGACAg                  >  2:4991051/1‑150 (MQ=255)

                                                                                                                                     cGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCAGCCAAGGTGACAgc                 >  2:3392128/1‑150 (MQ=255)

                                                                                                                                     cGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGAAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCAGCCAAGGTGACAgc                 >  1:2456916/1‑150 (MQ=255)

                                                                                                                                        aaaCCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCAGCCAAGGTGACAgcggc              <  1:4254566/150‑1 (MQ=255)

                                                                                                                                          aCCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCAGCCAAGGTGACAGCGGCGt            >  1:601518/1‑150 (MQ=255)

                                                                                                                                           gcGATGGCAAGGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGGTGCGGCTCGGCTGCAGTCAGGGCGTCGAGCCAAGGTGACAGCGGCGTg           >  2:856688/2‑150 (MQ=255)

                                                                                                                                           gcGATGGCAAGGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCAGCCAAGGTGACAGCGGCGTg           >  1:314820/2‑150 (MQ=255)

                                                                                                                                           ccGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCAGCCAAGGTGACAGCGGCGTg           <  1:4456154/150‑1 (MQ=255)

                                                                                                                                            cGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCAGCCAAGGTGACAGCGGCGTgg          >  2:4850559/1‑150 (MQ=255)

                                                                                                                                                ggCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCAGCCAAGGTGACAGCGGCGTGGTATc      <  2:2290373/150‑1 (MQ=255)

                                                                                                                                                 gCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCAGCCAAGGTGACAGCGGCGTGGTATCg     >  2:2097387/1‑150 (MQ=255)

                                                                                                                                                   cTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCAGCCAAGGTGACAGCGGCGTGGTATCGCa   <  1:3078929/150‑1 (MQ=255)

                                                                                                                                                    ttGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCAGCCAAGGTGACAGCGGCGTGGTATCGCAg  >  1:1288376/1‑150 (MQ=255)

                                                                                                                                                    ttGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCAGCCAAGGTGACAGCGGCGTGGTATCGCAg  >  1:1288348/1‑150 (MQ=255)

                                                                                                                                                    ttGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCAGCCAAGGTGACAGCGGCGTGGTATCGCAg  <  1:4269206/150‑1 (MQ=255)

                                                                                                                                                    ttGCCATCGCGCCAGATAGCGCGGGCGGAGATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCAGCCAAGGTGACAGCGGCGTGGTATCGCAg  <  1:4816363/150‑1 (MQ=255)

GCGCAACAAGGCCTGCTCCAGGTGCCGGGCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCAGCCAAGGTGACAGCGGCGTGGTATCGCAG  >  NC_002947/1321369‑1321666

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