BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	433,847	C→G	7.2%	V452V (GTG→GTC)	glcB ←	malate synthase G

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	433,847	0	C	G	7.2%	682.8 / 32.3	208	V452V (GTG→GTC)	glcB	malate synthase G
Reads supporting (aligned to +/- strand): ref base C (95/98); new base G (10/5); total (105/103)
Fisher's exact test for biased strand distribution p-value = 2.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

Only 100 of 215 total aligned reads displayed.

CCGGTGGCCAGGCCAACGTCGACGTTGTTGTTTTCGTAGGCGCCGATCCACTTCTCGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGCCACTTC  >  NC_002947/433698‑433996

 cGGTGGCCAGGCCAACGTCGACGTTGTTGTTTTCGTAGGCGCCGATCCACTTCTCGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAaca                                                                                                                                                      <  1:4869307/150‑1 (MQ=255)

         aggccaACGTCGACGTTGTTGTTTTCGTAGGCGCCGATCCACTTCTCGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACgcgc                                                                                                                                              >  1:1146725/1‑150 (MQ=255)

           gccaACGTCGACGTTGTTGTTTTCGTAGGCGCCGATCCACTTCTCGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTc                                                                                                                                            >  1:2770906/1‑150 (MQ=255)

            ccaACGTCGACGTTGTTGTTTTCGTAGGCGCCGATCCACTTCTCGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTcg                                                                                                                                           <  1:1587019/150‑1 (MQ=255)

             caACGTCGACGTTGTTGTTTTCGTAGGCGCCGATCCACTTCTCGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTcgg                                                                                                                                          >  2:3877409/1‑150 (MQ=255)

               aCGTCGACGTTGTTGTTTTCGTAGGCGCCGATCCACTTCTCGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTcggcg                                                                                                                                        <  2:878499/150‑1 (MQ=255)

                cGTCGACGTTGTTGTTTTCGTAGGCGCCGATCCACTTCTCGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTcggcgg                                                                                                                                       <  1:3614679/150‑1 (MQ=255)

                 gTCGACGTTGTTGTGTTCGTAGGCGCCGATCCACTTCTCGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCggcggc                                                                                                                                      >  2:3380417/1‑150 (MQ=255)

                      cGTTGTTGTTTTCGTAGGCGCCGATCCACTTCTCGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCt                                                                                                                                 <  1:3981639/150‑1 (MQ=255)

                         tgttgttTTCGTAGGCGCCGATCCACTTCTCGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGa                                                                                                                              <  2:321171/150‑1 (MQ=255)

                             gttTTCGTAGGCGCCGTTCGACTTCTCGTTCTTCATGGCTCCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTTGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTTTTGATGAACACCACGCGCTCGGCTGCTGCCTTGATGCa                                                                                                                          <  1:4384162/150‑1 (MQ=255)

                              ttttCGTAGGCGCCGATCCACGTCTCGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAg                                                                                                                         <  1:3521968/150‑1 (MQ=255)

                                    tAGGCGCCGATCCACTTCTCGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTg                                                                                                                   <  1:2654988/150‑1 (MQ=255)

                                     aGGCGCCGATCCACTTCTCGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGa                                                                                                                  <  2:1121044/150‑1 (MQ=255)

                                       gcgcCGATCCACTTCTCGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCtt                                                                                                                                >  1:1079679/1‑134 (MQ=255)

                                       gcgcCGATCCACTTCTCGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCtt                                                                                                                                <  2:1079679/134‑1 (MQ=255)

                                        cgcCGATCCACGTCTCGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGt                                                                                                               <  2:4265317/150‑1 (MQ=255)

                                         gcCGATCCACTTCTCGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGa                                                                                                                  >  1:1001194/1‑146 (MQ=255)

                                         gcCGATCCACTTCTCGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGa                                                                                                                  <  2:1001194/146‑1 (MQ=255)

                                         gcCGATCCACTTCTCGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGtt                                                                                                              <  1:3646897/150‑1 (MQ=255)

                                              tCCACTTCTCGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCg                                                                                                         >  1:2868330/1‑150 (MQ=255)

                                                 aCTTCTCGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTgg                                                                                                      >  2:5251070/1‑150 (MQ=255)

                                                 aCTTCTCGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTgg                                                                                                      >  1:4825029/1‑150 (MQ=255)

                                                 aCTTCTCGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATAGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTgg                                                                                                      <  1:2265192/150‑1 (MQ=255)

                                                      tcGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACgg                                                                                                 >  2:3634091/1‑150 (MQ=255)

                                                         ttcttcATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGa                                                                                                            <  1:1381307/136‑1 (MQ=255)

                                                         ttcttcATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGt                                                                                              >  2:3062710/1‑150 (MQ=255)

                                                         ttattcATGGCACCTGTGCGCACCACGGCGCCGGCTTCCGGCGAGGTGTGGCTTTCGTCGGCAGTGAGGTCAAGGAGGCCGGTGGGGGTGAACACCACGCGCTCGGCGGCGGCCTTGGTGGAGGACGTGAGGTGGGCCGGGGTGGGGCGt                                                                                              >  2:4453984/4‑150 (MQ=255)

                                                           cttcATGTCGCCCTTGCGGATCATCGGGCCAGCTTCCATCGAGGTATGGATTTCGTCGCCAGTGCGGTCAAGGAAGCCGGTATTGATGAAGACGACGCGTTCCTTGGCAGCGCGGATCGCTTCCTTGAGGTTAACCGTGGTGCGGCGTTc                                                                                            <  1:5436363/150‑1 (MQ=255)

                                                              cATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGa                                                                                                                              <  1:970711/113‑1 (MQ=255)

                                                              cATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTc                                                                                         >  1:2964071/1‑150 (MQ=255)

                                                                 ggCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTc                                                                                      <  1:1240329/150‑1 (MQ=255)

                                                                 ggCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTc                                                                                      <  1:1811072/150‑1 (MQ=255)

                                                                  gCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTcc                                                                                     <  2:5243634/150‑1 (MQ=255)

                                                                    aCCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCat                                                                                   >  2:1931543/1‑150 (MQ=255)

                                                                      cTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCatga                                                                                 <  1:706354/150‑1 (MQ=255)

                                                                      cTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCatga                                                                                 <  1:5400323/150‑1 (MQ=255)

                                                                      cTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCatga                                                                                 >  1:447917/1‑150 (MQ=255)

                                                                       tttGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCatgat                                                                                >  2:4944688/1‑150 (MQ=255)

                                                                          gcgcACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGcc                                                                             <  1:4944659/150‑1 (MQ=255)

                                                                            gcACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGa                                                                           <  1:4307324/150‑1 (MQ=255)

                                                                              accacGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGAcc                                                                         >  1:4087711/1‑150 (MQ=255)

                                                                               ccacGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCt                                                                        <  2:1822831/150‑1 (MQ=255)

                                                                                  cGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGa                                                                     <  2:4087711/150‑1 (MQ=255)

                                                                                     cgcCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACgg                                                                  <  1:2028319/150‑1 (MQ=255)

                                                                                      gcCAGCTTCCATCGAGGTATGGATTTCGTCGCCAGTGCGGTCAAGGAAGCCGGTATTGATGAAGACGACGCGTTCCTTGGCAGCGCGGATCGCTTCCTTGAGGTTAACCGTGGTGCGGCGTTc                                                                                            >  2:4385899/1‑123 (MQ=255)

                                                                                       ccGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGgtg                                                                <  1:5251049/150‑1 (MQ=255)

                                                                                        cGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGgtgt                                                               >  2:4939896/1‑150 (MQ=255)

                                                                                         ggCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGtt                                                              <  2:944505/150‑1 (MQ=255)

                                                                                            ttCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGTCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTgcg                                                           >  1:382709/1‑150 (MQ=255)

                                                                                             tCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCCCGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTgcgc                                                          <  2:659807/150‑1 (MQ=255)

                                                                                             tCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTgcgc                                                          <  2:1501567/150‑1 (MQ=255)

                                                                                               cATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCgg                                                        >  2:5415529/1‑150 (MQ=255)

                                                                                                aTCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGc                                                       >  2:2776986/1‑150 (MQ=255)

                                                                                                   gAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATg                                                    >  2:4205425/1‑150 (MQ=255)

                                                                                                    aGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGc                                                   <  2:1379755/150‑1 (MQ=255)

                                                                                                      gtgtGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGcc                                                  >  2:533678/1‑149 (MQ=255)

                                                                                                      gtgtGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGcc                                                  <  1:533678/149‑1 (MQ=255)

                                                                                                       tgtgGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCa                                                                                                                          <  1:2113791/76‑1 (MQ=255)

                                                                                                       tgtgGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCa                                                                                                                          >  2:2113791/1‑76 (MQ=255)

                                                                                                        gtgGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCcag                                               <  2:5089750/150‑1 (MQ=255)

                                                                                                        gtgGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCcag                                               >  2:2492352/1‑150 (MQ=255)

                                                                                                           gATTTCGTCGCCAGTGCGGTCAAGGAAGCCGGTATTGATGAAGACGACGCGTTCCTTGGCAGCGCGGATCGCTTCCTTGAGGTTAACCGTGGTGCGGCGTTCCTCATCCATGATGCCGATTTTCAGTGTGTTGGGCTTCATACCaagcag                                            >  1:160502/1‑150 (MQ=255)

                                                                                                             tttCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAggc                                          >  2:1120062/1‑149 (MQ=255)

                                                                                                              ttCGTCGCCAGTGCGGTCAAGGAAGCCGGTATTGATGAAGACGACGCGTTCCTTGGCAGCGCGGATCGCTTCCTTGAGGTTAACCGTGGTGCGGCGTTc                                                                                            >  1:3288647/1‑99 (MQ=255)

                                                                                                               tCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCt                                        >  1:1189535/1‑150 (MQ=255)

                                                                                                                    gCCAGTGCGGTCAAGGAAGCCGGTATTGATGAAGACGACGCGTTCCTTGGCAGCGCGGATCGCTTCCTTGAGGTTAACCGTGGTGCGGCGTTCCTCATCCATGATGCCGATTTTCAGTGTGTTGGGCTTCATACCAAGCAGCTCTTCggt                                   >  2:2345917/1‑148 (MQ=255)

                                                                                                                     ccAGTGCGGTCAAGGAAGCCGGTATTGATGAAGACGACGCGTTCCTTGGCAGCGCGGATCGCTTCCTTGAGGTTAACCGTGGTGCGGCGTTCCTCATCCATGATGCCGATTTTCAGTGTGTTGGGCTTCATACCAAGCAGCTCTTCggtg                                  >  1:1863008/1‑147 (MQ=255)

                                                                                                                      cAGTGCGGTCAAGGAAGCCGGTATTGATGAAGACGACGCGTTCCTTGGCAGCGCGGATCGCTTCCTTGAGGTTAACCGTGGTGCGGCGTTCCTCATCCATGATGCCGATTTTCAGTGTGTTGGGCTTCATACCAAGCAGCTCTTCggtgc                                 >  1:4218319/1‑146 (MQ=255)

                                                                                                                        gTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTgg                               <  1:1469478/150‑1 (MQ=255)

                                                                                                                        gTGCGGTCAAGGAAGCCGGGGTTGCTGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTgg                               <  1:2581345/150‑1 (MQ=255)

                                                                                                                          gCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCt                             <  1:1451827/150‑1 (MQ=255)

                                                                                                                          gCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCt                             <  1:3478712/150‑1 (MQ=255)

                                                                                                                          gCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCt                             <  1:169951/150‑1 (MQ=255)

                                                                                                                            ggTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGa                           <  1:1234482/150‑1 (MQ=255)

                                                                                                                              tCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAg                         <  2:1705585/150‑1 (MQ=255)

                                                                                                                               cAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGCCTTCGACCTGGCTGAAGa                        <  1:1120062/150‑1 (MQ=255)

                                                                                                                                aaGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCatgat                                                                                <  2:3813256/93‑1 (MQ=255)

                                                                                                                                aaGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGAt                       >  1:576303/1‑150 (MQ=255)

                                                                                                                                  ggAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGAtct                     >  2:796865/1‑150 (MQ=255)

                                                                                                                                  ggAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGAtct                     >  1:2716422/1‑150 (MQ=255)

                                                                                                                                   gAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGAtctc                    >  1:2410380/1‑150 (MQ=255)

                                                                                                                                      gCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGt                 <  2:2754275/150‑1 (MQ=255)

                                                                                                                                       ccGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTg                <  1:4970501/150‑1 (MQ=255)

                                                                                                                                       ccGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTg                <  1:4581592/150‑1 (MQ=255)

                                                                                                                                        cGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACg                                                                   >  2:1768523/1‑98 (MQ=255)

                                                                                                                                        cGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTgg               <  2:4762382/150‑1 (MQ=255)

                                                                                                                                         ggTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGc              <  1:4004692/150‑1 (MQ=255)

                                                                                                                                           tgtTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGa            <  1:3897818/150‑1 (MQ=255)

                                                                                                                                               gatgaACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGc        <  1:1058041/150‑1 (MQ=255)

                                                                                                                                               gatgaACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGc        >  2:3656827/1‑150 (MQ=255)

                                                                                                                                               gatgaACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGc        >  2:2810356/1‑150 (MQ=255)

                                                                                                                                               gatgaACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGc        <  1:3064858/150‑1 (MQ=255)

                                                                                                                                                atgaACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGcc       >  2:4730033/1‑150 (MQ=255)

                                                                                                                                                 tgaACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGCCa      >  2:1523162/1‑150 (MQ=255)

                                                                                                                                                 tgaACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGCCa      >  1:5280573/1‑150 (MQ=255)

                                                                                                                                                 tgaACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGCCa      >  1:4721370/1‑150 (MQ=255)

                                                                                                                                                  gaACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGCCAc     >  2:559144/1‑150 (MQ=255)

                                                                                                                                                    acacCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGTTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGCCACtt   >  2:1689516/1‑150 (MQ=255)

                                                                                                                                                    acacCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGCCACtt   >  2:969283/1‑150 (MQ=255)

CCGGTGGCCAGGCCAACGTCGACGTTGTTGTTTTCGTAGGCGCCGATCCACTTCTCGTTCTTCATGGCACCTTTGCGCACCACGGCGCCGGCTTCCATCGAGGTGTGGATTTCATCGCCAGTGCGGTCAAGGAAGCCGGTGTTGATGAACACCACGCGCTCGGCGGCTGCCTTGATGCAGGACTTGAGGTTGACCGTGGTACGGCGTTCCTCGTCCATGATGCCGACCTTGACGGTGTTGCGCGGCATGCCCAGCAGGTCTTCGACCTGGCTGAAGATCTCGGTGGCGAAGGCCACTTC  >  NC_002947/433698‑433996

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