BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	450,290	G→A	100%	intergenic (+45/‑299)	bolA → / → tig	stationary‑phase morphogene, transcriptional repressor for mreB; also regulator for dacA, dacC, and ampC/peptidyl‑prolyl cis/trans isomerase (trigger factor)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	450,290	0	G	A	99.3%	1004.0 / ‑2.6	300	intergenic (+45/‑299)	bolA/tig	stationary‑phase morphogene, transcriptional repressor for mreB; also regulator for dacA, dacC, and ampC/peptidyl‑prolyl cis/trans isomerase (trigger factor)
Reads supporting (aligned to +/- strand): ref base G (2/0); new base A (171/126); total (173/127)
Fisher's exact test for biased strand distribution p-value = 5.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic base substitution creates a homopolymer stretch.

Only 100 of 302 total aligned reads displayed.

CATGCGCTGGCTCTGCATACTTACACTATTAAGGAGTGGGAAGGGTTGCAGGACACCGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTG‑TTTTGT‑CTGG‑ATTTTGCGCTTTTTGCCCAGCATTCAGA‑CG‑AAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAACCGTTCTCGACTCATAAAAGTGATGCC  >  CP009273/450144‑450435

cATGCGCTGGCTCTGCATACTTACACTATTAAGGAGTGGGAAGGGTTGCAGGACACCGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAAt                                                                                                                                                      >  2:2381395/1‑149 (MQ=255)

 aTGCGCTGGCTCTGCATACTTACACTATTAAGGAGTGGGAAGGGTTGCAGGACACCGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAAtt                                                                                                                                                     >  2:128418/1‑149 (MQ=255)

        ggCTCTGCATACTTACACTATTAAGGAGTGGGAAGGGTTGCAGGACACCGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGaa                                                                                                                                              <  1:4418613/149‑1 (MQ=255)

           tctGCATACTTACACTATTAAGGAGTGGGAAGGGTTGCAGGACACCGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGaaaaa                                                                                                                                           <  1:1358895/149‑1 (MQ=255)

             tGCATACTTACACTATTAAGGAGTGGGAAGGGTTGCAGGACACCGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACg                                                                                                                                         >  2:1236450/1‑149 (MQ=255)

             tGCATACTTACACTATTAAGGAGTGGGAAGGGTTGCAGGACACCGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACg                                                                                                                                         >  2:2809499/1‑149 (MQ=255)

               cATACTTACACTATTAAGGAGTGGGAAGGGTTGCAGGACACCGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGc                                                                                                                                       <  2:296314/149‑1 (MQ=255)

                aTACTTACACTATTAAGGAGTGGGAAGGGTTGCAGGACACCGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGcc                                                                                                                                      >  2:3072148/1‑149 (MQ=255)

                 tACTTACACTATTAAGTAGTGGGAAGGGTTGCAGGACACCGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCATTATGTTGCTAAAAATTTTATGAAAAACGGCCt                                                                                                                                     <  2:1955554/149‑1 (MQ=255)

                 tACTTACACTATTAAGGAGTGGGAAGGGTTGCAGGACACCGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCt                                                                                                                                     >  2:3245587/1‑149 (MQ=255)

                 tACTTACACTATTAAGGAGTGGGAAGGGTTGCAGGACACCGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCt                                                                                                                                     >  2:3673103/1‑149 (MQ=255)

                  aCTTACACTATTAAGGAGTGGGAAGGGTTGCAGGACACCGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTg                                                                                                                                    >  2:4280096/1‑149 (MQ=255)

                    ttACACTATTAAGGAGTGGGAAGGGTTGCAGGACACCGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCg                                                                                                                                  >  1:3765341/1‑149 (MQ=255)

                       cacTATTAAGGAGTGGGAAGGGTTGCAGGACACCGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGc                                                                                                                               <  2:1000917/149‑1 (MQ=255)

                       cacTATTAAGGAGTGGGAAGGGTTGCAGGACACCGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTAGGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGc                                                                                                                               >  2:4115394/1‑149 (MQ=255)

                         cTATTAAGGAGTGGGAAGGGTTGCAGGACACCGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCg                                                                                                                             >  1:3059539/1‑149 (MQ=255)

                           aTTAAGGAGTGGGAAGGGTTGCAGGACACCGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCgtt                                                                                                                           <  1:2460106/149‑1 (MQ=255)

                            ttAAGGAGTGGGAAGGGTTGCAGGACACCGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCgttt                                                                                                                          >  2:4555130/1‑149 (MQ=255)

                             tAAGGAGTGGGAAGGGTTGCAGGACACCGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACACATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCgggtt                                                                                                                         >  1:505909/1‑145 (MQ=255)

                              aaGGAGTGGGAAGGGTTGCAGGACACCGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCgtggtt                                                                                                                        >  1:574411/1‑145 (MQ=255)

                                 gAGTGGGAAGGGTTGCAGGACACCGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCgttttg‑ttt                                                                                                                    <  1:3449037/149‑1 (MQ=255)

                                   gTGGGAAGGGTTGCAGGACACCGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGttttg‑ttttg                                                                                                                  >  2:3191570/1‑149 (MQ=255)

                                   gTGGGAAGGGTTGCAGGACACCGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGttttg‑ttttg                                                                                                                  >  2:2292317/1‑149 (MQ=255)

                                   gTGGGAAGGGTTGCAGGACACCGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGtgttg‑ttttg                                                                                                                  >  1:296271/1‑149 (MQ=255)

                                     gggAAGGGTTGCAGGACACCGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCgttttg‑ttttttg                                                                                                                >  1:2249649/1‑146 (MQ=255)

                                     gggAAGGGTTGCAGGACACCGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTTTTG‑TTTTGT‑c                                                                                                               <  2:4749524/149‑1 (MQ=255)

                                             ttGCAGGACACCGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTTTTG‑TTTTGT‑CTGG‑Atttt                                                                                                      <  1:182277/149‑1 (MQ=255)

                                             ttGCAGGACACCGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTGTTT‑TTTTTT‑GTGG‑Ttttt                                                                                                      >  1:3764201/1‑149 (MQ=255)

                                                 aGGACACCGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTGgtg‑gttttt‑gtgg‑tttttttgg                                                                                                  >  1:1224554/1‑145 (MQ=255)

                                                 aGGACACCGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTGTTG‑TTTTGT‑GTGG‑Atttttggc                                                                                                  >  1:4691712/1‑145 (MQ=255)

                                                  ggACACCGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTTTTG‑TTTTGT‑CTGG‑ATTTTGCGCt                                                                                                 >  2:636453/1‑149 (MQ=255)

                                                  ggACACCGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTTTTG‑TTTTGT‑CTGG‑ATTTTGCGCt                                                                                                 <  2:1052359/149‑1 (MQ=255)

                                                  ggACACCGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTTTTG‑TTTTGT‑CTGG‑ATTTTGCGCt                                                                                                 <  1:2194668/149‑1 (MQ=255)

                                                    acacCGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTGTTG‑TTTTGT‑CTGG‑ATTTTTCGCttt                                                                                               >  2:4336231/1‑149 (MQ=255)

                                                     cacCGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTTTTG‑TTTTGT‑CTGG‑ATTTTGCGCtttt                                                                                              <  2:2162125/149‑1 (MQ=255)

                                                      acCGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTGTTG‑TTTTGT‑CTGG‑ATTTTGCGCttttt                                                                                             >  2:1905955/1‑149 (MQ=255)

                                                      acCGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTGTTG‑TTTTGT‑CTGG‑ATTTTGAGCttttt                                                                                             >  2:3614457/1‑149 (MQ=255)

                                                        cGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTGTTG‑TTTTTT‑CTGG‑ATTTTGCGCTTTTTgg                                                                                           >  2:45826/1‑148 (MQ=255)

                                                          tCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTTTTG‑TTTTGT‑CTGG‑ATTTTGCGCTTTTTGccc                                                                                         <  1:2089979/149‑1 (MQ=255)

                                                          tCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTGTTG‑TTTTTTGGTGT‑TTTTTGGGTTTTTTGcc                                                                                          >  1:1678037/1‑149 (MQ=255)

                                                          tCTTTGCCTCTCCTCCCTGTCGTGGAGAATTAAGCAGCGCGTAAAAACGCATTTGAAACTGTCGGCTCTTTTCCAGTATGTTTCTAAATATTTTATGAAAAACGGCCTGCGGGCCGTTTTG‑TTTTGT‑CTGG‑ATTTTGCGCTTTTTGccc                                                                                         <  1:4097297/149‑1 (MQ=255)

                                                            tttGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGGATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTGTTG‑TTTTGG‑CTGG‑AGTTTGCGCTTTTTTccccg                                                                                       >  1:1059186/1‑147 (MQ=255)

                                                               gCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTGTTG‑TTTTGT‑CTGG‑GTTTTGCGCTTTTTGCCCAGCaa                                                                                    >  2:1439466/1‑148 (MQ=255)

                                                               gCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTGTTG‑TTTTGT‑CTGG‑ATTTTTCGCTTTTTGCCCAGCAt                                                                                    >  2:1530144/1‑149 (MQ=255)

                                                               gCCTCTCCGCCCTGGCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTGGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTTTTG‑TTTTGT‑CTGG‑ATTTTGCGCTTTTTGCCCAGCAt                                                                                    <  1:4001559/149‑1 (MQ=255)

                                                                  tctcCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTTTTG‑TTTTGT‑CTTG‑ATTTTGCGCTTTTTGCCCCGCATTCa                                                                                 >  2:1051124/1‑149 (MQ=255)

                                                                    tcctccCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTTTTG‑TTTTGT‑CTGG‑ATTTTGCGCTTTTTGCCCAGCATTCAGa                                                                               <  1:3730800/149‑1 (MQ=255)

                                                                        cccTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTGTTG‑TTTTTT‑TTTG‑GTTTTGGGGTTTTTTGCCACCAGTCaaaaaca                                                                           >  1:4584307/1‑143 (MQ=255)

                                                                         ccTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTTTTG‑TTTTGT‑CTGG‑ATTTTGCGCTTTTTGCCCAGCATTCAGA‑CG‑aaa                                                                        <  1:2381395/149‑1 (MQ=255)

                                                                           ggTCGTGGAGCAGGACGCAGCGCGTAAAAACGCATGTGCAACTGTCGGCGCTTTTCCAGTATGTGGCGAAAAATTTTATGAAAAACGGCCTGCGGGCCGTTTTG‑TTTTGT‑CGGG‑ATTTTGCGCTTTTTGCCCAGCATTCAGA‑CG‑AAAAt                                                                      <  2:3440744/148‑1 (MQ=255)

                                                                                 ggAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTTTTG‑TTTTGT‑CTGG‑ATTTTGCGCTTTTTGCCCAGCATTCAGA‑CG‑AAAATTGCCCg                                                                <  2:2243378/149‑1 (MQ=255)

                                                                                  gAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTGTTG‑TTTTTT‑GTGG‑TTTTTTGGGTTTTTTGCCACCAATCATG‑CG‑AAAATTGCCCgg                                                               >  1:3353739/1‑149 (MQ=255)

                                                                                    gCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTGTTG‑GTTTGT‑CTGG‑ATTTTGCGCTTTTTGCCCAGCATTCAGA‑CG‑AAAATTGCCCgggc                                                             >  2:4195553/1‑148 (MQ=255)

                                                                                      aGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTTTTG‑TTTTGT‑CTGG‑ATTTTGCGCTTTTTGCCCAGCATTCAGA‑CG‑AAAATTGCCCGGGAAt                                                           <  2:610009/149‑1 (MQ=255)

                                                                                       ggAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTTTTG‑TTTTGT‑CTGG‑TTTTTTCGCTTTTTGCCCAGCATTCAGA‑CG‑AAAATTGCCCGGGAAtt                                                          >  2:2175951/1‑149 (MQ=255)

                                                                                       ggAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTTTTG‑TTTTGT‑CTGG‑ATTTTGCGCTTTTTGCCCAGCATTCAGA‑CG‑AAAATTGCCCGGGAAtt                                                          <  2:2869336/149‑1 (MQ=255)

                                                                                        gAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTTTTG‑TTTTGT‑CTGG‑ATTTTGCGCTTTTTGCCCAGCATTCAGA‑CG‑AAAATTGCCCGGGAATtg                                                         <  1:1070368/149‑1 (MQ=255)

                                                                                          aGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTTTTG‑TTTTTT‑CTGG‑TTTTTTGGCTTTTTGGCCAACATTCAGA‑CG‑AAAATTTCCCGGGGAttttt                                                       >  2:1843961/1‑146 (MQ=255)

                                                                                          aGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTGTTT‑TTTTTT‑GTGG‑TTTTTGGGGTTTTTGGCCACCAATCAGA‑CGAAAAATTACCCgggcttgtg                                                        >  1:3313055/1‑143 (MQ=255)

                                                                                          aGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTGTTG‑TTTTTT‑GTGG‑GTTTTTGGCTTTTTTGCCCCCAATCAGA‑AC‑CAAAATGCCCCGGACGtgtg                                                       >  1:1365145/1‑149 (MQ=255)

                                                                                          aGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGGTTTG‑TTTTTTGGTGT‑GTTTTG‑GCTTTTTTGCCAGCCATAAGA‑CGAAAAATGGCCCGGGAAtgtg                                                        >  1:3313016/1‑146 (MQ=255)

                                                                                             aTCGCGTAAAAACGCATTTGCAACTGTCGGCGCTGTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTTTTG‑TTTTGT‑CTGG‑ATTTTGCGCTTTTTGCCCAGCATTCAGA‑CG‑AAAATTGCCCGGGAATTGTGaaa                                                    <  2:4608178/149‑1 (MQ=255)

                                                                                                  gTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTTTTG‑TTTTGT‑CTGG‑ATTTTTCGCTTTTTGCCCAGCATTCAGA‑CG‑AAAATTGCCCGGGAATTGTGAAAAAATa                                               >  2:3552937/1‑149 (MQ=255)

                                                                                                  gTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTGTTG‑GTTTGT‑CTGG‑ATTTTGCGCTTTTTGCCCAGCAGTCAGA‑CG‑CAAATTGCCCGGGACATGTGAAAAAATa                                               >  2:4652179/1‑149 (MQ=255)

                                                                                                   tAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTTTTG‑TTTTGT‑CTGG‑ATTTTGCGCTTTTTGCCCAGCATTCAGA‑CG‑AAAATTGCCCGGGAATTGTGAAAAAATAc                                              <  2:3784540/149‑1 (MQ=255)

                                                                                                   tAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTGTTG‑TTTTTT‑GTGG‑GTTTTGGGCTTTTTTACCCCCAAACAGACCA‑AAAATTGCCCGGGAA‑TGTGTAAAAATaa                                              >  1:682401/1‑148 (MQ=255)

                                                                                                   tAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTGTTG‑TTTTGT‑CTGG‑ATTTTGTGCTTTTTGCCCAGCATTCAGA‑CC‑AAAATTTCCCGGGACTTGTGAAAAAATaa                                              >  2:2470210/1‑148 (MQ=255)

                                                                                                   tAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTACTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTGTTG‑TTTTGT‑GTGG‑TGTTTGCGCTTTTTTCCCCGCAATAAGA‑CG‑AAAAATGCCCGGGGAGTGTGAAAAAATaa                                              >  1:2375194/1‑148 (MQ=255)

                                                                                                    aaaaaCGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTGGTGTTTTTTT‑GTGG‑TTTTTTGGTTTTTTTGCCCCCATAAGAA‑GG‑AAAAATGCCCGGGAGTTGTTaaaaataaa                                              >  1:3593714/1‑145 (MQ=255)

                                                                                                    aaaaaCGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGGTGTG‑TTTTTT‑GTGGTTTTTTGGGGTTTTTG‑GCAGCAATCAGA‑CG‑CAAAATTCCCGGGGATTTTGTAAAAATaag                                             >  1:3257658/1‑147 (MQ=255)

                                                                                                      aaaCGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTGTTG‑TTTTTT‑CGGG‑ATTTTGCGGTTTTTTCCCCCCCATCAGA‑CC‑AAAAATGCCCGGGCAATGTGgaaaaaaaaccc                                           >  1:4308194/1‑143 (MQ=255)

                                                                                                       aaCGCATTTGCAACTGTCGGCGCTTTTCCTGTATGTTGATAAAAATTTTATGAAAAACGGCCTGCGGGCCGTTTTG‑TTTTGT‑CTGG‑ATTTTGCGCTTTTTGCCCAGCATTCAGA‑CG‑AAAATTGCCCGGGAATTGTGAAAAAATACGCGa                                          <  1:1723732/149‑1 (MQ=255)

                                                                                                       aaCGCATTTGCAACTGTCGGCGCTTTTCCAGTTTGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTTTTG‑TTTTGT‑CTGG‑ATTTTGCGCTTTTTGCCCAGCATTCAGA‑CG‑AAAATTGCCCGGGAATTGTGAAAAAATACGCGa                                          <  1:1723864/149‑1 (MQ=255)

                                                                                                         cgaaTTTGAAAGTGTCGGCGCTTTTCCTTTTTGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTTTTG‑TTTTGT‑CTGG‑ATTTTGCGCTTTTTGCCCAGCATTCAGA‑CG‑AAAATTGCCCGGGAATTGTGAAAAAATACGCGACa                                        <  1:40083/146‑1 (MQ=255)

                                                                                                         cGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTGTTG‑TTTTGT‑CTGG‑ATTTTGCGGTTTTTGCCCAGCAATCCTA‑CG‑AAAATTGCCCGGGAATTGTGAAAAAATACGCGACa                                        >  2:910322/1‑149 (MQ=255)

                                                                                                           cATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAATCGGCCTGCGGGCCGTTTTG‑TTTTGT‑CTGG‑ATTTTGCGCTTTTTGCCCAGCATTCAGA‑CG‑AAAATTGCCCGGGAATTGTGAAAAAATACGCGACCgc                                      >  2:4320117/1‑149 (MQ=255)

                                                                                                           cATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTTTTG‑TTTTTT‑CTGG‑ATTTTTCGCTTTTTTCCCAGCATTCAGA‑CG‑AAAAATGCCCGGGACTTGTGAAAAAATACGCGACAgc                                      >  2:1929959/1‑149 (MQ=255)

                                                                                                             tttGCAAGTGTCGGGGCTTTTCCTGTTTGTTGCTAAAATTTTTATGAAAAACGGCCTGCGGGCCGTTTTG‑TTTTGT‑CTGG‑ATTTTGCGCTTTTTGCCCAGCATTCAGA‑CG‑AAAATTGCCCGGGATTTGTGAAAAAATACGCGACAgcgc                                    <  1:4216126/149‑1 (MQ=255)

                                                                                                             tttGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTTTTG‑TTTTTT‑CTGG‑ATTTTGCGCTTTTTTACCCGCAATCAGA‑CG‑AAAATTGCCCGGGGATTGTGAAAAAATAAGCGACAgcgc                                    >  2:2138320/1‑149 (MQ=255)

                                                                                                               tGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTTTTG‑TTTTGT‑GTGG‑ATTTTGCGCTTTTTGCCCAGCATTCAGA‑CG‑AAAATTGCCCGGGAATTTTGAAAAAATAAGCGACAgcgcgc                                  >  2:175303/1‑149 (MQ=255)

                                                                                                                 cAACTGTCGGCGCTTTTCCTTTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTTTTG‑TTTTGT‑CTGG‑ATTTTGCGCTTTTTGCCCAGCATTCAGA‑CG‑AAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCaa                                <  1:2852225/149‑1 (MQ=255)

                                                                                                                 cAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTTTTG‑TTTTGT‑CTGG‑ATTTTGCGCTTTTTGCCCAGCATTCAGA‑CG‑AAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCaa                                <  2:2712624/149‑1 (MQ=255)

                                                                                                                 cAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGGTTTT‑TTTTTT‑CTGG‑CTTTTTCGCTTTTTTCCCCGCATTCAGA‑CG‑AAAATTTCCCGGGCATAGTGTAAAAATACGCGAAAGCGCGCac                                >  2:2078476/1‑148 (MQ=255)

                                                                                                                      gTCGGCGCTTTGCCTGTATGTTGCTAAAATTTTTATGAAAAACGGCCTGCGGGCCGTTTTT‑TTTTGT‑CTGG‑ATTTTGCGCTTTTTGCCCAGCATTCAGA‑CG‑AAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAAcc                           <  1:4647864/149‑1 (MQ=255)

                                                                                                                         ggCGCTTTTCCAGTTTGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTTTTG‑TTTTGT‑CTGG‑ATTTTGCGCTTTTTGCCCAGCATTCAGA‑CG‑AAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAACCGtt                        <  2:3908404/149‑1 (MQ=255)

                                                                                                                         ggCGCTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTTTTG‑TTTTGT‑CTGG‑ATTTTGCGCTTTTTGCCCAGCATTCAGA‑CG‑AAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAACCGtt                        <  2:4494641/149‑1 (MQ=255)

                                                                                                                           ggcTTTTCCTTTATGTTTCTAAAAATTTTTTGAAAAACGGCCTGCGGGCCGTTTTG‑TTTTGT‑CTGG‑ATTTTGCGCTTTTTGCCCAGCATTCAGA‑CG‑AAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAACCGTtct                      <  2:4056026/148‑1 (MQ=255)

                                                                                                                           cgcTTTTCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTTTTG‑TTTTGT‑CTGG‑ATTTTGCGCTTTTTGCCCAGCATTCAGA‑CG‑AAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAACCGTtct                      <  2:115810/149‑1 (MQ=255)

                                                                                                                            gcTTTTCCATTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTTTTG‑TTTTGT‑CTGG‑ATTTTGCGCTTTTTGCCCAGCATTCAGA‑CG‑AAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAACCGTtctc                     <  2:740563/149‑1 (MQ=255)

                                                                                                                              ttttCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTTTTT‑TTTTTT‑GTGG‑ATTTTTGGCTTTTTTCCCAGCATTCAGA‑AG‑AAAAATTCCCCGGAAATGTGAAAAAAAAAGCGACACCGCGCCAAAAACCTTCTCGa                   >  1:1126464/1‑149 (MQ=255)

                                                                                                                               tttCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTTTTG‑TTTTGT‑CTGG‑ATTTTGCGCTTTTTGCCCAGCATTCAGA‑CG‑AAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAACCGTTCTCGAc                  <  1:725324/149‑1 (MQ=255)

                                                                                                                                ttCCAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTTTTG‑TTTTGT‑CTGG‑ATTTTGCGCTTTTTGCCCAGCATTCAGA‑CG‑AAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAACCGTTCTCGACt                 <  1:2741930/149‑1 (MQ=255)

                                                                                                                                  ccAGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTTTTG‑TTTTGT‑CTGG‑ATTTTGCGCTTTTTGCCCAGCATTCAGA‑AG‑AAAATTTCCCGGGAATTGTGAAAAAATAAGCGACAGCGCGCAATAACCGTTCTAGACTCa               >  2:3951256/1‑149 (MQ=255)

                                                                                                                                    aGTATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTTTTG‑TTTTGT‑CTGG‑ATTTTGCGCTTTTTGCCCAGCATTCAGA‑CG‑AAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAACCGTTCTCGACTCATa             <  1:165480/149‑1 (MQ=255)

                                                                                                                                      tATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTTTTG‑TTTTTT‑CTGG‑ATTTTTCGCTTTTTTCCCCGCATTCATA‑AG‑AAAATTTGCCGGGCATTGTGAAAAAATACCCGCCAGCGCGCAAAAACCGTTGTAGACACCTaaa           >  2:3725661/1‑149 (MQ=255)

                                                                                                                                      tATGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTGTTG‑TTTTTT‑CTTG‑ATTTTTCGCTTTTTTCCCCCCATTAATA‑CG‑AAAATTGCCCCGGAAATTTGAAAAAATACGCGACAGCCCGCCAAAACCCATCTCGACACAaaaa           >  2:4445905/1‑149 (MQ=255)

                                                                                                                                       aTGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTTTTG‑TTTTGT‑CTGG‑ATTTTGCGCTTTTTGCCCAGCATTCAGA‑CG‑AAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAACCGTTCTCGACTCATaaaa          <  2:3369623/149‑1 (MQ=255)

                                                                                                                                       aTGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTTTTG‑TTTTGT‑CTGG‑ATTTTGCGCTTTTTGCCCAGCAATCAGA‑CG‑AAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAAAAACCGTTCTCGACACATaaaa          >  2:621837/1‑149 (MQ=255)

                                                                                                                                        tGTTGCTAAAAATTTTATGAAAAACGGCCTGCGGGCCGTTTTG‑TTTTGT‑CTGG‑ATTTTGCGCATTTTGCCCCCCAATAACA‑CG‑AAAAATGCCCCGGCAGTTTTAAAAAATACGCGACAGCCCGCAAAAACCCTTCTCGACTCATAAAAg         >  2:1368192/1‑149 (MQ=255)

                                                                                                                                               aaaaATTTTATGAAAAACGGCCTGCGGGCCGTTTTG‑TTTTGT‑CTGG‑ATTTTGCGCTTTTTGCCCAGCATTCAGA‑CG‑AAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAACCGTTCTCGACTCATAAAAGTGATGcc  <  2:2600583/149‑1 (MQ=255)

CATGCGCTGGCTCTGCATACTTACACTATTAAGGAGTGGGAAGGGTTGCAGGACACCGTCTTTGCCTCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTG‑TTTTGT‑CTGG‑ATTTTGCGCTTTTTGCCCAGCATTCAGA‑CG‑AAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAACCGTTCTCGACTCATAAAAGTGATGCC  >  CP009273/450144‑450435

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