BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,294,586	A→T	5.3%	L246H (CTC→CAC)	PP_4663 ←	transmembrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,294,586	0	A	T	5.3%	848.4 / 10.1	288	L246H (CTC→CAC)	PP_4663	transmembrane protein
Reads supporting (aligned to +/- strand): ref base A (161/94); new base T (3/17); total (164/124)
Fisher's exact test for biased strand distribution p-value = 4.89e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 289 total aligned reads displayed.

AGGCCGCGCTGGACGTGCACTCGGCCCTTCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCACCCAAGGCATCCTGCCCG  >  NC_002947/5294438‑5294735

aGGCCGCGCTGGACGTGCACTCGGCCCTTCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAg                                                                                                                                                      <  1:5327319/150‑1 (MQ=255)

 ggCCGCGCTGGACGTGCACTCGGCCCTTCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAgg                                                                                                                                                     >  1:1464073/1‑150 (MQ=255)

  gCCGCGCTGGACGTGCACTCGGCCCTTCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAggg                                                                                                                                                    >  2:4153389/1‑150 (MQ=255)

    cgcgCTGGACGTGCACTCGGCCCTTCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGggcg                                                                                                                                                  >  1:2647619/1‑150 (MQ=255)

    cgcgCTGGACGTGCACTCGGCCCTTCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGggcg                                                                                                                                                  <  2:5406121/150‑1 (MQ=255)

    cgcgCTGGACGTGCACTCGGCCCTTCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGggcg                                                                                                                                                  <  1:1746149/150‑1 (MQ=255)

    cgcgCTGGACGTGCACTCGGCCCTTCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGggcg                                                                                                                                                  >  2:153126/1‑150 (MQ=255)

    cgcgCTGGACGTGCACTCGGCCCTTCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGggcg                                                                                                                                                  >  2:2286112/1‑150 (MQ=255)

    cgcgCTGGACGTGCACTCGGCCCTTCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGggcg                                                                                                                                                  <  2:3154322/150‑1 (MQ=255)

        cTGGACGTGCACTCGGCCCTTCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCa                                                                                                                                              <  1:4964086/150‑1 (MQ=255)

         tGGACGTGCACTCGGCCCTTCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAg                                                                                                                                             <  1:1545947/150‑1 (MQ=255)

         tGGACGTGCACTCGGCCCTTCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCCTCGGCGGCCac                                                                                                                                             <  1:3646374/150‑2 (MQ=255)

          ggACGTGCACTCGGCCCTTCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCacc                                                                                                                                            <  1:2274590/150‑3 (MQ=255)

           gACGTGCACTCGGCCCTTCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACGGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCacct                                                                                                                                           <  2:58121/150‑4 (MQ=255)

           gACGTGCACTCGGCCCTTCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGc                                                                                                                                           <  1:2693004/150‑1 (MQ=255)

                gCACTCGGCCCTTCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATg                                                                                                                                      >  2:5277550/1‑150 (MQ=255)

                  aCTCGGCCCTTCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGaa                                                                                                                                    >  2:3503049/1‑150 (MQ=255)

                    tCGGCCCTTCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGaaaa                                                                                                                                  >  1:1020134/1‑150 (MQ=255)

                    tCGGCCCTTCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGaaaa                                                                                                                                  >  1:1124808/1‑150 (MQ=255)

                    tCGGCCCTTCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGaaaa                                                                                                                                  >  1:4414520/1‑150 (MQ=255)

                    tCGGCCCTTCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGaaaa                                                                                                                                  >  1:3995895/1‑150 (MQ=255)

                     cGGCCCTTCCGCACATGGCCGCCTTCTTCTAAAAGCGCATGGACTACCTGGAGCGCTATTTCGATATCCAGCGGGTCGACATACTCAGCGCTTCGTCCGATAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGaaaaa                                                                                                                                 <  1:3261274/150‑1 (MQ=255)

                     cGGCCCTTCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGaaaaa                                                                                                                                 <  2:374114/150‑1 (MQ=255)

                     cGGCCCTTCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGaaaaa                                                                                                                                 <  2:2117476/150‑1 (MQ=255)

                         ccTTCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCaa                                                                                                                             >  2:850075/1‑150 (MQ=255)

                         ccTTCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCaa                                                                                                                             >  1:2061765/1‑150 (MQ=255)

                           ttCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACg                                                                                                                           >  2:1582221/1‑150 (MQ=255)

                            tCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGc                                                                                                                          <  2:4804447/150‑1 (MQ=255)

                               gCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTg                                                                                                                       >  1:5415708/1‑150 (MQ=255)

                                 acaTGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGaa                                                                                                                     >  2:4549595/1‑150 (MQ=255)

                                 acaTGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGaa                                                                                                                     >  1:2978680/1‑150 (MQ=255)

                                      gccgccTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAAcagca                                                                                                                >  1:2084636/1‑150 (MQ=255)

                                      gccgccTTCTTCGAAAAGAGAAGGTACTACCTGTAGCGGTATTTCTATATCCAGCTGGTCGACATGCTCAGGGATTCGTCCGCTAACCTTTGATCAGCTGGTGGCCGCCGTGGGGTGCCAGGCCGATGAAAAAGCAACGCTTGAAcagca                                                                                                                <  2:2620746/150‑1 (MQ=255)

                                           cttcttCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAgcgcg                                                                                                           >  1:965268/1‑150 (MQ=255)

                                           cttcttCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAgcgcg                                                                                                           >  2:4648930/1‑150 (MQ=255)

                                            ttcttcGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAgcgcgc                                                                                                          >  1:2453880/1‑150 (MQ=255)

                                            ttcttcGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAgcgcgc                                                                                                          >  1:1763356/1‑150 (MQ=255)

                                            ttcttcGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGCCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAgcgcgc                                                                                                          >  1:2918856/1‑150 (MQ=255)

                                               ttcGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGa                                                                                                       >  1:1834914/1‑150 (MQ=255)

                                               ttcGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGa                                                                                                       >  2:1811926/1‑150 (MQ=255)

                                                 cGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATc                                                                                                     >  1:1926576/1‑150 (MQ=255)

                                                  gAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAg                                                                                                                                >  1:162258/1‑122 (MQ=255)

                                                   aaaaGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCg                                                                                                   >  1:839648/1‑150 (MQ=255)

                                                    aaaGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGc                                                                                                  >  2:4056392/1‑150 (MQ=255)

                                                     atgcgcATGGACTACCTGGAGCGCTATTTCGATATCCCGCGGGTCGACATGCTCCGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGTGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCt                                                                                                 <  1:1293527/148‑1 (MQ=255)

                                                           aTGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCa                                                                                           >  2:4909874/1‑150 (MQ=255)

                                                                 tACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAgg                                                                                                                                                     >  1:300456/1‑86 (MQ=255)

                                                                   ccTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTgg                                                                                   >  1:4245954/1‑150 (MQ=255)

                                                                   ccTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTgg                                                                                   >  2:5582527/1‑150 (MQ=255)

                                                                   ccTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTgg                                                                                   >  2:654857/1‑150 (MQ=255)

                                                                   ccTGGAGCGCGAGTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGAGCCGCTGACGCAGCCCCTgg                                                                                   >  2:1028681/1‑150 (MQ=255)

                                                                              aTTTCTATACCCAGCGGGACTAAATCCTCAGCGATACGTCCGCTAACCTTTCAACAGCCGGGGCCCGCCGCGGGCGCCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCAt                                                                        <  1:690223/150‑1 (MQ=255)

                                                                              aTTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCAt                                                                        >  2:1236611/1‑150 (MQ=255)

                                                                               tttCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATg                                                                       >  2:563238/1‑150 (MQ=255)

                                                                               tttCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGAAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCACGCCAGCaag                                                                       >  2:5196629/1‑148 (MQ=255)

                                                                                 tCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGcc                                                                     >  2:3314977/1‑150 (MQ=255)

                                                                                  cGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCg                                                                    >  1:3780833/1‑150 (MQ=255)

                                                                                   gATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTACCCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGa                                                                   <  2:3932712/150‑1 (MQ=255)

                                                                                   gATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGa                                                                   >  2:2925753/1‑150 (MQ=255)

                                                                                    atatCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAATGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAg                                                                  >  2:4318742/1‑150 (MQ=255)

                                                                                    atatCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAg                                                                  <  2:1442977/150‑1 (MQ=255)

                                                                                       tCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCaa                                                               >  1:3913231/1‑150 (MQ=255)

                                                                                         cAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAgg                                                             >  2:1347410/1‑150 (MQ=255)

                                                                                         cAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAgg                                                             >  1:4303717/1‑150 (MQ=255)

                                                                                            cGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTTATCAGCTGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCg                                                          <  2:827311/150‑1 (MQ=255)

                                                                                              ggTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGgc                                                        >  2:1926033/1‑150 (MQ=255)

                                                                                              ggTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGgc                                                        <  2:2986261/150‑1 (MQ=255)

                                                                                                tCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCCGCTGGGGGCCGCCGCGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGgcgc                                                      <  1:3383743/150‑1 (MQ=255)

                                                                                                  gACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCCGCAGGTGGCCGCCGCGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCa                                                    <  1:879219/150‑1 (MQ=255)

                                                                                                  gACATGCTCAGCGCTTCGTCCGATAACCGTTGATCAGCTGGTGGCCGCCGTGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCa                                                    <  1:3512719/150‑1 (MQ=255)

                                                                                                   aCATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAg                                                   >  2:5275732/1‑150 (MQ=255)

                                                                                                       gCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATc                                               >  1:1513097/1‑150 (MQ=255)

                                                                                                        cTCCGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATca                                              <  2:1653692/150‑1 (MQ=255)

                                                                                                           aGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATcagca                                           >  1:2364793/1‑150 (MQ=255)

                                                                                                            gcgcTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATcagcag                                          >  1:360794/1‑150 (MQ=255)

                                                                                                            gcgcTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATcagcag                                          >  1:3407738/1‑150 (MQ=255)

                                                                                                                   gTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGcc                                   <  2:5034798/150‑1 (MQ=255)

                                                                                                                     ccGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGcccc                                 <  2:1464073/150‑1 (MQ=255)

                                                                                                                       gCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCa                               >  1:4145239/1‑150 (MQ=255)

                                                                                                                       gCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCAAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCa                               >  2:2802663/1‑150 (MQ=255)

                                                                                                                        cTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAg                              >  2:3007851/1‑150 (MQ=255)

                                                                                                                            ccGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGcc                          >  2:1745250/1‑150 (MQ=255)

                                                                                                                             cGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCt                         >  1:2011040/1‑150 (MQ=255)

                                                                                                                              gTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTg                        >  2:2505057/1‑150 (MQ=255)

                                                                                                                              gTTGATCAGCAGGAGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTAGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCAAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTg                        >  2:1914137/1‑150 (MQ=255)

                                                                                                                               ttGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTgg                       <  1:5130018/150‑1 (MQ=255)

                                                                                                                                     agcagGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCAccc                 >  2:1170299/1‑150 (MQ=255)

                                                                                                                                     agcagGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCAccc                 >  1:4903128/1‑150 (MQ=255)

                                                                                                                                      gcagGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCACCCa                >  1:1885849/1‑150 (MQ=255)

                                                                                                                                      gcagGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCACCCa                >  1:4734722/1‑150 (MQ=255)

                                                                                                                                       ggcctGGCCGCCCTCGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGACCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCACCCaa               >  2:2226389/5‑150 (MQ=255)

                                                                                                                                       cagGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCCCCCaa               >  2:1069271/1‑150 (MQ=255)

                                                                                                                                       cagGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCACCCaa               >  2:4357921/1‑150 (MQ=255)

                                                                                                                                         ggTGGCCGCCGAGGTCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCACCCAAgg             <  2:5529034/150‑1 (MQ=255)

                                                                                                                                         ggTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCACCCAAgg             <  2:2708099/150‑1 (MQ=255)

                                                                                                                                         cctGGCCGCCCTCGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTgg                       >  2:4272249/3‑140 (MQ=255)

                                                                                                                                              ccgccgAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCACCCAAGGCATcc        >  2:3345359/1‑150 (MQ=255)

                                                                                                                                               cgccgAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCACCCAAGGCATCCt       <  1:4647777/150‑1 (MQ=255)

                                                                                                                                               cgccgAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCACCCAAGGCATCCt       >  2:3136904/1‑150 (MQ=255)

                                                                                                                                                    aGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGAGCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCACCCAAGGCATCCTGCCCg  >  2:2034077/1‑150 (MQ=255)

AGGCCGCGCTGGACGTGCACTCGGCCCTTCCGCACATGGCCGCCTTCTTCGAAAAGCGCATGGACTACCTGGAGCGCGATTTCGATATCCAGCGGGTCGACATGCTCAGCGCTTCGTCCGCTAACCGTTGATCAGCAGGTGGCCGCCGAGGGCGGCCAGGCCGATGAAAAAGCAACGCTTGAACAGCAGCGCGCTGATCCGCTGACGCAGCCACTGGCCCGCCAGCATGCCGAGCAAGGCCGGCGCCAGCATCAGCAGCGACGCCCCCAGCGCCTGGCCACCCAAGGCATCCTGCCCG  >  NC_002947/5294438‑5294735

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