BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,163,456	C→G	9.3%	A158G (GCG→GGG)	PP_3661 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,163,456	0	C	G	9.3%	489.0 / 6.9	180	A158G (GCG→GGG)	PP_3661	membrane protein
Reads supporting (aligned to +/- strand): ref base C (79/82); new base G (12/7); total (91/89)
Fisher's exact test for biased strand distribution p-value = 3.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 188 total aligned reads displayed.

ATCGTGTCGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCAC  >  NC_002947/4163308‑4163603

 tCGTGTCGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGgcgc                                                                                                                                                    >  2:3425540/1‑150 (MQ=255)

 tCGTGTCGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGgcgc                                                                                                                                                    >  1:2133090/1‑150 (MQ=255)

 tCGTGTCGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGgcgc                                                                                                                                                    >  1:3914610/1‑150 (MQ=255)

    tgtCGGTGCTGGTAGTCACGAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGGTCGACCGCGATACCGCCTTGCTGTCCGCCGCCGGCAGCGCCGTCTGGGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCTgg                                                                                                                                                 >  2:1943967/1‑149 (MQ=255)

     gtCGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTgcg                                                                                                                                                <  2:3444678/150‑1 (MQ=255)

        ggTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCggcgcggcgccg                                                                                                                                             >  2:4171334/1‑148 (MQ=255)

          tGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAgcg                                                                                                                                           >  2:3629161/1‑150 (MQ=255)

          tGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAgcc                                                                                                                                           >  2:4774286/1‑149 (MQ=255)

          tGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCcgcg                                                                                                                                           >  2:1158893/1‑150 (MQ=255)

            cTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCTGCGCAGCGcc                                                                                                                                         >  2:5848350/1‑150 (MQ=255)

             tGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGGTGCGCAGCGccc                                                                                                                                        >  1:267916/1‑150 (MQ=255)

             tGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGccc                                                                                                                                        <  2:1603436/150‑1 (MQ=255)

             tGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGccc                                                                                                                                        <  2:2133090/150‑1 (MQ=255)

             tGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGGTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCAGTCGACTCGGCGCGGCGCCGCGccc                                                                                                                                        >  2:6840431/1‑150 (MQ=255)

               gTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGGTGCGCAGCGccccc                                                                                                                                      <  2:3505479/150‑1 (MQ=255)

               gTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGccccc                                                                                                                                      >  2:850635/1‑150 (MQ=255)

                tAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGcccccc                                                                                                                                     >  2:2489353/1‑150 (MQ=255)

                tAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGcccccc                                                                                                                                     <  2:4944778/150‑1 (MQ=255)

                tAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGcccccc                                                                                                                                     >  2:5255990/1‑150 (MQ=255)

                  gTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGccccccca                                                                                                                                   >  1:6098316/1‑148 (MQ=255)

                  gTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGccccccca                                                                                                                                   >  1:4975442/1‑148 (MQ=255)

                   tCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCTGGGGAGCGcccccccaa                                                                                                                                  >  1:2318050/1‑147 (MQ=255)

                     aCCAGTACCCTGCTGTTGGGCCTGTGGGGTGGCATGAATGTGTTCAAGCTCGACGGCGATACCGCGTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCTTTGGATTGGGGGCTGGGCAGCGCCCCCCATaaa                                                                                                                                <  2:2901717/150‑1 (MQ=255)

                     aCCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCTGCGCACCGcccccccgaaa                                                                                                                                >  1:4498268/1‑145 (MQ=255)

                     aCCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGcccccccaaaa                                                                                                                                >  1:3814296/1‑145 (MQ=255)

                     aCCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCAaaaa                                                                                                                                >  1:3635985/1‑150 (MQ=255)

                        tgTACCCTGCTGTTGGGCGTGTGGTGTGGCTTTAATGTGTTCAAGCTCGACCGGGTTACGGCCTTGCTGACCGCGGCCGGCTGCGCCATCTGGGGCGCGGCCGCGGTGTTGGATTTCGAATGGGCGCTGCGCAGCGCCCCCCATAAAAgc                                                                                                                             <  1:7329640/149‑1 (MQ=255)

                        cgTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCTATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCTGCGCCCCCCATAAAAgc                                                                                                                             <  2:2326689/149‑1 (MQ=255)

                        aGTACCCTGCTGTTCGGCCTGTGGTGTGGAATTAATTTTTTCATGCTCTACCGCTATACCGCCTTGCTGACGGCCGCCGGCAGCGCCATCTGCGGCGCGGCCGCCGTGCTGGCATTCGAATCGGCGGTGCGCAGCGCCCCCCATAAAAgc                                                                                                                             <  2:1380873/150‑1 (MQ=255)

                            cccTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCTGCGCAATCTGCGGCGCCGCCGCGGTGCTGGCTTTCTAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCa                                                                                                                        <  2:1404965/150‑1 (MQ=255)

                              ctgctgATCGGCCTGTGGTGTGGCATGAAAGTGTTCATGCTCGACCGCGATACGGCCTTGCTGACCGCCGCCGGCTGCGCCATCTGGGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCATg                                                                                                                      <  2:3737344/150‑1 (MQ=255)

                              ctgctgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCCGCGCAGACCCCCCCAAAAAAAAG‑CCaag                                                                                                                      >  1:2924944/1‑148 (MQ=255)

                               tgctgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCCCCGCGCCCCCCCAAAAAGCG‑cccggg                                                                                                                     >  2:5372753/1‑146 (MQ=255)

                                  tgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAg                                                                                                                  >  2:2430906/1‑150 (MQ=255)

                                   gATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCCAAAAAGCG‑CCCAGGCAGt                                                                                                                 >  1:4298301/1‑150 (MQ=255)

                                    aTCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCACAAAAACG‑CCAAGGCAGTc                                                                                                                >  1:456809/1‑150 (MQ=255)

                                        gCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAGTCGGCGCGGCGCAGCGCCCCCCCAAAAACCC‑CCCAGGCAGTCggca                                                                                                            >  2:7097632/1‑148 (MQ=255)

                                               gtgtggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGGGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGt                                                                                                     <  2:6877969/150‑1 (MQ=255)

                                                tgtggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTg                                                                                                    <  2:4925751/150‑1 (MQ=255)

                                                  tggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCCCCCCCCCCCAAAAAGCC‑CCCAGGCCGACGGAAGCGGAGTggg                                                                                                  >  1:5739799/1‑148 (MQ=255)

                                                    gCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCCCCCCCCCAAAAAACC‑CCAGGGCCCTCGGGAGCCGAGGGCTGt                                                                                                >  1:3714218/1‑150 (MQ=255)

                                                     cATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGGAGTGCTGtt                                                                                               >  2:5652301/1‑150 (MQ=255)

                                                      aTGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCCAAAAAGCG‑CCCTGGCAGTCGGCAGCGGAGTGCTGTTc                                                                                              >  1:6825209/1‑150 (MQ=255)

                                                      aTGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTc                                                                                              >  2:4236727/1‑150 (MQ=255)

                                                       tGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCGCCCCCCCCCAAAAAACG‑CCATGGCAGTCGGGAGCGTAGTTCTGttgg                                                                                             >  1:898564/1‑148 (MQ=255)

                                                        taatgttttCAAGCTCGACCGCTATACCGCTTTGCTGACCGCCGCCGGCAGCGCCATCGGGGGCGCCGCCGCGGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCgg                                                                                            <  1:5510084/143‑1 (MQ=255)

                                                         aatttgtTCACGCCCGCCCGCGACCCGGCCTTGCGGCCCGCCGCGGGCCGCGCCAGGGGCGGCGCCGCCGCGGGGGGGGATTTCTCCGGGGCGCGGCGCCGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGt                                                                                           <  1:2190671/146‑1 (MQ=255)

                                                         aaaGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAAGCGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGt                                                                                           <  2:6598199/150‑1 (MQ=255)

                                                          aaGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCCTGGCAGTCGGAAGCGTAGTGCTGTTCGGTa                                                                                          >  2:7362587/1‑150 (MQ=255)

                                                              gtTCACGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCTTCCTCACCGGCGCTGCGCAGCGCCCCCCATACACGCG‑CCAGGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCtt                                                                                      <  2:1464239/150‑1 (MQ=255)

                                                              gtTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCACCCACCCCCAGAAAAACG‑CCAAGGCAGTCGGAAGCGTAGTGCAGTTCGGTAccat                                                                                      >  2:266447/1‑148 (MQ=255)

                                                               ttCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCCCCGCCCCCCCGAAAAGCG‑CCCTGGCAGGGGGAAGCGTAGTGGTGGTCGGTGCCGTg                                                                                     >  1:5231814/1‑150 (MQ=255)

                                                                tCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGGCGCAGCGCCCCCCCCAAAAGCG‑CCAAGGCACTCGGAAGCCTAGTGGTGGTGGGGAccgagt                                                                                    >  1:4381733/1‑146 (MQ=255)

                                                                 cAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTc                                                                                   <  2:7409895/150‑1 (MQ=255)

                                                                 cAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCAGAAAAGCG‑CCATGGCAGTCGGAAGCGGAGTGCTGTTAGGTACCTTGTc                                                                                   >  2:6208744/1‑150 (MQ=255)

                                                                 cAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCACCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTc                                                                                   <  2:5110140/150‑1 (MQ=255)

                                                                  aaGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGGGCCGCGCCCCCCCTAAAAGCG‑CCCAGGCAGGCGGAAGCGTAGTGCTGTTGGGTACCTGGTCg                                                                                  >  1:1139982/1‑150 (MQ=255)

                                                                      tCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGTTGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGt                                                                              <  2:3709672/150‑1 (MQ=255)

                                                                      tCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGt                                                                              <  2:5359548/150‑1 (MQ=255)

                                                                                  ccGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGGTGCGCCGCGCCCCCCCGGAAATCGCCCATGGC‑GTAGGAAGGGTTGTGCGGGTGGGTACGGTGTGGGTGGTGCTTTAGCCAc                                                                  >  1:4579514/1‑150 (MQ=255)

                                                                                   cGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGGCGCAGCGCCCCCCCCAAAAAAC‑CCAACGCCGTAGGGAGCGGAGTGGGGTGCGGTGCCGTGTGGGTGTTTCGTTAGCCAca                                                                 >  1:785019/1‑149 (MQ=255)

                                                                                    gCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGGGGCGCCGCGCCCCCCCAAAAAACG‑CCAAGGCAGTCGGAAGCGGAGTGGGGTTGGGTAGCGTGTGGGTGTTGCTTTTACCAcag                                                                >  1:3150630/1‑148 (MQ=255)

                                                                                    gCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGGCGCCGCGCCCCCCCCAAAAACG‑CCCAGGCAGTCGGGAGCGTACTGGTGGTGGGGACCATGTGGATGTTGCTTTTACCAcag                                                                >  1:3150584/1‑148 (MQ=255)

                                                                                     ccTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACAgg                                                               >  2:6121612/1‑150 (MQ=255)

                                                                                      cgtGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCGGGCATTCGACTCGGCGGTGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGc                                                              <  2:4647703/148‑1 (MQ=255)

                                                                                        tGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCGGTGCTGGCTTTCGAATGGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTGGTCGATGTTCCTTTACCCACTGGCGa                                                            <  2:4460377/150‑1 (MQ=255)

                                                                                         gCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCCGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGTCGAt                                                           >  2:6703587/1‑150 (MQ=255)

                                                                                          cTGACCGCCGCCGGCAGCGCCATCTGGGGCGCCGCCGCGGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATc                                                          <  2:1511820/150‑1 (MQ=255)

                                                                                           tGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCAAAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGGTCGGTACCATGTCGATGTTCCTTTAACCACTGGCGATCa                                                         >  2:1289363/1‑150 (MQ=255)

                                                                                               cgccgccgGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCAGGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAAttc                                                     >  2:1221463/1‑148 (MQ=255)

                                                                                                    ccgGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCTTTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCggct                                                <  2:550890/150‑1 (MQ=255)

                                                                                                      ggCTGCGCCTTGTGCGGCGCCGCCGCCGTGCTTTAATTCGAATCGGCGGTGCGCCGCCCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCggctgg                                              <  1:6213527/150‑1 (MQ=255)

                                                                                                       gCAGCGCCTTCTGCGGCGCCGCCGCCGTGCTGGCATTCGAAGCGGCGCTGCGCAGCGCCCCCCATCAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCggctggc                                             <  2:6851200/150‑1 (MQ=255)

                                                                                                        ccgcgcCATCGGCGGCGCCGCCGCCGTGCTGGCATTCGACTCGGCGCGGCGCCGCGCCCCCCATAACAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCggctggct                                            <  2:3113739/148‑1 (MQ=255)

                                                                                                          gcgcCAGCGGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGc                                          <  1:4250760/150‑1 (MQ=255)

                                                                                                               cacggcggcggcgCCGCAGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCAAAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGCTGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTg                                     >  2:5809631/5‑150 (MQ=255)

                                                                                                               cacggcggcggcgCCGCAGTGCTGGCATTCGAATCGGCGCAGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACAGGCGATCAATGCCGGCTGGCTGCATCTg                                     >  2:5810080/5‑150 (MQ=255)

                                                                                                               aTCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCCGCGCCCCCCAAAAAAGCG‑CCATGGCCGTCGGAAGCGTAGTGCTGTTGGGGTCCATGTCGATGTTACTTTAACCACTGGCGATAAAATCCGGCTGGCTGCAACTg                                     >  1:1087008/1‑150 (MQ=255)

                                                                                                                tCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCGGCATCTgg                                    >  2:6029916/1‑150 (MQ=255)

                                                                                                                       cgccgccgcCGTGCTGGCTTTCGAATCGGCGCGGCGCCGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCAt                             <  1:7054345/150‑1 (MQ=255)

                                                                                                                       cgccgccgcCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCAt                             <  2:7222835/150‑1 (MQ=255)

                                                                                                                           gccgccgTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCAAAAAAGCG‑CCAAGGCAGTCGGAAGCGTAGTGCGGTTCGGTAACATGTCGTTGTTCCGTTAACCACAGGCGATCAAATCCGGCTGGCTGCATCGGGATAACATGGgc                         >  1:4126520/1‑150 (MQ=255)

                                                                                                                           gccgccgTGCTGGCATTCGAATCGGCGCGGCGCAGCGCCCCCCAAAAAAGCG‑CCATTGCAGTCGGGAGCGTAGTGCTGTTCGGTACCTTGTCGATTTTACTTTACCCACAGGCGGTCAAATCCCGCAGGCTGCATCTGGATACCCTGGgc                         >  1:6077098/1‑150 (MQ=255)

                                                                                                                            ccgccgTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCAAGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACACACAGGCGATCAAAGCCCGCTGGCTGCATCTGGGTACCATGGgcg                        >  2:270002/1‑150 (MQ=255)

                                                                                                                             cgccgTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCAAAAAAGCG‑CCATGGCAGTAGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTAACCACAGGCGAACAAAGACCGCTGGCTGCAACTGGATACCAAGGgcgc                       >  1:5450163/1‑150 (MQ=255)

                                                                                                                             cgccgTGCTGGCAGTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCGGTTCGGTACCTTGTCGATGTTCCGTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGgcgc                       >  2:3074416/1‑150 (MQ=255)

                                                                                                                             cgccgTGCTGGAATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGgcgc                       <  1:1750474/150‑1 (MQ=255)

                                                                                                                                ggTGCTGGTATTCGAGGCGGGGCTGCGCCGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCgg                    <  1:689843/149‑1 (MQ=255)

                                                                                                                                cgTGCTGGCATTCGATTCGGCGCTGCGCAGCGCCCCCCATAACAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCgg                    <  2:5679698/150‑1 (MQ=255)

                                                                                                                                      ggCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGgctgct              >  2:2254707/1‑150 (MQ=255)

                                                                                                                                        cATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGgctgctgc            <  2:5059128/150‑1 (MQ=255)

                                                                                                                                         aTTCGACTCGGCGCGGCGCAGCGCCCCCCCGAAAAGCG‑CCATGGCAGTCGGAAGCGTAGGGCGGGTCGGTACCGTGTCGATGTTCCGTTACCCACAGGCGATCAAAGCCGGCGGGCGGCATCTGGATACCATGGGCGCCGGgcggctgct           >  2:937915/1‑150 (MQ=255)

                                                                                                                                          ttggAATCGGCGCGGGGCTGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTGACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGctgctgctg          <  2:456809/147‑1 (MQ=255)

                                                                                                                                          ttCGAATCGGCGCGGCCCAGCGCCCCCCAGAAAAGCG‑CCATGGCAGTCGGAACCGAAGTGCTGTTCCGTACCAGGTAGATGTTCCTATAACCCCAGGCGACCAAAGCCGGCTGGCTGAATCAGGATACCAAGGGCGCCGGGctgctgctg          >  2:1451066/1‑150 (MQ=255)

                                                                                                                                             tacggggCGCTGCGCCGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGggc       <  1:3271739/145‑1 (MQ=255)

                                                                                                                                             gatgcggcgcggcgCCGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGGACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGggc       <  1:4774286/146‑1 (MQ=255)

                                                                                                                                             gatgcGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGggc       <  1:6511783/146‑1 (MQ=255)

                                                                                                                                              aagcGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGggcg      <  2:3914610/147‑1 (MQ=255)

                                                                                                                                                 cggcgcggcgCCGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCa   <  2:3677890/150‑1 (MQ=255)

                                                                                                                                                  ggCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGACCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCAc  >  2:2502483/1‑150 (MQ=255)

ATCGTGTCGGTGCTGGTAGTCACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTGCGCAGCGCCCCCCATAAAAGCG‑CCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCTGCTGGGCGGCAC  >  NC_002947/4163308‑4163603

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