BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,480,369	A→T	11.9%	K34M (AAG→ATG)	PP_2172 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,480,369	0	A	T	11.9%	342.7 / 3.5	126	K34M (AAG→ATG)	PP_2172	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (86/18); new base T (3/13); total (89/37)
Fisher's exact test for biased strand distribution p-value = 7.89e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 127 total aligned reads displayed.

CAACCCGCACTTTAGCTACCCACTCGCGCCTCCGGAGCCGAACAAAAAATGAAAAAACCTTTGGTTGCTGTCTTTGCTTGTGCATTGGCACTCGGGATGACGGCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAGGGTGGTGGTGGTGGTGGTGGTGGTGGTGGCGGTGGCGGTGGTCATGGCGGTGGAAACGGTGCCGGACATGGCGGTGGCATGGGAGGGCATAGTGGAAGCGGGAAGGGCCTTGGAAGTGATCACGCTGGAAAGGCGACACGCGATCAT  >  NC_002947/2480221‑2480517

cAACCCGCACTTTAGCTACCCACTCGCGCCTCCGGAGCCGAACAAAAAATGAAAAAACCTTTGGTTGCTGTCTTTGCTTGTGCATTGGCACTCGGGATGACGGCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAg                                                                                                                                                     >  2:3265277/1‑150 (MQ=255)

cAACCCGCACTTTAGCTACCCACTCGCGCCTCCGGAGCCGAACAAAAAATGAAAAAACCTTTGGTTGCTGTCTTTGCTTGTGCATTGGCACTCGGGATGACGGCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAg                                                                                                                                                     >  1:3977220/1‑150 (MQ=255)

cAACCCGCACTTTAGCTACCCACTCGCGCCTCCGGAGCCGAACAAAAAATGAAAAAACCTTTGGTTGCTGTCTTTGCTTGTGCATTGGCACTCGGGATGACGGCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAg                                                                                                                                                     >  1:4428010/1‑150 (MQ=255)

 aaCCCGCACTTTAGCTACCCACTCGCGCCTCCGGAGCCGAACAAAAAATGAAAAAACCTTTGGTTGCTGTCTTTGCTTGTGCATTGGCACTCGGGATGACGGCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAGg                                                                                                                                                    <  2:4139372‑M2/150‑2 (MQ=255)

  aCCCGCACTTTAGCTACCCACTCGCGCCTCCGGAGCCGAACAAAAAATGAAAAAACCTTTGGTTGCTGTCTTTGCTTGTGCATTGGCACTCGGGATGACGGCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAGgg                                                                                                                                                   <  1:3748411‑M2/150‑3 (MQ=255)

    ccGCACTTTAGCTACCCACTCGCGCCTCCGGAGCCGAACAAAAAATGAAAAAACCTTTGGTTGCTGTCTTTGCTTGTGCATTGGCACTCGGGATGACGGCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAGggtg                                                                                                                                                 <  1:3271146‑M2/150‑5 (MQ=255)

       cACTTTAGCTACCCACTCGCGCCTCCGGAGCCGAACAAAAAATGAAAAAACCTTTGGTTGCTGTCTTTGCTTGTGCATTGGCACTCGGGATGACGGCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAGggtggtg                                                                                                                                              >  1:2240001‑M2/1‑143 (MQ=255)

          tttAGCTACCCACTCGCGCCTCCGGAGCCGAACAAAAAATGAAAAAACCTTTGGTTGCTGTCTTTGCTTGTGCATTGGCACTCGGGATGACGGCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAGggtggtggtg                                                                                                                                           <  1:4050645‑M2/150‑11 (MQ=255)

           ttAGCTACCCACTCGCGCCTCCGGAGCCGAACAAAAAATGAAAAAACCTTTGGTTGCTGTCTTTGCTTGTGCATTGGCACTCGGGATGACGGCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAGggtggtggtgg                                                                                                                                          >  2:3283166‑M2/1‑139 (MQ=255)

              gCTACCCACTCGCGCCTCCGGAGCCGAACAAAAAATGAAAAAACCTTTGGTTGCTGTCTTTGCTTGTGCATTGGCACTCGGGATGACGGCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAGggtggtggtggtgg                                                                                                                                       <  1:1356962‑M2/150‑15 (MQ=255)

                 aCCCACTCGCGCCTCCGGAGCCGAACAAAAAATGAAAAAACCTTTGGTTGCTGTCTTTGCTTGTGCATTGGCACTCGGGATGACGGCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAGggtggtggtggtggtgg                                                                                                                                    <  1:961208‑M2/150‑18 (MQ=255)

                  cccACTCGCGCCTCCGGAGCCGAACAAAAAATGAAAAAACCTTTGGTTGCTGTCTTTGCTTGTGCATTGGCACTCGGGATGACGGCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAGggtggtggtggtggtggt                                                                                                                                   >  1:5019028‑M2/1‑132 (MQ=255)

                  cccACTCGCGCCTCCGGAGCCGAACAAAAAATGAAAAAACCTTTGGTTGCTGTCTTTGCTTGTGCATTGGCACTCGGGATGACGGCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAGggtggtggtggtggtggt                                                                                                                                   >  2:1525121‑M2/1‑132 (MQ=255)

                       tCGCGCCTCCGGAGCCGAACAAAAAATGAAAAAACCTTTGGTTGCTGTCTTTGCTTGTGCATTGGCACTCGGGATGACGGCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAGggtggtggtggtggtggtggtgg                                                                                                                              <  2:2477470‑M2/150‑24 (MQ=255)

                        cgcgCCTCCGGAGCCGAACAAAAAATGAAAAAACCTTTGGTTGCTGTCTTTGCTTGTGCATTGGCACTCGGGATGACGGCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAGggtggtggtggtggtggtggtggt                                                                                                                             >  2:1187986‑M2/1‑126 (MQ=255)

                         gcgcCTCCGGAGCCGAACAAAAAATGAAAAAACCTTTGGTTGCTGTCTTTGCTTGTGCATTGGCACTCGGGATGACGGCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAGggtggtggtggtg                                                                                                                                        >  1:694324‑M2/1‑125 (MQ=255)

                         gcgcCTCCGGAGCCGAACAAAAAATGAAAAAACCTTTGGTTGCTGTCTTTGCTTGTGCATTGGCACTCGGGATGACGGCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAGggtggtggtggtg                                                                                                                                        <  2:694324‑M2/138‑14 (MQ=255)

                         gcgcCTCCGGAGCCGAACAAAAAATGAAAAAACCTTTGGTTGCTGTCTTTGCTTGTGCATTGGCACTCGGGATGACGGCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAGggtggtggtggtg                                                                                                                                        >  1:694161‑M2/1‑125 (MQ=255)

                         gcgcCTCCGGAGCCGAACAAAAAATGAAAAAACCTTTGGTTGCTGTCTTTGCTTGTGCATTGGCACTCGGGATGACGGCA‑TTC‑TTGCCCGCCGATATCTCACCGATCAGTTCCGCCTATGCCAAGggtggtggtggtggtggtggtggtg                                                                                                                            <  1:3860368‑M2/150‑26 (MQ=255)

                           gcCTCCGGAGCCGAACAAAAAATGAAAAAACCTTTGGTTGCTGTATTTGCTTGTGCATTGGCACTCGGGATGACGGCA‑TTC‑TTGCCCGCCTATATCTCCCCGATCAGTTCCGCCTATGCCAAGggtggtggtggtggtgttggtggtggt                                                                                                                          <  2:2379966/150‑1 (MQ=255)

                            ccTCCGGAGCCGAACAAAAAATGAAAAAACCTTTGGTTGCTGTCTTTGCTTGTGCATTGGCACTCGGGATGACGGCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAGggtggtggtggtggtggtggtggttgtg                                                                                                                         >  1:2872381‑M2/1‑122 (MQ=255)

                             cTCCGGAGCCGAACAAAAAATGAAAAAACCTTTGGTTGCTGTCTTTGCTTGTGCATTGGCACTCGGGATGACGGCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCACGGGtggtggtggtggtggtggtggtggtgg                                                                                                                        <  1:3501880/150‑1 (MQ=255)

                               ccGGAGCCTAACAAAAAATTAAAAAACCTTTGGTTGCTGTCTTTGCTTGTGCATTGGCACTCGGGATGACGGCA‑TTC‑TTGCCCGCCTATATCTCCCCGATCAGTTCCGCCTATGCCAAGGGTGGTGGTGGTGGTGGTGGTGGTggtggcg                                                                                                                      <  2:3081404/150‑1 (MQ=255)

                                 ggACCCGAACAAAAAATCAAAAAACCTTTGGTTGCTGTCTTGGCTTTTGATTTGGCACTGGGGATGACGGAA‑TTC‑TTGCCCGCCTATATCGCCCCTATCCTTACCGCCTAGGCCATGGTTGGTGGTGTTGTTGGTGGTGTTggtggcggt                                                                                                                    <  2:3124353/150‑1 (MQ=255)

                                  gAGCCTAAAAAAAAATTAAAAAACCTTTGTTTGCTTTCTTTGTTTGTGCATTGGCACTCGGGTTGAGGGAT‑TTC‑TTGCCCGCCTATATCTCCCCGATCAGTTCCGCCTTTGCCATGggtggtggtggtggtggtggtggtggcggtggcg                                                                                                                   >  1:2983285‑M2/1‑116 (MQ=255)

                                       gAACAAAAAATGAAAAAACCTTTGGTTGGTGTCTTTGCTTGTGCATTGGCACTCGGGATGACGGCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAGGGTGGTGGTGGTGGTGGTGGTGGTGGTggcggtgggggt                                                                                                              >  2:2408942/1‑150 (MQ=255)

                                       gAACAAAAAATGAAAAAACCTTTGGTTGCTGTCTTTGCTTGTGCATTGGCACTCGGGATGACGGCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAGGGTGGTGGTGGTGGTGGTGGTGGTGGTggcggtggcggt                                                                                                              >  1:1609673/1‑150 (MQ=255)

                                       gAACAAAAAATGAAAAAACCTTTGGTTGCTGTCTTTGCTTGTGCATTGGCACTCGGGATGACGGCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAGGGTGGTGGTGGTGGTGGTGGTGGTGGTggcggtggcggt                                                                                                              >  1:2722731/1‑150 (MQ=255)

                                        aaCAAAAAATCAAAAAACCTTGGGTGGGTGTGTTTGTTGGTG‑ATTGGGCCCGGGGGTTGCGGCATTTG‑TTGCCCCCGGTTACCGCCCCTACCGGTCCGGCCTATGCCAGGGGTGGTGGTGGTGGTGGTGGTGGTGGTGgcggtggcggtg                                                                                                             <  2:313064/150‑1 (MQ=255)

                                         aCAAAAAATGAAAAAACCTTTGGTTGCTGTCTTTGCTTGTGCATTGGCACTCGGGATGACGGCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAGGGTGGTGGTGGTGGTGGTGGTGGTGGTGGcggtggcggtgg                                                                                                            >  2:3760813/1‑150 (MQ=255)

                                          cAAAAAATGAAAAAACCTTTGGTTGCTGTCTTTGCTTGTGCATTGGCACTCGGGATGACGGCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAGGGTGGTGGTGGTGGTGGTGGTGGTGGTGGCGGTGGCggtggt                                                                                                           >  1:5010304/1‑150 (MQ=255)

                                           aaaaaaTGAAAAAACCTTTGGTTGCTGTCTTTGCTTGTGCATTGGCACTCGGGATGACGGCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTGCGCCTATGCGATGGGTGGTGGTGGTGGTGGTGGTGGTGGTGGGGGTGGCGGTGGTc                                                                                                          >  2:2163614/1‑150 (MQ=255)

                                           aaaaaaTGAAAAAACCTTTGGTTGCTGTCTTTGCTTGTGCATTGGCACTCGGGATGACGGCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAGggtggtggtggtggtggtggtggtggcggtggcggtggtcatg                                                                                                          <  1:4406398‑M2/150‑44 (MQ=255)

                                           aaaaaaTGAAAAAACCTTTGGTTGCTGTCTTTGCTTGTGCATTGGCACTCGGGATGACGGCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAGggtggtggtggtggtggtggt                                                                                                                                <  1:3981603‑M2/128‑22 (MQ=255)

                                           aaaaaaTGAAAAAACCTTTGGTTGCTGTCTTTGCTTGTGCATTGGCACTCGGGATGACGGCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAGggtggtggtggtggtggtggt                                                                                                                                >  2:3981603‑M2/1‑107 (MQ=255)

                                           aaaaaaTGAAAAAACCTTTGGTTGCTGTCTTTGCTTGTGCATTGGCACTCGGGATGACGGCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAGGGTGGTGGTGGTGGTGGTGGTGGTGGTGGCGGTGGCGGTGGTc                                                                                                          >  1:1607109/1‑150 (MQ=255)

                                           aaaaaaTGAAAAAACCTTTGGTTGCTGTCTTTGCTTGTGCATTGGCACTCGGGATGACGGCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAGGGTGGTGGTGGTGGTGGTGGTGGTGGTGGCGGTGGCGGTGGTc                                                                                                          >  1:512010/1‑150 (MQ=255)

                                             aaaaTGAAAAAACCTTTGGTTGCTGTCTTTGCTTGTGCATTGGCACTCGGGATGACGGCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAGggtggtggtggtggtggtggtggtggcggtggcggtggtcatggc                                                                                                        <  1:3067407‑M2/150‑46 (MQ=255)

                                             aaaaTGAAAAAACCTTTGGTTGCTGTCTTTGCTTGTGCATTGGCACTCGGGATGACGGCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCGATGCCAAGGGTGGTGGTGGTGGTGGTGGTGGTGGTGGCGGTGGCGGTGGTCAt                                                                                                        >  2:4471768/1‑150 (MQ=255)

                                              aaaTGAAAAAACCTTTGGTTGCTGTCTTTGCTTGTGCATTGGCACTCGGGATGACGGCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAGGGTGGTGGTGGTGGTGGTGGTGGTGGTGGCGGTGGCGGTGGTCATg                                                                                                       >  2:4868980/1‑150 (MQ=255)

                                              aaaTGAAAAAACCTTTGGTTGCTGTCTTTGCTTGTGCATTGGCACTCGGGATGACGACA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAGGGTGGTGGTGGTGGTGGTGGTGGTGGTGGCGGTGGCGGTGGTCATg                                                                                                       >  1:1530650/1‑150 (MQ=255)

                                               aaTGAAAAAACCTTTGGTTGCTGTCTTTGCTTGTGCATTGGCACTCGGGATGACGGCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAGGGTGGTGGTGGTGGTGGTGGTGGTGGTGGCGGTGGCGGTGGTCATgg                                                                                                      >  1:2043618/1‑150 (MQ=255)

                                                    aaaaaCCTTTGGTTGCTGTCTTTGCTTGTGCATTGGCACTCGGGATGACGGCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAGGGTGGTGGTGGTGGTGGTGGTGGTGGTGGCGGTGGCGGTGGTCATGGCGGTg                                                                                                 >  1:3095738/1‑150 (MQ=255)

                                                     aaaaCCTTTGGTTGCTGTCTTTGCTTGTGCATTGGCACTCGGGATGACGGCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAGGGTGGTGGTGGTGGTGGTGGTGGTGGTGGCGGTGGCGGTGGTCATGGCGGTgg                                                                                                >  2:4393915/1‑150 (MQ=255)

                                                           tttGGTTGCTTTCTTGGGTTGTGATTTGGCAGGGGGTTTGCCGGCT‑TTC‑TGGCCCGCCGATATCGCCCCGAGCCGTGCCGCATAGGCCATGGGTGGTGGTGGGGGTGGTGGTGGGGGGGGCGGTGGCGGTGGTCATGGCGGTGGAAACgg                                                                                          <  2:2043618/150‑1 (MQ=255)

                                                           tttGGTTGCTGTCTTTGCTTGTGCATTGGCACTCGGGATGACGGGA‑TTC‑TTGCCCGCCGATATCTCCCCGTTCAGTTCCGCCTATGGCAAGGGTGGTGGTGGTGGTGGTGGTGGTGGTGGCGGTGGCGGTGGTCATGGCGGTGGAAACgg                                                                                          >  2:1559587/1‑150 (MQ=255)

                                                           tttGGTTGCTGTCTTTGCTTGTGCATTGGCACTCGGGATGACGGCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAGGGTGGTGGTGGTGGTGGTGGTGGTGGTGGCGGTGGCGGTGGTCATGGCGGTGGAAACgg                                                                                          >  2:513946/1‑150 (MQ=255)

                                                           tttGGTTGCTGTCTTTGCTTGTGCATTGGCACTCGGGATGACGGCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAGGGTGGTGGTGGTGGTGGTGGTGGTGGTGGCGGTGGCGGTGGTCATGGCGGTGGAAACgg                                                                                          >  2:1318757/1‑150 (MQ=255)

                                                           tttGGTTGCTGTCTTTGCTTGTGCATTGGCACTCGGGATGACGGCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAGGGTGGTGGTGGTGGTGGTGGTGGTGGTGGCGGTGGCGGTGGTCATGGCGGTGGAAACgg                                                                                          >  2:2077932/1‑150 (MQ=255)

                                                                   cTGTCTTTGCTTGTGCATTGGCACTCGGGATGACGGCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAGGGTGGTGGTGGTGGTGGTGGTGGTGGTGGCGGTGGCGGTGGTCATGGCGGTGGAAACGGTGCCGGAc                                                                                  >  2:1413300/1‑150 (MQ=255)

                                                                      tCTTTGCTTGTGCATTGGCACTCGGGATGACGGCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAGGGTGGTGGTGGTGGTGGTGGTGGTGGTGGCGGTGGCGGTGGTCATGGCGGTGGAAACGGTGCCGGACATg                                                                               >  1:1793776/1‑150 (MQ=255)

                                                                        tttGCTTGTGCATTGGCACTCGGGATGACGGCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAGGGTGGTGGTGGTGGTGGTGGTGGTGGTGGTGGTGGCGGTGGTCATGGCGGTGGAAACGGTGCCGGACATGGc                                                                             >  1:4521623/1‑150 (MQ=255)

                                                                           gCTTGTGATTTGGCACTGGGTTTGAGGGCT‑TTG‑TGGCCCGCCTATATCTCCCCGATCAGTCCCGCCTAGGCCACGGGTGGTGTGGGTGGTGGTGGTGGTGGTGGCGGTGGCGGTGGTCATGGCGGTGGAAACGGTGCCGGACATGGCGGt                                                                          <  2:5510/150‑1 (MQ=255)

                                                                               gtgCATTGGCACTCGGGATGACGGCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAGGGTGGTGGTGGTGGTGGTGGTGGTGGTGGCGGTGGCGGTGGTCATGGCGGTGGAAACGGTGCCGGACATGGCGGTGGCa                                                                      >  1:4385554/1‑150 (MQ=255)

                                                                               gtgCATTGGCACTCGGGATGACGGCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAGGGTGGTGGTGGTGGTGGTGGTGGTGGTGGCGGTGGCGGTGGTCATGGCGGTGGAAACGGTGCCGGACATGGCGGTGGCa                                                                      >  1:4385550/1‑150 (MQ=255)

                                                                                 gCATTGGCACTCGGGATGACGGCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAGGGTGGTGGTGGTGGTGGTGGTGGTGGTGGCGGTGGCGGTGGTCATGGCGGTGGAAACGGTGCCGGACATGGCGGTGGCagg                                                                    >  2:3719361/1‑148 (MQ=255)

                                                                                  ctttGGCAGTGGGGTTGACTGCT‑TTC‑TTGCCCGCCGATTTCGCCCCGATCAGTTCCGCCTATGCCACGGGTGGGGGTGGTGGTGGTGGTGGTGGTGGCGGTGGCGGTGGTCATGGCGGTGGAAACGGTGCCGGACATGGCGGTGGCATgg                                                                   <  2:4271843/148‑1 (MQ=255)

                                                                                  cATTGGCACTCGGGATGACGGCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAGGGTGGTGGTGGTGGTGGTGGTGGTGGTGGCGGTGGCGGTGGTCATGGCGGTGGAAACGGTGCCGGACATGGCGGTgcaggga                                                                   >  1:1587159/1‑144 (MQ=255)

                                                                                  cATTGGCACTCGGGATGACGGCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAGGGTGGTGGTGGTGGTGGTGGTGGTGGTGGCGGTGGCGGTGGTCATGGCGGTGGAAACGGTGCCGGACATGGCGGTGGCATgg                                                                   >  2:877047/1‑150 (MQ=255)

                                                                                   aTTGGCACTCGGGATGACGGCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAGggtggtggtggtggtggtggtggtggcggtggcggtggtcatggcggtggaaacggtgccggacatggcggtggcatgggagg                                                                  <  1:5062550‑M2/150‑84 (MQ=255)

                                                                                   aTTGGCACTCGGGATGACGGCA‑TTC‑TTGCCCGCCGATATCTCACCGATCAGTTACGCCTATGCCAAGGGTGGTGGTGGTGGTGGTGGTGGTGGTGGCGGTGGCGGTGGTCATGGCGGTGGAAACGGTGCCGGACATGGCGGTGGCATggg                                                                  >  2:544070/1‑150 (MQ=255)

                                                                                      ggCACTGGGGTTGAGGGAT‑TTC‑TTGCCCGCCTATATTGCCCCGATCTTTTCCGCCTATGCCATGGGTGGTGGTGGTGGTGGTGGTGGTGGTGGCGGTGGCGGTGGTCATGGCGGTGGAAACGGTGCCGGACATGGCGGTGGCATGGGAgg                                                               <  2:622000/150‑1 (MQ=255)

                                                                                         aCTCGGGATGACGGCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAGGGTGGTGGTGGTGGTGGTGGTGGTGGTGGCGGTGGCGGTGGTCATGGCGGTGGAAACGGTGCCGGACATGGCGGTGGCATGGGAGGGCa                                                            >  1:1203698/1‑150 (MQ=255)

                                                                                         aCTCGGGATGACGGCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAGGGTGGTGGTGGTGGTGGTGGTGGTGGTGGCGGTGGCGGTGGTCATGGCGGTGGAAACGGTGCCGGACATGGCGGTGGCATGGGAGGGCa                                                            >  1:3410543/1‑150 (MQ=255)

                                                                                         aCTCGGGATGACGGCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAGGGTGGTGGTGGTGGTGGTGGTGGTGGTGGCGGTGGCGGTGGTCATGGCGGTGGAAACGGTGCCGGACATGGCGGTGGCATGGGAGGGCa                                                            >  2:3437403/1‑150 (MQ=255)

                                                                                           tCGGGATGACGGCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAGGGTGGTGGTGGTGGTGGTGGTGGTGGTGGCGGTGGCGGTGGTCATGGCGGTGGAAACGGTGCCGGACATGGCGGTGGCATGGGAGGGCATa                                                          >  1:328745/1‑150 (MQ=255)

                                                                                             ggtttGACGGAA‑TCT‑TGGCCCGCCGATTTCTCCCCGATCCGTTCCGCATATGCCATGGGTGGTGGTGGTGGTGGTGGTGGTGGTGGCGGTGGCGGTGGTCATGGCGGTGGAAACGGTGCCGGACATGGCGGTGGCATGGGAGGGCATAGt                                                        <  2:4104227/146‑1 (MQ=255)

                                                                                                  gACGGCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAGGGTGGTGGTGGTGGTGGTGGTGGTGGTGGCGGTGGCGGTGGTCATGGCGGTGGAAACGGTGCCGGACATGGCGGTGGCAGGGGAGGGCATAGTGGAAg                                                   >  1:5118665/1‑150 (MQ=255)

                                                                                                     ggCA‑TTCGTTGCCCGCGTATCCCCCCCCGTTCAGTGCGGCATTGCCCTGGGGTGGGGGGGGGGGTGGGGGGGGGGGGGGGGGTGGCGGTGGTCATGGCGGTGGAAACGGTGCCGGACATGGCGGTGGCATGGGAGGGCATAGTGGAAGCg                                                 <  1:592218/150‑1 (MQ=255)

                                                                                                     ggCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAGGGTGGTGGTGGTGGTGGTGGTGGTGGTGGCGGTGGCGGTGGTCCTGGCGGTGGACACGGTGCCGGACATGGCGGTGGCAGGGGAGGGCATAGTGGAAGCgg                                                >  2:1564372/1‑150 (MQ=255)

                                                                                                     ggCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAGGGTGGTGGTGGTGGTGGTGGTGGTGGTGGCGGTGGCGGTGGTCATGGCGGTGGAAACGGTGCCGGACATGGCGGTGGCATGGGAGGGCATAGTGGAAGCgg                                                >  1:2135872/1‑150 (MQ=255)

                                                                                                     ggCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAGGGTGGTGGTGGTGGTGGTGGTGGTGGTGGCGGTGGCGGTGGTCATGGCGGTGGAAACGGTGCCGGACATGGCGGTGGCATGGGAGGGCATAGTGGAAGCgg                                                >  2:675927/1‑150 (MQ=255)

                                                                                                     ggCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAGGGTGGTGGTGGTGGTGGTGGTGGTGGTGGCGGTGGCGGTGGTCATGGCGGTGGAAACGGTGCCGGACATGGCGGTGGCATGGGAGGGCATAGTGGAAGCgg                                                >  2:3969742/1‑150 (MQ=255)

                                                                                                     ggCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAGGGTGGTGGTGGTGGTGGTGGTGGTGGTGGCGGTGGCGGTGGTCATGGCGGTGGAAACGGTGCCGGACATGGCGGTGGCATGGGAGGGCATAGTGGAAGCgg                                                >  2:3969741/1‑150 (MQ=255)

                                                                                                       cA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAGggtggtggtggtggtggtggtggtggcggtggcggtggtcatggcggtggaaacggtgccggacatggcggtggcatgggagggcatagtggaagcgggaagg                                              <  1:1686394‑M2/150‑104 (MQ=255)

                                                                                                                gCCCGCCGATATCCCCCCGACCAGTGCCGCCATGGCCATGgggggtggtggtggtggtggtgggggcggtggcggtggtcatggcggtggaaacggtgccggacatggcggtggcatgggagggcatagtggaagcgggaagggccttgg                                       >  2:5062550‑M2/1‑40 (MQ=255)

                                                                                                                 cccGCCGATATCTCCCCTACCCGTCCCCCCTAGGCCATGGGTGGTGGGGGTGGGGGTGGTGGTGGTGGGGGTGGCGGTGGTCATGGCGGTGGAAACGGTGCCGGACATGGCGGTGGCATGGGAGGGCATAGTGGAAGCGGGAAGGGCCtt                                      <  2:4428010/150‑1 (MQ=255)

                                                                                                                  ccgcctATCTCTCCCCTACCAGGTCCGCCTAGGCCACGGGTGGTGGTGGTGGGGGTGGTGGTGGTGGCGGTGGCGGTGGTCATGGCGGTGGAAACGGTGCCGGACATGGCGGTGGCATGGGAGGGCATAGTGGAAGCGGGAAGGGCCTTg                                     <  2:3275063/150‑1 (MQ=255)

                                                                                                                   cgccgATATCTCCCCGATCAGTTCCGCCTATGCCAAGGGTGGTGGTGGTGGTGGTGGTGGTGGTGGCGGTGGCGGTGGTCATGGCGGTGGAAACGGTGCCGGACATGGCGGTGGCATGGGAGGGCATAGTGGAAGCGGGAAGGGCCTTgg                                    >  1:1373809/1‑150 (MQ=255)

                                                                                                                           tctcCCCGTTCTGTTCCGCCTTTGCCTTGGGTGGTGGTGGTGGTGGTGGTGGTGGTGGCGGTGGCGGTGGTCATGGCGGTGGAAACGGTGCCGGACATGGCGGTGGCATGGGAGGGCATAGTGGAAGCGGGAAGGGCCTTGGAAGTGATc                            <  2:1793776/150‑1 (MQ=255)

                                                                                                                           tctcCCCGATCAGTTCCGCCTATGCCAAGGGTGGTGGTGGTGGTGGTGGTGGTGGTGGCGGTGGCGGTGGTCATGGCGGTGGAAACGGTGCCGGACATGGCGGTGGCATGGGAGGGCATAGTGGAAGCGGGAAGGGCCTTGGAAGTGATc                            >  1:518014/1‑150 (MQ=255)

                                                                                                                             tcCCCGCTCAGTTCCGCCTATGCCAAGGGTGGTGGTGGTGGTGGTGGTGGTGGTGGCGGTGGCGGTGGTCATGGCGGTGGAAACGGTGCCGGACATGGCGGTGGCATGGGAGGGCATAGTGGAAGCGGGAAGGGCCTTGGAAGTGATCAc                          >  1:3961485/1‑150 (MQ=255)

                                                                                                                             nccccGATCAGTTCCGCCTATGCCAAGGGTGGTGGTGGTGGTGGTGGTGGTGGTGGCGGTGGCGGTGGTCATGGCGGTGGAAACGGTGCCGGACATGGCGGTGGCATGGGAGGGCATAGTGGAAGCGGGAAGGGCCTTGGAAGTGATCAc                          >  1:3771054/2‑150 (MQ=255)

                                                                                                                              ccccGATCAGTTCCGCCTATGCCAAGGGTGGTGGTGGTGGTGGTGGTGGTGGTGGCGGTGGCGGTGGTCATGGAGGTGGAAACGGTGCCGGACATGGCGGTGGCATGGGAGGGCATAGTGGAAGCGGGAAGGGCCTTGGAAGTGATCACg                         >  1:4569220/1‑150 (MQ=255)

                                                                                                                               cccGTTCTGTCCCGCCTATGCCATGGGGGGTGGGGGTGGTGGGGGTGGTGGTGGCGGTGGCGGTGGTCATGGCGGTGGAAACGGTGCCGGACATGGCGGTGGCATGGGAGGGCATAGTGGAAGCGGGAAGGGCCTTGGAAGTGATCACGc                        <  2:5301574/150‑1 (MQ=255)

                                                                                                                               cccGATCAGTTCCGCCTATGCCAAGGGTGGTGGTGGTGGTGGTGGTGGTGGTGGCGGTGGCGGTGGTCATGGCGGTGGAAACGGTGCCGGACATGGCGGTGGCATGGGAGGGCATAGTGGAAGCGGGAAGGGCCTTGGAAGTGATCACGc                        >  1:4231621/1‑150 (MQ=255)

                                                                                                                               cccGATCAGTTCCGCCTATGCCAAGGGTGGTGGTGGTGGTGGTGGTGGTGGTGGCGGTGGCGGTGGTCATGGCGGTGGAAACGGTGCCGGACATGGCGGTGGCATGGGAGGGCATAGTGGAAGCGGGAAGGGCCTTGGAAGTGATCACGc                        >  1:4266541/1‑150 (MQ=255)

                                                                                                                                 cGATCAGTTCCGCCTATGCCAAGGGTGGTGGTGGTGGTGGTGGTGGTGGTGGCGGTGGCGGTGGTCATGGCGGTGGAAACGGTGCCGGACATGGCGGTGGCATGGGAGGGCATAGTGGAAGCGGGAAGGGCCTTGGAAGTGATCACGCTg                      >  1:3383271/1‑150 (MQ=255)

                                                                                                                                 cGATCAGTTCCGCCTATGCCAAGGGTGGTGGTGGTGGTGGTGGTGGTGGTGGCGGTGGCGGTGGTCATGGCGGTGGAAACGGTGCCGGACATGGCGGTGGCATGGGAGGGCATAGTGGAAGCGGGAAGGGCCTTGGAAGTGATCACGCTg                      >  2:3386060/1‑150 (MQ=255)

                                                                                                                                  gATCAGTTCCGCCTATGCCAAGGGTGGTGGTGGTGGTGGTGGTGGTGGTGGCGGTGGCGGTGGTCATGGCGGTGGAAACGGTGCCGGACATGGCGGTGGCATGGGAGGGCATAGTGGAAGCGGGAAGGGCCTTGGAAGTGATCACGCTgg                     >  1:4609300/1‑150 (MQ=255)

                                                                                                                                     cAGTTCGGCCTAGGCCATGGGGGGTGGGGGTGGTGGTGGTGGGGGGGGCGGTGGCGGTGGTCATGGCGGTGGAAACGGTGCCGGACATGGCGGTGGCATGGGAGGGCATAGTGGAAGCGGGAAGGGCCTTGGAAGTGATCACGCTGGaaa                  <  2:3748411/150‑1 (MQ=255)

                                                                                                                                       gTTCCGCCTATGCCAAGGGTGGTGGTGGTGGTGGTGGTGGTGGTGGCGGTGGCGGTGGTCATGGCGGTGGAAACGGTGCCGGACATGGCGGTGGCATGGGAGGGCATAGTGGAAGCGGGAAGGGCCTTGGAAGTGATCACGCTGGAAAgg                >  1:2881738/1‑150 (MQ=255)

                                                                                                                                        ttCCGCCTATGCCAAGGGTGGTGGTGGTGGTGGTGGTGGTGGTGGCGGTGGCGGTGGTCATGGCGGTGGAAACGGTGCCGGACATGGCGGTGGCATGGGAGGGCATAGTGGAAGCGGGAAGGGCCTTGGAAGTGATCACGCTGGAAAGGc               >  2:3964015/1‑150 (MQ=255)

                                                                                                                                             ccTATGCCAAGGGTGGTGGTGGTGGTGGTGGTGGTGGTGGCGGTGGCGGTGGTCATGGCGGTGGAAACGGTGCCGGACATGGCGGTGGCATGGGAGGGCATAGTGGAAGCGGGAAGGGCCTTGGAAGTGATCACGCTGGAAAGGCGacac          >  2:410015/1‑150 (MQ=255)

                                                                                                                                               tATGCCAAGGGTGGTGGTGGTGGTGGTGGTGGTGGTGGCGGTGGCGGTGGTCATGGCGGTGGAAACGGTGCCGGACATGGCGGTGGCATGGGAGGGCATAGTGGAAGCGGGAAGGGCCTTGGAAGTGATCACGCTGGAAAGGCGACAcgc        >  1:3068322/1‑150 (MQ=255)

                                                                                                                                                 tGCCATGGGGGGTGGTGGTGGTGTTGGTGGTGGTGGCGGTGGCGGTGGTCATGGCGGTGGAAACGGTGCCGGACATGGCGGTGGCATGGGAGGGCATAGTGGAAGCGGGAAGGGCCTTGGAAGTGATCACGCTGGAAAGGCGACACGCGa      <  2:961208/150‑1 (MQ=255)

                                                                                                                                                   ccAAGGGTGGTGGTGGTGGTGGTGGTGGTGGTGGCGGTGGCGGTGGTCATGGCGGTGGAAACGGTGCCGGACATGGCGGTGGCATGGGAGGGCATAGTGGAAGCGGGAAGGGCCTTGGAAGTGATCACGCTGGAAAGGCGACACGCGATc    >  1:1610844/1‑150 (MQ=255)

                                                                                                                                                    cAAGGGTGGTGGTGGTGGTGGTGGTGGTGGTGGTGGTGGCGGTGGTCATGGCGGTGGAAACGGTGCCGGACATGGCGGTGGCATGGGAGGGCATAGTGGAAGCGGGAAGGGCCTTGGAAGTGATCACGCTGGAAAGGCGACACGCGATCa   >  1:2318552/1‑150 (MQ=255)

                                                                                                                                                    cAAGGGTGGTGGTGGTGGTGGTGGTGGTGGTGGCGGTGGCGGTGGTCATGGCGGTGGAAACGGTGCCGGACATGGCGGTGGCATGGGAGGGCATAGTGGAAGCGGGAAGGGCCTTGGAAGTGATCACGCTGGAAAGGCGACACGCGATCa   >  2:2831738/1‑150 (MQ=255)

                                                                                                                                                     aaGGGTGGTGGTGGTGGTGGTGGTGGTGGTGGCGGTGGCGGTGGTCATGGCGGTGGAAACGGTGCCGGACATGGCGGTGGCATGGGAGGGCATAGTGGAAGCGGGAAGGGCCTTGGAAGTGATCACGCTGGAAAGGCGACACGCGATCAt  >  1:1853833/1‑150 (MQ=255)

CAACCCGCACTTTAGCTACCCACTCGCGCCTCCGGAGCCGAACAAAAAATGAAAAAACCTTTGGTTGCTGTCTTTGCTTGTGCATTGGCACTCGGGATGACGGCA‑TTC‑TTGCCCGCCGATATCTCCCCGATCAGTTCCGCCTATGCCAAGGGTGGTGGTGGTGGTGGTGGTGGTGGTGGCGGTGGCGGTGGTCATGGCGGTGGAAACGGTGCCGGACATGGCGGTGGCATGGGAGGGCATAGTGGAAGCGGGAAGGGCCTTGGAAGTGATCACGCTGGAAAGGCGACACGCGATCAT  >  NC_002947/2480221‑2480517

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