BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,856,777	C→T	5.2%	A235V (GCC→GTC)	PP_5133 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,856,777	0	C	T	5.2%	295.1 / 3.1	96	A235V (GCC→GTC)	PP_5133	membrane protein
Reads supporting (aligned to +/- strand): ref base C (53/38); new base T (3/2); total (56/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.75e-01

CTCACCGCGCTGGGCATGCTACTGGTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGGTTATGCCTGT  >  NC_002947/5856638‑5856923

cTCACCGCGCTGGGCATGCTACTGGTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCcga                                                                                                                                           >  5:321280/1‑149 (MQ=255)

 tCACCGCGCTGGGCATGCTACTGGTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCcgac                                                                                                                                          >  3:191895/1‑149 (MQ=255)

            ggCATGCTACTGGTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGAt                                                                                                                               >  5:58930/1‑149 (MQ=255)

              cATGCTACTGGTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATcc                                                                                                                             >  7:184803/1‑149 (MQ=255)

                   tACTGGTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCgcgc                                                                                                                         >  3:275903/1‑148 (MQ=255)

                   tACTGGTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCgcgc                                                                                                                         >  6:282781/1‑148 (MQ=255)

                     cTGGTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAg                                                                                                                      >  6:209353/1‑149 (MQ=255)

                      tGGTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAgc                                                                                                                     >  6:110023/1‑149 (MQ=255)

                        gTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAg                                                                                                                                      >  5:36865/1‑130 (MQ=255)

                        gTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAg                                                                                                                                      <  6:36865/130‑1 (MQ=255)

                        gTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAgcgg                                                                                                                   >  1:173694/1‑149 (MQ=255)

                         tCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAgcggc                                                                                                                  <  1:286965/149‑1 (MQ=255)

                           ttGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGa                                                                                                                >  6:285586/1‑149 (MQ=255)

                           ttGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGa                                                                                                                >  6:198795/1‑149 (MQ=255)

                           ttGCTGGCGACCCGGCACTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGa                                                                                                                >  1:136064/1‑149 (MQ=255)

                                   gACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGc                                                                                                        <  5:285586/149‑1 (MQ=255)

                                      ccGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAg                                                                                                     >  5:91489/1‑149 (MQ=255)

                                      ccGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAg                                                                                                     >  1:176706/1‑149 (MQ=255)

                                      ccGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAg                                                                                                     >  8:175963/1‑149 (MQ=255)

                                         gcgcTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCCCCAGCGGCGACAGCCTGCCAGccc                                                                                                  <  5:110023/149‑1 (MQ=255)

                                           gcTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTg                                                                                                >  3:145838/1‑149 (MQ=255)

                                           gcTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTg                                                                                                <  8:184803/149‑1 (MQ=255)

                                             tATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGa                                                                                              >  5:202547/1‑149 (MQ=255)

                                              aTGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCcgc                                                                                                                           >  4:145478/1‑119 (MQ=255)

                                              aTGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCcgc                                                                                                                           <  3:145478/119‑1 (MQ=255)

                                              aTGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGAc                                                                                             >  8:261152/1‑149 (MQ=255)

                                                       tCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGc                                                                                    >  4:14541/1‑149 (MQ=255)

                                                         tGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGc                                                                                                                                               <  3:254964/88‑1 (MQ=255)

                                                         tGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGc                                                                                                                                               >  4:254964/1‑88 (MQ=255)

                                                         tGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTgg                                                                                               >  5:167221/1‑136 (MQ=255)

                                                         tGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTgg                                                                                               <  6:167221/136‑1 (MQ=255)

                                                         tGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGc                                                                                    >  7:243194/1‑147 (MQ=255)

                                                           ggAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGcc                                                                                   <  4:159945/146‑1 (MQ=255)

                                                           ggAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAgg                                                                                >  4:144690/1‑149 (MQ=255)

                                                            gAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGc                                                                               >  1:118177/1‑149 (MQ=255)

                                                             aaaCCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGcc                                                                                                   <  3:79169/128‑1 (MQ=255)

                                                             aaaCCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGcc                                                                                                   >  4:79169/1‑128 (MQ=255)

                                                               accaccTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTg                                                                            <  5:198795/149‑1 (MQ=255)

                                                               accaccTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTg                                                                            <  3:144690/149‑1 (MQ=255)

                                                                ccaccTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCTAGCACTGGGCGCCCTCCCGTCGTTGCTGTGCTTCGCCTTGCCCGACGAAGCAATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTgg                                                                           <  8:209939/149‑1 (MQ=255)

                                                                ccaccTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCGCGCCAGGCCTgg                                                                           <  7:245563/149‑1 (MQ=255)

                                                                ccaccTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTgg                                                                           <  4:275903/149‑1 (MQ=255)

                                                                  accTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGc                                                                                                            >  5:184885/1‑114 (MQ=255)

                                                                  accTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGc                                                                                                            <  6:184885/114‑1 (MQ=255)

                                                                   ccTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGc                                                                         >  5:37166/1‑148 (MQ=255)

                                                                    cTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCa                                                                       <  6:91489/149‑1 (MQ=255)

                                                                     ttGCTCGCCGCCGACCCGTTCATTAACCTGACCCAAGCACTGGGCTCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAg                                                                      >  7:266291/1‑149 (MQ=255)

                                                                            ccgccgACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAgctggctg                                                               >  3:54739/1‑149 (MQ=255)

                                                                                cgACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCg                                                           >  1:31511/1‑149 (MQ=255)

                                                                                 gACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCgg                                                          >  3:86371/1‑149 (MQ=255)

                                                                                    ccGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCgt                                                       >  2:65700/1‑149 (MQ=255)

                                                                                     cGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCgtg                                                      <  7:175963/149‑1 (MQ=255)

                                                                                      gTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCgtgg                                                     >  1:85954/1‑149 (MQ=255)

                                                                                      gTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCgtgg                                                     >  6:271427/1‑149 (MQ=255)

                                                                                        tcattcaCCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGCCGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCgtggtg                                                   >  2:306274/1‑149 (MQ=255)

                                                                                          attcaCCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCt                                                 >  2:140926/1‑149 (MQ=255)

                                                                                             caCCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTgg                                               <  5:271427/148‑1 (MQ=255)

                                                                                             caCCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTgg                                               <  2:176706/148‑1 (MQ=255)

                                                                                              aCCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAgcggc                                                                                                                  <  4:306976/80‑1 (MQ=255)

                                                                                              aCCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAgcggc                                                                                                                  >  3:306976/1‑80 (MQ=255)

                                                                                              aCCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGt                                              <  4:145838/148‑1 (MQ=255)

                                                                                               gcTGACCCAAGCCCTCGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCGCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGtat                                            <  8:243194/148‑1 (MQ=255)

                                                                                                cTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGtata                                           >  3:170716/1‑149 (MQ=255)

                                                                                                     ccAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTgg                                               <  3:18828/140‑1 (MQ=255)

                                                                                                     ccAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTgg                                               >  4:18828/1‑140 (MQ=255)

                                                                                                     ccAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGcc                                      >  6:4033/1‑149 (MQ=255)

                                                                                                     ccAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGcc                                      >  1:148353/1‑149 (MQ=255)

                                                                                                            cTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGcc                               >  5:310879/1‑149 (MQ=255)

                                                                                                             tGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGTCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGccc                              >  7:26138/1‑149 (MQ=255)

                                                                                                               ggCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGcccc                             <  7:235479/148‑1 (MQ=255)

                                                                                                               ggCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCg                            >  1:105406/1‑149 (MQ=255)

                                                                                                               ggCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCg                            >  1:101544/1‑149 (MQ=255)

                                                                                                                gcgcCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGc                           <  2:148353/149‑1 (MQ=255)

                                                                                                                gcgcCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGc                           <  2:31511/149‑1 (MQ=255)

                                                                                                                  gcCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCgtg                                                      <  3:178839/120‑1 (MQ=255)

                                                                                                                  gcCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCgtg                                                      >  4:178839/1‑120 (MQ=255)

                                                                                                                        ccGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCctgctgct                   >  3:23552/1‑149 (MQ=255)

                                                                                                                         cgtcgtTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCctgctgctg                  >  7:139880/1‑149 (MQ=255)

                                                                                                                             gtTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGcc               <  2:101544/148‑1 (MQ=255)

                                                                                                                             gtTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGcc               <  2:105406/148‑1 (MQ=255)

                                                                                                                              ttGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCgg             <  4:170716/149‑1 (MQ=255)

                                                                                                                              ttGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCgg             <  2:118177/149‑1 (MQ=255)

                                                                                                                              ttGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCgg             <  4:86371/149‑1 (MQ=255)

                                                                                                                                gCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGGt           >  2:91859/1‑149 (MQ=255)

                                                                                                                                  tGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGCCCTGCTGCCCCGCCTGCTGCTGGCCggg            <  2:66171/146‑1 (MQ=255)

                                                                                                                                  tGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGCCCTGCTGCCCCGCCTGCTGCTGGCCggg            >  1:66171/1‑146 (MQ=255)

                                                                                                                                   tggCTTCGTCGGGCACGGCGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCt                                                                             >  2:230847/2‑80 (MQ=38)

                                                                                                                                   tggCTTCGTCGGGCACGGCGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCt                                                                             <  1:230847/79‑1 (MQ=38)

                                                                                                                                    ggCTTCGTCGGGCACGGCGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGGTTATg       >  4:32157/1‑149 (MQ=255)

                                                                                                                                    ggCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGGTTa         <  7:261152/147‑1 (MQ=255)

                                                                                                                                      cTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGGTTATGc      <  7:176999/148‑1 (MQ=255)

                                                                                                                                      cTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGGTTATGc      >  8:176999/1‑148 (MQ=255)

                                                                                                                                       ttCGTCGGGCACGGCGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGGTTATGc      >  6:315940/1‑147 (MQ=255)

                                                                                                                                       ttCGTCGGGCACGGCGAAGCCATGATCCGCGCCACCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGGTTATGc      <  5:315940/147‑1 (MQ=255)

                                                                                                                                          gCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGGTTATGCCtgt  <  2:173694/148‑1 (MQ=255)

                                                                                                                                          gCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGGTTATGCCtgt  <  1:91859/148‑1 (MQ=255)

CTCACCGCGCTGGGCATGCTACTGGTCTTGCTGGCGACCCGGCGCTATGGCTTTGTCTGGGAAACCACCTTGCTCGCCGCCGACCCGTTCATTCACCTGACCCAAGCACTGGGCGCCCTCCCGTCGTTGCTGGGCTTCGCCGTGCCCGACGAAGCCATGATCCGCGCCAGCGGCGACAGCCTGCCAGCCCTGGACGTGGCCCGCCAGGCCTGGGCCAGCTGGCTGCTCGGCGTGGTGCTGGTATACGGCCTGCTGCCCCGCCTGCTGCTGGCCGGGTTATGCCTGT  >  NC_002947/5856638‑5856923

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: `‑`