BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,661,519	C→A	8.6%	A50D (GCC→GAC)	PP_1456 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,661,519	0	C	A	8.6%	489.4 / 27.0	151	A50D (GCC→GAC)	PP_1456	membrane protein
Reads supporting (aligned to +/- strand): ref base C (77/61); new base A (9/4); total (86/65)
Fisher's exact test for biased strand distribution p-value = 3.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

Only 100 of 155 total aligned reads displayed.

ATGTTCTATCTGATTGCCTACATCAGCAGCGTAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGACCGCCT  >  NC_002947/1661371‑1661666

aTGTTCTATCTGATTGCCTACATCAGCAGCGTAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGcc                                                                                                                                                    <  1:5403910/150‑1 (MQ=255)

   ttCTATCTGATTGCCTACATCAGCAGCGTAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTg                                                                                                                                                 >  1:1020990/1‑150 (MQ=255)

      tATCTGATTGCCTACATCAGCAGCGTAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTc                                                                                                                                              >  2:3216431/1‑150 (MQ=255)

            aTTGCCTACATCAGCAGCGTAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATg                                                                                                                                        >  1:2656258/1‑150 (MQ=255)

                ccTACATCAGCAGCGTAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTg                                                                                                                                                 >  2:4632375/1‑137 (MQ=255)

                ccTACATCAGCAGCGTAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTATGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCtgg                                                                                                                                    >  1:5304235/1‑150 (MQ=255)

                 cTACATCAGCAGCGTAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCtggt                                                                                                                                   >  2:1213121/1‑150 (MQ=255)

                 cTACATCAGCAGCGTAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCtggt                                                                                                                                   >  1:875839/1‑150 (MQ=255)

                       cagcagCGTAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACt                                                                                                                             >  1:5309218/1‑150 (MQ=255)

                       cagcagCGTAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACt                                                                                                                             >  2:4863341/1‑150 (MQ=255)

                          cagcGTAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTGGTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGt                                                                                                                          >  2:3926767/1‑150 (MQ=255)

                           agcGTAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTg                                                                                                                         <  1:1416442/150‑1 (MQ=255)

                           agcGTAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTg                                                                                                                         <  1:2283823/150‑1 (MQ=255)

                               tAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGt                                                                                                                     >  1:124812/1‑150 (MQ=255)

                               tAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGt                                                                                                                     >  1:124490/1‑150 (MQ=255)

                                aGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTc                                                                                                                    <  2:124812/150‑1 (MQ=255)

                                aGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTc                                                                                                                    <  1:4573632/150‑1 (MQ=255)

                                  tGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCt                                                                                                                  <  1:671016/150‑1 (MQ=255)

                                  tGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTGTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCt                                                                                                                  >  2:687085/1‑150 (MQ=255)

                                   gCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTc                                                                                                                                              >  2:3372655/1‑121 (MQ=255)

                                   gCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTc                                                                                                                                              <  1:3372655/121‑1 (MQ=255)

                                    cTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTAt                                                                                                                <  1:4714536/150‑1 (MQ=255)

                                       aTCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTc                                                                                                             <  1:4863341/150‑1 (MQ=255)

                                       aTCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTc                                                                                                             <  1:4863427/150‑1 (MQ=255)

                                       aTCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTc                                                                                                             >  1:4240404/1‑150 (MQ=255)

                                        tCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCa                                                                                                            >  1:2811304/1‑150 (MQ=255)

                                        tCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCa                                                                                                            >  2:1251811/1‑150 (MQ=255)

                                          aaCTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCAcc                                                                                                          >  1:213016/1‑150 (MQ=255)

                                             tACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCg                                                                                                       <  2:3937516/150‑1 (MQ=255)

                                             tACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGACCAGGGCGCCATga                                                                                                                                       <  1:4417108/118‑2 (MQ=255)

                                             tACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGACCAGGGCGCCATga                                                                                                                                       >  2:4417108/1‑117 (MQ=255)

                                                   ttCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCg                                                                                                 >  1:4750272/1‑150 (MQ=255)

                                                   ttCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCg                                                                                                 >  2:2189335/1‑150 (MQ=255)

                                                    tctcCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGt                                                                                                              >  1:3811234/1‑136 (MQ=255)

                                                      tcCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGcc                                                                                              >  2:2398776/1‑150 (MQ=255)

                                                        cAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCAt                                                                                            >  2:5228571/1‑150 (MQ=255)

                                                         aGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATc                                                                                           <  2:4668552/150‑1 (MQ=255)

                                                           tGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGc                                                                                         <  2:1948363/150‑1 (MQ=255)

                                                           tGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGc                                                                                         <  2:1948344/150‑1 (MQ=255)

                                                              cccGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCt                                                                                      >  1:3469218/1‑150 (MQ=255)

                                                              cccGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCt                                                                                      >  2:4838183/1‑150 (MQ=255)

                                                              cccGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGAGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCt                                                                                      <  2:1887845/150‑1 (MQ=255)

                                                               ccGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGATGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTc                                                                                                             >  1:1876208/1‑126 (MQ=255)

                                                               ccGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGATGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTc                                                                                                             <  2:1876208/126‑1 (MQ=255)

                                                                         aCATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCgg                                                                                                <  1:4794678/129‑1 (MQ=255)

                                                                         aCATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCgg                                                                                                >  2:4794678/1‑129 (MQ=255)

                                                                              atcTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCg                                                                      <  2:5309218/150‑1 (MQ=255)

                                                                              atcTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCg                                                                      <  2:5310601/150‑1 (MQ=255)

                                                                               tcTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGt                                                                     <  2:1020990/150‑1 (MQ=255)

                                                                                  ggTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGc                                                                                                                                            >  2:4877614/1‑76 (MQ=255)

                                                                                   gTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGttctt                                                                 <  1:5386156/150‑1 (MQ=255)

                                                                                      tGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCAt                                                              >  1:4900228/1‑150 (MQ=255)

                                                                                       gCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCAtt                                                             >  1:1097370/1‑150 (MQ=255)

                                                                                         cTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTc                                                           >  1:2102040/1‑150 (MQ=255)

                                                                                         cTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTc                                                           <  1:208347/150‑1 (MQ=255)

                                                                                           ggggCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCg                                                         >  1:4933646/1‑150 (MQ=255)

                                                                                            gggCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGa                                                        <  2:3123456/150‑1 (MQ=255)

                                                                                            gggCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGa                                                        <  2:5415818/150‑1 (MQ=255)

                                                                                              gcggcCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGc                                                      <  2:4728433/150‑1 (MQ=255)

                                                                                                 gcCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGa                                                   <  1:4068487/150‑1 (MQ=255)

                                                                                                  ccTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGAt                                                  >  1:4336411/1‑150 (MQ=255)

                                                                                                  ccTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGAt                                                  >  1:1590655/1‑150 (MQ=255)

                                                                                                    tGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCg                                                >  1:1632789/1‑150 (MQ=255)

                                                                                                         ttCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGActgg                                           >  1:3921682/1‑150 (MQ=255)

                                                                                                         ttCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGActgg                                           <  1:2484471/150‑1 (MQ=255)

                                                                                                          tCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGActggc                                          <  1:2006283/150‑1 (MQ=255)

                                                                                                          tCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCGGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGActggc                                          <  2:3823347/150‑1 (MQ=255)

                                                                                                            aTCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCt                                                                                                                       <  1:5366516/71‑1 (MQ=255)

                                                                                                            aTCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGActggctg                                        <  1:1849754/150‑1 (MQ=255)

                                                                                                            aTCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGActggctg                                        <  1:2353245/150‑1 (MQ=255)

                                                                                                            aTCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGActggctg                                        >  2:5251555/1‑150 (MQ=255)

                                                                                                             tCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGActggctgg                                       >  1:2985688/1‑150 (MQ=255)

                                                                                                              ccTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGgt                                      >  1:3112011/1‑150 (MQ=255)

                                                                                                              ccTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGgt                                      >  1:3112033/1‑150 (MQ=255)

                                                                                                                tGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGgtgt                                    >  2:179249/1‑150 (MQ=255)

                                                                                                                  cgcgACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTc                                  <  1:1829109/150‑1 (MQ=255)

                                                                                                                  cgcgACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTc                                  <  1:2381700/150‑1 (MQ=255)

                                                                                                                   gcgACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCa                                 <  1:1814723/150‑1 (MQ=255)

                                                                                                                   gcgACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCa                                 >  1:2965899/1‑150 (MQ=255)

                                                                                                                          ggTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCAc                          <  2:4933646/150‑1 (MQ=255)

                                                                                                                          ggTGCAGACCCGTTGGGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACGGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCAc                          <  2:2656697/150‑1 (MQ=255)

                                                                                                                           gTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCAcc                         <  1:4710136/150‑1 (MQ=255)

                                                                                                                           gTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCAcc                         <  2:3528150/150‑1 (MQ=255)

                                                                                                                            tGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGActgg                                           <  2:2982476/131‑1 (MQ=255)

                                                                                                                               aGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTc                     <  1:4648735/150‑1 (MQ=255)

                                                                                                                                gACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCg                    <  2:422849/150‑1 (MQ=255)

                                                                                                                                  cccGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGc                  >  2:4796736/1‑150 (MQ=255)

                                                                                                                                   ccGTTTTGGTCGTGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCCTCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGcc                 >  2:1641733/1‑150 (MQ=255)

                                                                                                                                   ccGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGcc                 >  2:763047/1‑150 (MQ=255)

                                                                                                                                       caccagcATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTg             >  2:1729580/7‑150 (MQ=255)

                                                                                                                                         ccagcATGGCGACCAGGGCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGc                                                                       >  1:2987401/5‑90 (MQ=38)

                                                                                                                                            tCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTg             >  1:3772372/1‑145 (MQ=255)

                                                                                                                                            tCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTg             <  2:3772372/145‑1 (MQ=255)

                                                                                                                                             cATGGCGACCAGGGCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGtt                                                                    >  1:2690523/1‑89 (MQ=255)

                                                                                                                                             cATGGCGACCAGGGCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGAc       >  1:3395467/1‑150 (MQ=255)

                                                                                                                                               tGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGACCg     >  1:2365849/1‑150 (MQ=255)

                                                                                                                                               tGGCGACCAGGGCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGACCg     >  2:2055497/1‑150 (MQ=255)

                                                                                                                                               tGGCGACCAGGGCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGACCg     >  2:2055602/1‑150 (MQ=255)

                                                                                                                                                  cgcCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGACCGCCt  <  2:2790571/150‑1 (MQ=255)

                                                                                                                                                   gcCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTgc            <  2:3077669/139‑1 (MQ=255)

ATGTTCTATCTGATTGCCTACATCAGCAGCGTAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGACCGCCT  >  NC_002947/1661371‑1661666

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